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Abstract 



Several aspects of classical and quantum mechanics applied to a class of strongly chaotic 
systems are studied. The latter consists of single particles moving without external forces 
on surfaces of constant negative Gaussian curvature whose corresponding fundamental 
groups are supplied with an arithmetic structure. 

It is shown that the arithmetical features of the considered systems lead to exceptional 
properties of the corresponding spectra of lengths of closed geodesies (periodic orbits). 
The most significant one is an exponential growth of degeneracies in these geodesic length 
spectra. Furthermore, the arithmetical systems are distinguished by a structure that ap- 
pears as a generalization of geometric symmetries. These pscudosymmetries occur in the 
quantization of the classical arithmetic systems as Hecke operators, which form an infinite 
algebra of self-adjoint operators commuting with the Hamiltonian. 

The statistical properties of quantum energies in the arithmetical systems have previ- 
ously been identified as exceptional. They do not fit into the general scheme of random 
matrix theory. It is shown with the help of a simplified model for the spectral form fac- 
tor how the spectral statistics in arithmetical quantum chaos can be understood by the 
properties of the corresponding classical geodesic length spectra. A decisive role is played 
by the exponentially increasing multiplicities of lengths. The model developed for the 
level spacings distribution and for the number variance is compared to the corresponding 
quantities obtained from quantum energies for a specific arithmetical system. 

Finally, the convergence properties of a representation for the Selberg zeta function 
as a Dirichlet series are studied. It turns out that the exceptional classical and quantum 
mechanical properties shared by the arithmetical systems prohibit a convergence of this 
important function in the physically interesting domain. 
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1 Introduction 



Two classes of theoretical activities in the field of natural sciences may be distinguished. 
The first one consists of efforts towards unraveling the kinematics of a given system and of 
finding the underlying dynamics governing the time evolution of the system; the second one 
is formed by the development of methods and techniques to obtain numerical data from 
the theory describing the dynamics in order to be able to compare them with empirical 
data. A sensible and satisfying theory of natural phenomena then has to meet these two 
requirements: it should be logically consistent and ought produce testable predictions. 
Both of these aspects have to be well distinguished. One consequence of this consideration 
is that a theory containing deterministic dynamics need not be predictive with respect to 
its time evolution. This realization lies at the heart of what nowadays is widely known as 
deterministic chaos. 

To be more specific we now consider Hamiltonian dynamical systems (with a finite 
number N of degrees of freedom) in classical mechanics. These may be described on N- 
dimensional configuration manifolds M with Riemannian metrics ds'^ = gijdq^dq^ defined 
on them, (g = (g^, . . . ,q^) are local coordinates on M.) The dynamics is specified by 
providing M with a Lagrangian function L(q,q) = |(^)^ — V{(]). The equations of mo- 
tion for this system may then be obtained by Hamilton's principle as the Euler-Lagrange 
equations 

d dL dL , , 

Specifying initial values go and g'g at a time to then uniquely fixes the time evolution q{t) 
of the system for all times t > to. This is the manifestation of determinism in classical 
mechanics. As yet, nothing is said, however, about predictability. Since in practice it is 
never possible to prepare a system at an initial time to to be in a definite state (go, go), one 

• — * • — * 

has to allow the initial values to be taken from [% — e,% + e\ x [% — 3,% + S] for some 

— * 

small uncertainties e and 6. Predictability would now require the uncertainties to grow 
only modestly under the time evolution dictated by (1.1). By this we mean an increase of 
l^l and \6\ at most like a power of t. It is, however, by no means clear from the information 
provided so far that this will be the case. 

The class of dynamical systems possessing the most regular kind of time evolution is 
given by the integrable ones. For them there exist N independent constants of motion with 
pairwise vanishing Poisson brackets. If one considers their time evolution in phase space 
(i.e. on the cotangent bundle T*M), this is found to take place on an iV-dimensional torus. 
The equations of motion can then be integrated by quadratures. Integrable systems are 
predictable in the sense just introduced. 

The other extreme is given by irregular systems sharing the property of ergodicity. 
Ahiiost all of their phase space trajectories fill the (2A^ — 1) -dimensional hypersurface of 
constant energy in phase space densely. The probability of finding an arbitrary phase 
point in some bounded region of the hypersurface of constant energy is proportional to the 
volume of that region. This means that the trajectories of ergodic systems are uniformly 
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distributed in phase space. Among such systems one can find a hierarchy of even higher 

irregularities: mixing systems, Anosov-systems, K-systems, . . .; see e.g. [37] as a reference. 

In between the two extremes of integrable and ergodic systems there exists any kind of 
intermediate behaviour. These systems, however, wiU not be deah with in the course of 
the present investigation. 

It now appears to be useful to reformulate the setting a little bit. Suppose the system 
under study is in a state of energy E. According to the conservation of energy it can 
visit only those parts of M during its time evolution, where V{q) < E is fulfilled. On 
this domain of M one introduces the Jacobi metric dS"^ := [E — V{q)]ds'^. Maupertuis' 
principle of least action, being equivalent to the equations of motion (1.1), is now also 
equivalent to the statement that the trajectories of the system on M are geodesies in the 
Jacobi metric dS"^, see e.g. [1] for details. Hence every Hamiltonian dynamics can be 
viewed as the geodesic fiow on some Riemannian manifold. Hopf [56] has shown that a 
negative curvature associated with the Jacobi metric is a sufficient condition to render the 
system ergodic. 

A guiding principle in the present study will be to keep things as simple as possible, 
without giving up the essential structures that determine the properties of a -hopefully- 
general enough class of systems. Hence we agree to restrict our attention to the following 
kind of examples. They will have two degrees of freedom, since this is the minimal dimen- 
sionality required for a Hamiltonian system to be non-integrable; the reason for this being 
that conservation of energy renders every system of one degree of freedom integrable. Hav- 
ing in mind the above remark on the Jacobi metric and its curvature, we choose geodesic 
flows on two dimensional Riemannian surfaces of constant negative Gaussian curvature as 
our prime examples. We are then going to study the interplay between the classical and 
quantum dynamics for these systems. 

A central issue of this investigation will be to identify fingerprints of the classical 
properties of a given dynamical system in its quantum version. Since the time evolution 
in (non-relativistic) quantum mechanics, which is governed by the Schrodinger equation 



for the time dependent wave function ip{t,q), is linear, no chaotic phenomena like an 
exponential sensibility to initial conditions can occur. Thus the question arises, how a 
potential irregularity of the classical system can show up in the semiclassical limit ^ ^ 0. 
Although this limit is not "smooth" in that at the value h = typical quantum mechanical 
structures and quantities no longer exist, it is expected that for small values of h one should 
be able to detect classical properties of the system. One such sign of the structure of the 
classical phase space in quantum mechanics is the applicability of semiclassical quantization 
rules such as the WKB-method (for one-dimensional systems) or the EBK-method (for 
multi-dimensional systems), which is restricted to the integrable case. Therefore a first and 
obvious question would be to find a semiclassical quantization procedure e.g. for strongly 
chaotic classical systems. 
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An answer to this question is provided by Gutzwiller's periodic- orbit theory [42, 44], 
in which the spectrum of the quantum Hamihonian H is determined in a semiclassical 
approximation by the set of actions evaluated along the classical periodic trajectories {pe- 
riodic orbits) of the system. Plane billiards seem to be well suited for an explanation of 
the general procedure. Let therefore D C be some connected domain on the euclidean 
plane. Its boundary dD is assumed to be piecewise smooth. A plane billiard then consists 
of a particle moving freely inside D and being elastically reflected on dD, which will be 
chosen such that the classical dynamics is strongly chaotic. This should mean that the sys- 
tem is ergodic and that all periodic orbits are isolated in phase space and are unstable. The 
action along a periodic orbit 7 is then given by = pl^, where p denotes the particle's 
momentum and 1^ is the length of the orbit 7. The quantum mechanical Hilbert space for 
such a billiard is given by the space of square integrable functions ijj{x, y) on D that vanish 
on dD. In units where h = 1 = 2m the quantum Hamiltonian is H = —Ae = — (<9^ + dy) 
and possesses a purely discrete spectrum < Ei < E2 < . . ., = p^. Weyl's famous 
law for the number N{E) of energy eigenvalues not exceeding E states that asymptotically 
for £■ — > 00 

N{E)r~.^^E. (1.3) 

Thus En ~ Q^YQSiii^D) ^ — ^ 00, and hence the resolvent operator {H — E) ^ is not of trace 
class. Introducing, however, a suitable smearing allows to investigate trace class operator 
valued functions of H. In [88] one can find a regularized version of Gutzwiller's trace 
formula. Asymptotically in the semiclassical limit it reads 



In this formula h{p) is an even function, holomorphic in the strip \Imp\ <r — l + e, £>0, 
that decreases faster than for \p\ 00. g{x) := ^e^^^h{p), and d{p) is a mean 



spectral density expressed as a function of momentum p = y E. Under these assumptions 
the integral and the sums involved in (1.4) are absolutely convergent, x-y = (~1)''^) where 
is the number of reflections on dD when traversing 7 once. A-y = u^/l^ is the Lyapunov 
exponent of 7 and is its stability exponent, denotes the sign of the trace of the 
monodromy matrix for 7. The topological entropy t describes the proliferation of the 
number Mil) of periodic orbits of lengths not exceeding / by 



pTl 



N{1) ~ — , / ^ 00 . (1.5) 

tI 

A denotes the asymptotic average of the A^'s and is also known as the metric entropy. 

The trace formula (1.4) can be viewed as an identity in the sense of distributions 
for the spectral density d{E) = J2'^=i S{E — En) (expressed by the variable p) of the form 
d(p) ~ d(p)-\-dfi{p). dfi{p) is an oscillatory correction to the mean density d(p), determined 
by the periodic-orbit sum in (1.4). The conditions on h{p) stated after (1.4) then define 
the space of test functions. 
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It is possible to use Gutzwiller's trace formula in order to derive semiclassical quan- 
tization rules that yield approximations for the quantum energies {E^} in terms of the 
actions of periodic orbits in the corresponding classical system, see e.g. [44] and the ref- 
erences therein. In some sense the so-obtained quantization rules are a substitute for the 
EBK-quantization of classically integrable systems. A fundamental difference between the 
semiclassical quantization schemes for classically integrable and chaotic systems is provided 
by the computational effort to be spent in order to resolve some quantum energy E^. In 
the integrable case this is independent of the energy, whereas (1.5) requires for chaotic sys- 
tems an exponentially increasing number of periodic orbits to be taken into account when 
trying to compute higher and higher energies. This observation reflects the high degree of 
irregularity and complexity of chaotic systems. 

Historically, the first chaotic systems under investigation were of mathematical nature 
and mainly served as examples to develop the mathematical theory of dynamical systems. 
Hadamard [46] studied the geodesic flow on closed surfaces endowed with Riemannian met- 
rics of constant negative Gaussian curvature, which also go under the notion of hyperbolic 
surfaces, and discovered an instability of the trajectories of such flows. Since the surfaces 
he studied cannot be realized as being embedded in M^, they were considered as purely 
mathematical examples. For a specific hyperbolic surface Artin could later prove [2] the 
ergodicity of the geodesic flow on it. Moreover, at this occasion he introduced a symbolic 
dynamics for this system. This was the first time that the property of ergodicity could ever 
be rigorously demonstrated for a dynamical system. One could thus view Artin's paper as 
the foundation of the modern theory of dynamical systems. Only after the development of 
powerful computers that allow for a quantitative analysis of irregular dynamical systems 
these subjects arose interest among a considerable number of physicists. Starting with 
the series of papers [42] by Gutzwiller a discussion of the quantization of chaotic classical 
systems was rendered possible. Again it was Gutzwiller who noticed [43] that his trace 
formula turned into an exact identity when applied to geodesic fiows on surfaces of con- 
stant negative curvature. In mathematics this identity had long before been introduced 
by Selberg [84], who intended to understand the Riemann zeta function and the Riemann 
hypothesis with the help of his trace formula. 

Since Gutzwiller's observation hyperbolic surfaces have been intensively studied in the 
context of the quantization of classically chaotic systems. The first numerical results on 
quantum energies for a hyperbolic triangle obtained from a solution of the stationary 
Schrodinger equation by Schmit have been presented in [31, 17]. Aurich and Steiner [10] 
first calculated lengths of periodic orbits and also quantum energies [9, 11] on a specific 
hyperbolic surface, the so-called regular octagon. The lengths were used as an input to 
evaluate the Selberg trace formula [9, 11]. By chance, the regular octagon that had been 
chosen for the numerical studies turned out to be very specific: its fundamental group is 
of an arithmetical origin. It was realized [10, 4] that the arithmetic structure underlying 
the regular octagon leads to exponentially increasing multiplicities of lengths of periodic 
orbits. As was later noticed [12, 13], this property is exceptional and has remarkable 
consequences. One of these is the occurrence of unexpected statistical properties of the 
related quantum energy spectra. Although the classical system is as chaotic as it could be. 
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the quantum spectral statistics are more alike those of classically integrable systems than 
those expected for generic classically chaotic ones. Namely, due to a conjecture of Bohigas, 
Giannoni and Schmit [30] which was mainly based on numerical observations, generic 
chaotic systems excel by quantum energy spectra that can be well described by eigenvalues 
of large random matrices. The resulting spectral statistics then differ drastically from 
those in the integrable case as obtained before by Berry and Tabor [24]. Somewhat later 
it was shown [12], however, that the energy fluctuations of 30 non-arithmetic octagons arc 
in agreement with the predictions of random matrix theory. At the time the unexpected 
statistical properties of the eigenvalues for the regular octagon were first observed it was, 
however, not clear that these follow from the arithmetic properties of the fundamental 
group. This was later clearly expressed to hold for all arithmetical systems in the two 
simultaneously appearing papers [27] and [32]. 

The objective of the present investigation now is to study the class of dynamical systems 
arising from geodesic flows on hyperbolic surfaces with arithmetic fundamental groups both 
from a classical as well as from a quantum mechanical point of view. The ultimate goal 
then will be to understand the exceptional statistics of the quantum energy spectra. The 
problem of studying the wave functions arising on the quantum mechanical side, however, 
will not be treated here. It seems that the cigcnfunctions of the arithmetical systems are 
not as exceptional as the eigenvalue spectra, see [54, 79, 16]. 

The organization of this paper is as follows. Chapter 2 reviews some general properties 
of (discrete) quantum energy spectra and introduces the necessary means and notions for 
their investigation in order to prepare for studying the spectral statistics of arithmetical 
systems. Since in order to be ready for the latter the classical aspect of the problem has 
to be understood in quite some detail, chapter 3 is devoted to an investigation of geodesic 
length spectra on arithmetic hyperbolic surfaces. A central result, obtained in section 
3.4, is the observation of exponentially growing degeneracies in arithmetic length spectra. 
Another peculiarity arising from the arithmcticity of a fundamental group is the existence 
of infinitely many pseudo symmetries on the corresponding surface. It will be attempted 
to give a geometric picture of this structure in section 3.5. Chapter 3 will be closed by an 
investigation of fluctuations of the lengths of closed geodesies on both arithmetic and non- 
arithmetic hyperbolic surfaces. These fluctuations are closely related to the fluctuation 
properties of the corresponding quantum energies. Using the latter in order to learn about 
the former is also known as inverse quantum chaology. 

Chapter 4 then deals with quantum energy spectra of arithmetical systems. At first 
Hecke operators are discussed as representations of pseudosymmetries on the wave functions 
and it is found that this structure results in constraints on the eigenvalue spectra. This 
realization is taken as a first hint towards exceptional spectral statistics. Thereafter the 
spectral form factor is reviewed and its properties in the arithmetical case, derived from 
the exponentially increasing multiplicities of lengths of periodic orbits, are employed to 
introduce a simplified model form factor. The latter can then be used to derive a model 
for the level spacings distribution as well as for the number variance in arithmetical chaos. 
The model will also be compared to existing numerical data in order to test its quality. 
Finally, convergence properties of the Selberg zeta function [84], which arises from the 
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Selberg trace formula and plays an important role for obtaining quantization rules, are 
studied. It is observed that the spectral statistics influence the convergence properties of 
this function considerably, and that the arithmetic systems suffer from a lack of convergence 
in the physically interesting domain. 

Chapter 5 finally summarizes the findings of the preceding investigations. Furthermore, 
three appendices are included. The first one reviews the theory of the Riemann zeta 
function, which is in many respects similar to the Selberg zeta function. The analogy of 
these two functions is particularly useful for inverse quantum chaology. Appendix B treats 
the desymmetrization of a class of hyperbolic surfaces explicitly in order to illustrate the 
general procedure of getting rid of unwanted symmetries on hyperbolic surfaces. The final 
appendix collects the definitions of 0-, o-, and Jl-estimates for remainder terms that occur 
at several instants in the main body of the text. 

2 Some General Remarks on Discrete Quantum En- 
ergy Spectra 

Let a classical dynamical system be given with a finite number of degrees of freedom. 
After quantization its Hamiltonian H should have a discrete spectrum Q < Ei < E2 < 
E2 . . .. A quantization rule then is the specification of a function f : IN ^ IR such that 
En — /(^)- The exphcit knowledge of such a quantization function completely solves the 
problem of determining the energy spectrum of the system. If one is merely interested in 
energy eigenvalues and not in wavefunctions, a quantization rule is equivalent to solving 
the stationary Schrodinger equation Hifjn = Enipn- One possibility to obtain a quantization 
function is to study the spectral staircase 

N{E) := ^ 1 . (2.1) 

E„<E 

Defining then No(E) :— | lim£_^o [N{E + e) + N{E — e)] yields the quantization rule 

No{En) ^n-^, nelN , (2.2) 

which can be converted to the condition cos(7riVo(ii^)) = for E to be an energy eigenvalue, 
see [14, 7]. This consideration stresses the importance of studying the spectral staircase 
thoroughly if one is interested in the quantization of a dynamical system. 

There is, however, only a very restricted number of examples where an exact quan- 
tization function is explicitly known, among which are the typical textbook examples of 
a particle in a box, the harmonic oscillator, or the Coulomb potential. Semiclassically, 
quantization rules are provided by the WKB- or EBK-methods for all classically inte- 
grable systems. For these there exist canonical transformations of the classical phase space 
variables (p, q) to action-angle variables (/, Cu), such that the Hamiltonian function H(p, q) 
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transforms to one being only dependent on the actions /, H — H{l). Then 



Erh^H{I^), I^^m+^a, (2.3) 

is the desired (semiclassical) quantization rule. In (2.3) fn runs through and a is the 
vector of Maslov indices. 

It is one of the major goals of quantum chaology to obtain an analogue of (2.3) for 
strongly chaotic systems. Since there do not exist action-angle variables and even no rem- 
nants thereof for crgodic systems, the quantization problem has to be tackled in a totally 
different manner. A starting point one could think of would be to study the spectral stair- 
case in as much detail as possible, thereby keeping in mind the relation (2.2). For such an 
analysis it proves useful to split N{E) into a smooth part N{E), which is an approximation 
to N{E) in that the steps have been smeared out, and a remaining contribution Nfi(E). 
This then describes the fluctuations of N{E) about its mean N{E). A priori there is no 
preferred way to define the smoothing N{E). The splitting 

N{E) = N{E) + Nfi (E) (2.4) 

should only be subject to some general requirements. N{E) should be smooth without 
"too many" oscillations (ideally one would require it to be monotonically increasing), and 
it should be asymptotically identical to the true staircase, i.e. N{E) N{E)j for E oo. 
Nfi{E) then should describe fluctuations about N{E), i.e. Nfi{E) = N{E)-N{E) should 
fluctuate about zero in the limit E ^ oo. This requirement can be expressed in a variety 
of different ways, one of which is 

lim Y NfiiEr,) = , (2.5) 

i.e. an asymptotic vanishing of an arithmetic mean of Nfi{E). Another, although similar, 
requirement uses an integral instead of a sum, namely 

lim ^ dE Nfi{E) = . (2.6) 

L-^oo L Jo 

Provided there is a Gutzwiller trace formula available, like (1.4) for chaotic plane bil- 
hards, a natural sphtting (2.4) is yielded by the two contributions to the r.h.s. of the trace 
formula. In order to obtain this one introduces the complex variable s — ^ — ip, i.e. 

E{s) = = — (s — l)^. A regularized trace of the resolvent operator may be introduced 
by choosing the function h{j)) = ^^^|_^2 — e{J)-p'^ ' ^'^^ Res, Re a > r, which fulflUs the 
requirements for a test function to be used in (1.4). Inserting this into the trace formula 
yields 



E 

n=l 



1 



E{s) - E^ E{a) - E^ 



poo 

2 / dp p dip) 
Jo 



E{s)-p^ E{a)-p^ 
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EE: 

(2.7) 



1 x^/^ e-(''-t)«7 



tEE 



Choosing now -^((t) e iR, i.e. Re a and keeping Res > ^, i.e. ImE{s) > 0, 

taking the imaginary part on both sides of (2.7), and multiplying the result with — - yields 

(E=p^p>0) 

£ 6{E - E^) ~ -d{E) - Im (^^^1^1 , (2.8) 
^1 { 2«P7r ^(1 _ ip) J 

where the dynamical zeta function Z{s) has been introduced. It is for Res > t defined by 
the Euler product 

oo _ 

Z{s) n n (l - X7(7!;e-(^+"^^+^(^--^))'^) . (2.9) 

7 71=0 

Integrating both sides of (2.8) from zero to E yields (assuming argZ(|) = 0) 

N{E) ^N{E) + ^axgZ (^+i\f^ , (2.10) 

where N{E) = dE'd{E'). One has now obtained, in the semiclassical limit E ^ oo, 
that _ 

Nfi{E)^^a.rgz(^^ + i^^ . (2.11) 

This relation stresses the importance of dynamical zeta functions for the study of quantum 
spectral properties of chaotic dynamical systems. In addition, dynamical zeta functions 
turn out to be powerful tools to obtain quantization rules in periodic-orbit theory directly 
[89, 68, 95, 26, 63, 23, 13]. Their zeroes s„ on the critical line Res — ^ yield semiclassical 
approximations to the quantum energies via En — E{sn)- This result may be obtained by 
integrating both sides of (2.7). 

In [11, 14, 7] a periodic-orbit expression for the r.h.s. of (2.11) has been obtained 
that enabled the authors to compute Nfi{E) numerically in some approximation from the 
lengths of primitive periodic orbits for several chaotic systems. To do so one inserts the 
Gaussian h{p) = -^(e-^P-i)V^^ + e-^P+i^y^^), with g{x) = ^cos{qx) e-^'^'/^ into the trace 
formula (1.4). Then one integrates the resulting formula in q from zero to p, p^ = E. Since 
lim^^o h{p) = S{p — q) + 5{p + q), the l.h.s. of (1.4) yields for £ > a smoothing N^{E) of 
the spectral staircase N{E) with lim^^o -^£(-5') = N{E), 

limE7^/>^ {e-'^ + e-^) = tQ{E-En) = NiE). (2.12) 



^1 e^/^ Jo 
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On the r.h.s. of (1.4) one obtains Ns{E) + N^ji{E), with 

- - E Ej ,.^.7. _ M^kQ . (2.13) 

As long as £ > 0, the double sum on the r.h.s. of (2.13) converges absolutely, since then 
the Gaussian damping factor ^ 't/^ overcompensates the exponential proliferation (1.5) 
of periodic orbits. 

It will now be shown, using the smoothing N^ji{E) for the fluctuating part Nfi{E) of 
the spectral staircase, that the two conditions (2.5) and (2.6) are satisfied by the sphtting 
(2.10). Thus, in the semiclassical hmit, Nfi{E) indeed describes the fluctuations of N{E) — 
N{E) about zero, which is a major justification to identify N[E) as a mean spectral 
staircase. The essential part of demonstrating that the requirement (2.5) is met is to 
study, after inserting (2.13) into (2.5), 

J2 Mkl.Pn) (2.14) 



Nil) ^ 

for every term in the summation over periodic orbits 7. For the following the reasonable 
assumption has to be made that the momenta pi, . . . ,pn are linearly independent over Z 
for any N G IN. It is then known [61] (see also Bohigas' contribution in [41]) that the 
functions cospit, . . . , cospArt are statistically independent. If A is any Lebesgue-measurable 
set in M and 1{A) is its Lcbcsgue-measure, define the relative measure of A as lr{A) :— 
limT^oo n (-T, +T)}. If now 

[2" ^ 

XN{t) '■= E cospn^ , (2.15) 



n=l 



and AN{a,b) := {t; XN{t) e (a, b)} denotes the set of the t & M such that XN{t) takes its 
values in the interval (a, b), then [61] 



lim lr{AN{a,b)} = ^= dx e"^^" . (2.16) 



In other words, the values of the function x^it) approach a normal distribution with zero 
mean and unit variance in the limit N 00. If one then replaces cosp„t by sinp^t (for 
t = kl^), which does not change the above argument, the values of the function 

f^{kl^):^^= ^ sm{kl,pr,) (2.17) 

become, in the limit L — > 00, Gaussian distributed with zero mean and variance o-j = ^. 
Inserting (2.13) into (2.5) leads to a triple summation. The one over the energies En yields 
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(2.14). The remaining ones over periodic orbits 7 and their repetitions k sum /l(^) taken 
at t — klj times a prefactor, 

E N.AEn) = ^— E E t .mA ^k.-ku..,M kl,) . (2.18) 



Since fiikl-y) to be summed over is Gaussian distributed about zero, this is for kl^ — > 00 
an effectively bounded function. Therefore the sum on the r.h.s. of (2.18) converges due 
to the regularizing Gaussian damping as long as £ > 0. In the limit L — > 00 then the 
complete expression (2.18) vanishes like '^y^^- Weyl's law (1.3) thus determines the rate 

of vanishing to be of the order of 

In order to proof the validity of (2.6), when N^ji{E) from (2.13) is used instead of 
Nfi{E), one has to calculate 



— / dE sin(kLvE) — — dp p sin(pA;L) 
L Jo L Jo 



2 2 

cos{kl^\fL) + -r-r^ sm{kl^VL) , (2.19) 



kl^^/L ' ^ ' {kl^YL 

which vanishes for L ^ 00. The double sum over periodic orbits and their repetitions one 
is in analogy to (2.18) left with therefore also vanishes. It is interesting to notice that both 
mean values (2.5) and (2.6) extended over a finite interval of length L vanish as for 
L — > 00. For every e > the conditions (2.5) and (2.6) are thus fulfilled, and hence also 
in the limit e 0. It seems that any other reasonable requirement expressing the same 
statement on the fluctuations of N^ji{E) will also be satisfied, so that the splitting (2.10) 
seems to be a natural one for a semiclassical analysis of the spectral staircase N{E). This 
point may stress the importance and usefulness of employing Gutzwiller's trace formula 
and dynamical zeta functions. 

It should be remarked that completely analogous results can be obtained for the spectral 
density d{E) = ^ N{E). Splitting d{E) = d{E) + dfi{E) according to (2.8) and using the 
regularization (2.13) yields 

1 ^ ~ L e-='^'''/4 



4 fliE) = V E „ " . ^ /o COs{\fEkL) . (2.20) 



Again, the two averaging procedures (2.5) and (2.6) can be employed, leading to the result 
that deji{E) vanishes with respect to those mean values if one extends the interval (0, L) 
to be averaged over to infinity, i.e. for L — > oo. Compared to N^ji{E) these means of 
deji{E) vanish faster, namely like . 

As described above a major objective in quantum chaology is to explain and classify the 
statistical properties of quantum energy spectra for systems with a chaotic classical limit, 
and one way to achieve this is to describe the spectral staircase as explicitly as possible. 
Since N{E) is smooth and rather "harmless", the more interesting object is Nfi{E). To 
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be able to compare systems with different mean behaviours of their respective spectral 
staircases one conventionally introduces an unfolding of spectra. Define x := N{E), Xn = 
N{En). The unfolded spectrum has a spectral staircase that shall be denoted by 

N{x). By its very definition, N{x) = x and thus 

N{x) ^ X + Nfi{x) . (2.21) 

After being unfolded different spectra can be compared quite easily. They then only differ 
in the fluctuating parts of their respective spectral staircases. 

A conventional measure to investigate spectral statistics is the spectral rigidity ^.^{L; x), 
originally introduced by Dyson and Mchta [38] to analyze spectra of complicated atomic 
nuclei. It is defined as the integrated quadratic deviation of the spectral staircase N{x) 
from the best fitting straight line over an interval [x — L,x + L], 

As{L; x) :=< min ^ f^"^ dy [N{x + y) - A - Byf > . (2.22) 

{A,B) ZL J-L 

< . . . > denotes an averaging in x over an interval of length Ax, such that Ax <S x and 
Ax ^ 1. The latter condition means that the averaging should take place over many 
eigenvalues. (The mean distance of neighbouring eigenvalues is one due to N{x) = x.) At 
the same time the interval should be much smaller than x itself. This procedure presupposes 
that X ^ 1, i.e. it is only possible in the semiclassical regime, and is thus also referred to 
as a semiclassical averaging [20]. Notice that (2.22) differs from the definition in [20] by 
the replacement L — 2L. 

Introducing the splitting (2.21) in (2.22) one obtains the minimum on the r.h.s. for 

1 

1 + ^ dy y Nfi{x + y) . (2.23) 
We will soon see that A^in x and i^min ^ 1 in the limit L — > oo. Thus one finds that 
Aoo(x) := A3(oo; x) = lim < ^ j^"^ dy [N^i^x + y)f > , (2.24) 

L-^oo ZL J-L 

which means that Aoo(a;) is the quadratic mean of Nfi{x). This expression will also some- 
times be denoted as < Nfi{xY >. 

When considering the limit L ^ oo one first encounters the problem to define Nfi{x) 
for X < 0, a situation occurring in (2.22) for L > x. The only reasonable definition seems 
to be Nfi{x) = for a; < 0. Since N{E) = ioi E < 0, an obvious choice for N{E) seems 
to be N{E) = for £; < 0. Then also Nfi{E) = for £; < 0. This is consistent with the 
fact that s = I - i\^ e IR for E < 0, and thus argZ(| + iVE) = 0. We therefore cut 
off, in the limit L — > oo, the integrations on the r.h.s. of (2.23) at —x. 



A 
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As our whole discussion deals with two dimensional problems, Weyl's law (see e.g. (1.3)) 
yields N{E) ~ const. E, E ^ oo. Thus d{E) ^ d = const. , £■ — > oo. We are interested in 
the semiclassical limit and hence use x — N{E) dE from now on. Thus 



A possible way to treat the integrals on the r.h.s. of (2.25) is to use the smoothing (2.12). 
This means that Nsji{E) for e > from (2.13) has to be inserted into (2.25), and the limit 
e — > has to be taken after the evaluation of the resulting integrals. The vanishing of the 
r.h.s. of the first line of (2.25) for L ^ oo has already been demonstrated, see (2.19). The 
second fine can be treated similarly, leading to 



which also vanishes for L — > oo. Thus for any e > the smoothings of (2.25) vanish in the 
limit L — s> oo. From this we conclude that also after performing e — the integrals on the 
r.h.s. of (2.23) vanish in the limit L —>■ oo, leading to — > x and -B^m ~^ 1) as it has 
been claimed above. 

For finite values of L the spectral rigidity measures deviations of the spectral staircase 
from the best fitting straight line over an interval of length 2L. Therefore A3(L; x) indicates 
correlations in a quantum energy spectrum on a scale L. For L the fact that N{x) 
is a step function causes the rigidity to be j^L. For slightly larger values of L A^{L;x) 
begins to measure spectral correlations. A completely uncorrelated spectrum excels by a 
rigidity of A3(L;x) = j^L throughout the whole range of L-values, < L < oo. Such 
a spectrum is obtained by a Poisson process through Xj+i = Xi + Sj, where the Sj's arc 
independent outcomes of measurements of the random variable s, which is distributed 
according to the probability density P{s) = e^^. The Sj called the nearest- 

neighbour level spacings, or briefly level spacings, of the unfolded spectrum {xi}. Due to 
its universal behaviour for L — > the rigidity is, however, not the appropriate tool to 
investigate short-range correlations in a spectrum. It becomes useful only for medium- and 
long-range correlations corresponding to L ^ 1, i.e. on scales taking several levels into 
account. The level spacings are then suitable quantities to measure correlations on small 
scales L ^. 

The systems to be considered from now on shall be completely desymmetrized in order 
to avoid degeneracies in their spectra that superimpose the effects one is looking for. For 
such an arbitrary unfolded quantum energy spectrum one can then study the distribution 




(2.25) 
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of the level spacings of the Xn < x and ask the question, whether there will exist a limiting 
distribution P{s), 

hm #K x^<x s eia,b)} ^ 

=^-^0° N{x) Ja 

and what it will look like if it exists. The question for the existence of P{s) is a very 
subtle one and no mathematically rigorous answer is known in general. For some specific 
examples, where one does have a rigorous theory, see [90, 25]. Extensive numerical calcula- 
tions in numerous examples, however, indicate that the existence of a limiting distribution 
cannot be doubted. Therefore this question will not be pursued any further. 

It turns out that answering the question for the precise form of P{s) allows for a classi- 
fication of quantum systems according to the properties of their respective classical limits. 
Going further on one can also study medium- and long-range correlations in the spectra and 
try to extend this classification to the spectral rigidity As{L; x). Berry developed in [20] a 
semiclassical theory of the spectral rigidity based on periodic-orbit theory (see also Berry's 
contribution in [41]). He found a universal behaviour of A:^{L;x) for 1 ^ L ^ Lmax] 
Lmax is proportional to 7fr—, where T^i„ is the period of the shortest periodic orbit of 

^ min 

the classical system. (For biUiard systems Tj^in — so that Lj^ax ^ V^-) In the 

range 1 <^ L <^ L^nax classically integrable systems share a rigidity of ^^{L^x) = ^L, 
whereas for classically chaotic systems possessing a time-reversal invariance Berry found 
that A3(L; x) = logL + const.. For very large values of L, L ^ Lj^ax, he obtained for 
both cases a saturation of the rigidity at a value that is determined by the contributions of 
short periodic orbits to dfi{x). The energy dependence of the saturation value A^oix) now 
is again characteristic for the classical limit of the quantum system. If this is integrable. 
Berry can show that A^{x) ~ const, a/x, x ^ oo. Chaotic time-reversal invariant systems, 
however, are demonstrated to yield Aoo{x) ~ ^ logx, x — > oo. According to (2.24) Aoo{x) 
is just the mean square < Nfi{xY > of the fluctuation part of the spectral staircase. Hence 
Berry found, in the semiclassical limit, a universal behaviour for this quantity that only 
depends on whether the classical hmit of the considered system is integrable or chaotic. It 
should be mentioned that the characterization of systems as being time-reversal invariant 
means that generically the lengths of periodic orbits are twofold degenerate, i.e. the mul- 
tiplicities g{L/) of lengths should asymptotically, for oo, approach two. There are, 
however, classically completely chaotic time-reversal invariant systems with exponentially 
increasing multiplicities of lengths to which Berry's results do not apply. These systems are 
subsumed under the notion of arithmetical chaos and are the main objects of the present 
investigation. 

A question that might immediately arise is whether or not this universality of the spec- 
tral properties of quantum systems carries over to short-range correlations, i.e. to the level 
spacings distributions P{s). Berry and Tabor demonstrated [24] that generic classically in- 
tegrable systems show a Poissonian level spacing, P(s) = e~^. Thus their quantum energy 
spectra are close to totally uncorrelated ones, since the short-, medium-, and long-range 
correlations are Poissonian. The only difference occurs with the saturation of A^lL; x) for 
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very large L, L ^ Lmax- Concerning chaotic systems, however, there do not exist theoret- 
ical results for the level spacings. But it is generally believed that for these systems P{s) 
may be well described by the distribution of spacings of eigenvalues in an ensemble of large 
random matrices. In random matrix theory (RMT) the statistical properties of the eigen- 
values of random matrices from several ensembles have been thoroughly studied, see e.g. 
[70, 35, 45] and Bohigas' contribution in [41]. For time-reversal invariant chaotic systems 
the appropriate ensemble is the Gaussian orthogonal ensemble (GOE) of real symmetric 
matrices, whose level spacings distribution may be well approximated by Wigner's surmise 
P{s) = |se~'^** Dropping the requirement of time-reversal invariance one has to con- 
sider the Gaussian unitary ensemble (GUE) of complex hermitian matrices. Historically, 
the conjecture that quantum spectra may be described by RMT goes back to Wigner and 
Landau and Smorodinsky, who applied this to the resonance levels of complicated atomic 
nuclei, see [77] for a collection of the original contributions. Bohigas, Giannoni and Schmit 
extended the conjecture of an RMT-behaviour of quantum energy spectra to systems with 
only a few degrees of freedom possessing chaotic classical limits [30], supporting this by 
detailed numerical studies for several systems. In addition, the spectral rigidity A3(L;x) 
for 1 <^ L -C Lmax found by Berry agrees with the RMT-prediction. Again, as in the 
integrable case, the difference lies in the saturation of the rigidity of the actual spectra for 
L — > cxD, as described by Berry's periodic-orbit theory. 

In summary, the quantum spectral properties of classically integrable and classically 
chaotic systems differ in that for small level spacings, s — > 0, the former ones show P{s) ~ 
1 — s, whereas the latter ones behave as P{s) ~ |s. For classically integrable systems 
hence the phenomenon of level attraction occurs, whereas classically chaotic systems excel 
by a level repulsion. On larger scales L rather strong correlations, measured by the spectral 
rigidity A3(L; x), are found in the chaotic case leading to an only logarithmically increasing 
rigidity. In contrast, integrable systems possess nearly uncorrelated spectra; only on very 
large scales correlations do occur. 

Up to now several measures to study the statistical properties of quantum energy 
spectra have been introduced, most of them involving the fluctuating part Nfi{E) of the 
spectral staircase. In sections 3.6 and 4.5 then the strengths of spectral fluctuations as 
measured by Nfi{E) or 1^^{E) play an important role and decide on the applicabihty of 
certain methods of inverse quantum chaology. It hence seems to be worthwhile to study 
Nfi{E) more thoroughly to get hands on quantum energy spectra of classically chaotic 
systems. Thus it seems to be justified to devote some fraction of the following investigations 
of quantum energy spectra to a study of their fluctuations, and to the question how these 
are being influenced by the properties of the classical hmits of the considered systems. 

3 Classical Aspects of Arithmetical Chaos 

In this chapter the classical dynamics of the systems subsumed as showing arithmetical 
chaos shall be introduced. In order to make the whole presentation self-contained, the first 
section of this chapter is devoted to a review of important definitions and results in two 
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dimensional hyperbolic geometry. The second section then contains a discussion of geodesic 
length spectra of hyperbolic surfaces and introduces some useful relations among quantities 
referring to length spectra. Then, in the following two sections, arithmetic Fuchsian groups 
are defined and their geodesic length spectra are studied, leading to the final result on the 
exponential growth of the mean multiphcities of lengths. The results of sections 3.2 to 3.4 
have previously been published in [33] . The next section then is devoted to a discussion of 
pseudosymmetries on arithmetic surfaces. Finally, in section 3.6, statistical properties of 
length spectra are studied. 



3.1 A Brief Review of Hyperbolic Geometry 

The classical dynamical systems to be considered from now on are geodesic flows on hy- 
perbolic surfaces of finite area. In physics terms they are given by single point-particles of 
mass m moving on surfaces of constant negative Gaussian curvature without being subject 
to any external force. The absence of a potential causes the Jacobi metric to be propor- 
tional to the hyperbolic metric the surface has been endowed with. Thus the classical 
trajectories of a particle on such a surface are the geodesies of the hyperbolic metric. In 
the following the geometric setting shall be reviewed in a rather sketchy manner in order 
to recall the necessary notions and to introduce the notations used further on. A rather 
extensive review of the physical and mathematical aspects of the problem may be found 
in [17]. 

A convenient model for hyperbolic geometry in two dimensions is the upper complex 
half-plane Ti. = {z = x+iy; y > 0} endowed with the Poincare metric rfs^ = y'"^ {dx^ + dy"^) . 
In this setting the Gaussian curvature is = — 1 everywhere on Ti. This is dimensionless 
since the internal length scale R has been normahzed to R — 1. Otherwise the metric 
would have to be replaced by ds^ R^ ds^. The hyperbolic distance d{z, w) of two points 
z,w E 7i IS the infimum of the set of lengths of curves connecting z and w. measured 
with ds^. As Ti endowed with ds^ is geodesically complete the infimum is attained for the 
geodesic segment connecting z and w. Its length is given by 

I 1 2 

Z — W 

cosh, diz, w) = 1 -\ . (3.1) 

2 Im z Im w 

There exists an operation of the group PSL{2, M) — SL{2, 1R)/{±1.} on H by fractional 
linear transformations. If 7 = {^d) ^ SL(2, IR), and z E H, then 

az + b 

-fz = —— . 3.2) 

cz + d 

This operation is transitive and compatible with the group structure of SL{2,]R). The 
stabilizing subgroup oi zq = i E T-i is the maximal compact subgroup K = S0{2,]R) C 
SL{2,IR). Thus the hyperbohc plane may also be reahzed as the symmetric space Ti, = 
SL{2, IR)/S0{2, M). Since both matrices 7 and —7 from SL{2,]R) map a 2; G to the 
same image it is the projective group PSL{2, IR), where the centre {±1} of SL{2, IR) has 
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been factored out, that effectively operates on V.. Prom now on the distinction between 
matrices 7 G SL(2, IR) and classes [7] G PSL{2, M) will be dropped and an identification 
of 7 and —7 will be understood automatically. Matrices 7 G SL{2, M) will be chosen such 
that tr 7 > 0. This somewhat sloppy notation seems to be convenient and should not cause 
confusion. 

PSL{2, M) operating on H turns out to be the group of orientation preserving isometrics 
of the Riemannian space {Ti.ds^), i.e. 7 : 7i ^ 7i preserves the lengths of curves on H.. 
Three classes of elements of PSL{2, M) have to be distinguished according to their traces 
(77^1): 

1. Elliptic elements 7 have < tr 7 < 2. Such a 7 has one fixed point in the interior of 
n. It is conjugate within SL{2,R) to an element (4°^^ ^^f^), 1^ G [0,27r). 

2. Parabolic elements 7 have trj = 2. Such a 7 has one fixed point on the boundary 
OH of H {dH — IRU {ioo}). It is conjugate within SL{2, IR) to an element Q f^, 
X eM. 

3. Hyperbolic elements 7 have tr^ > 2. Such a 7 has two fixed points on dH. It is 
conjugate within SL{2, M) to an element 1 ^ . N > 1 is called the norm of 7. 

The geodesies of the Poincare metric on H are the half-circles and the straight half-lines 
perpendicular to the real axis. The two end points of a geodesic are the fixed points of 
a (unique) hyperbolic 7 G PSL{2, M). This 7 maps the geodesic onto itself, which is 
therefore also called the invariant geodesic of 7. 

A discrete subgroup F C PSL{2, M) gives rise to the orbit space r\H = {Tz; z G 7i}. 
This means that an equivalence relation on H is introduced by an identification of points 
that are related by a P-transformation. As the Poincare metric is P-invariant it induces a 
metric on the orbit space T\H, which will also be called Poincare (or hyperbolic) metric. 
If P contains no elliptic elements, i.e. no 7 G P has a fixed point on Ti, r\H is a regular 
surface. If, however, elliptic elements are present in P, P\7^ will be a regular surface outside 
the respective fixed points. Including these turns r\H into an orbifold. Despite this slight 
complication any P\7i, irrespective of a possible existence of orbifold-points, will be called 
a hyperbolic surface. If V\H has finite area (measured with the Poincare metric), P is 
called a Fuchsian group of the first kind, otherwise a Fuchsian group of the second kind. 
The latter ones will be excluded from the further investigations. If a Fuchsian group of the 
first kind contains no parabolic elements, T\H is compact. Therefore such a P is called 
cocompact. 

The abstract orbit space V\H may be realized explicitly as a fundamental domain T 
of F on 7i. This is a simply connected region JF c 7i such that the interior of T contains 
no two F-equivalent points and the union of all F-translates of JF covers all of ?i. The 
boundary dT of T consists of segments that have to be identified under the operation of 
P appropriately in order to yield a closed surface. If P is strictly hyperbolic, i.e. it contains 
besides the identity only hyperbolic elements, r\H is a compact surface of genus g > 2. 
Thus it is topologically equivalent to a sphere with g handles. It is then always possible 
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to find a fundamental domain T that is bordered by geodesic segments (a 4g'-gon). 
Since the Gaussian curvature is K = —1 the GauB-Bonnet theorem yields for the area of 
r\7i area(r\7i) = area(jF) = 47r((? — 1). If F is not cocompact the F-conjugacy classes of 
parabolic elements are called cusps. A fundamental domain then extends to infinity (i.e. 
to dH). For each cusp dT contains a point of dH — ML) {ioo}. 

A convenient presentation of a Fuchsian group F may be given in terms of generators 
related to a given fundamental domain J-", and relations among them. As generators one 
can take those transformations that identify pairs of geodesic segments of dJ-", together 
with their inverses. If F is strictly hyperbolic, dJ- has 4g, g > 2, components. Thus F 
has Ag generators which are conventionally denoted as ai,bi, . . . , ag,bg, , . . . ,bg^ . One 
frequently chooses the order of the generators such that the one relation they obey is 

oi^iar^^r^ . . . agbga~^b~^ = 1 . (3.3) 

This presentation seems to be well suited for exphcit calculations since the group may 
be constructed explicitly in terms of SL{2, iR)-matrices once a fundamental domain JF is 
given. (See e.g. [12] and appendix B for examples.) The procedure will then be as follows: 
draw an arbitrary geodesic ig-gon {g > 2) with area(jF) = 4:n{g — 1). It only has to share 
the further property that the geodesic segments have to come in pairs of identical lengths. 
Construct the 2g matrices (from SL{2, IR)) that identify the pairs of edges. Together 
with the 2g inverse matrices 4g generators for the desired Fuchsian group F have been 
found. F itself can then be obtained by forming all possible products of the generator 
matrices. (These are also called words in the generators.) In general the Ag generators 
will be "arbitrary" matrices, which means that it does not seem possible to obtain an 
explicit law that determines their matrix entries. Forming words in the generators will 
also yield matrices with seemingly arbitrary entries. In general there will thus not exist 
some (algebraic) structure for the matrix entries of the group elements. A consequence 
of this fact for practical purposes hence is that there seems to be no explicit enumeration 
scheme for the group elements other than the generator method just described. 

This observation is indeed true for "generic" Fuchsian groups. An exception to this 
rule, however, does exit. This is formed by the class of arithmetic Fuchsian groups. For 
these the arithmetic structure appearing in their definition determines their matrix entries. 
The basic and most prominent example for an arithmetic Fuchsian group is the modular 
group SL{2, Z). It consists of all 2 x 2 matrices with rational integers as entries and being 
of unit determinant. Leaving aside the determinantal condition the determination of the 
matrix entries could not be easier. For general arithmetic groups the definition is a bit 
more involved and will therefore be postponed to section 3.3. But still, the characterization 
of the set of traces of group matrices is rather simple compared to the generic case. 

For apphcations in the context of quantum chaos the knowledge of the set of traces 
of elements of a Fuchsian group is of major importance, since it is tightly connected to 
the knowledge of the geodesic length spectrum. The latter in turn is the input on the 
classical side of periodic-orbit sum rules. Since the set out of which the traces of an 
arithmetic Fuchsian group have to be taken is known (see section 3.4), it is possible to 
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develop an algorithm that allows to calculate the geodesic length spectrum (including 
multiplicities) for an arithmetic group up to a certain length completely, see [4, 81, 68, 73]. 
For non-arithmetic groups there seems to be no better way to calculate the length spectrum 
numerically other than forming words in the group generators up to a given number of 
generators per word. This method, however, does not yield the length spectrum completely 
up to some value of the length. 

Because of their importance for periodic-orbit theory the following section will be de- 
voted to a discussion of geodesic length spectra connected with Fuchsian groups. 

3.2 Geodesic Length Spectra of Hyperbolic Surfaces 

By the very construction of hyperbolic surfaces T\H, where F is a Fuchsian group of the first 
kind, the hyperbolic plane Ti. is the universal covering space of F\7i. A closed geodesic c on 
r\H lifts to a geodesic c on Ti., which is invariant under some hyperbolic 7c G F. For any 
7 e F the transformation 77c7~^ has some invariant geodesic c-y on Ti, that projects down 
to the same closed geodesic c on F\7^. There exists in fact a one-to-one correspondence 
between F-conjugacy classes of hyperbolic elements of F and closed geodesies on F\7i. The 
hyperbolic length / of c can be related to 7c as follows: the lift c of c (being a segment of c) 
on 7i connects two points z and z' on H which are F-equivalent. Thus there exists a 7 G F 
such that z' — 72;. Then I — d{z,z') = d{z,jz). Since c is a geodesic segment d{z,'jz) 
minimizes the lengths of all curves on H connecting z and 72;. Varying z continuously along 
c cannot change the choice of 7. Since all these z's and 72;'s lie on c, this has to be the 
invariant geodesic of 7, and hence 7 = 7c. To compute / = d{z, 7c-s) one can conjugate 7c in 

PSL{2, IR) to obtain the diagonal matrix 7' = 1 ^ = g'jcg~^, for some g e SL{2, IR) 

and N > 1. Then / = d{z,^cz) = d{z, g^^'-f'gz) = d{gz,^'gz). The invariant geodesic of 7' 
is obviously the imaginary axis, and choosing gz to be on it one obtains from (3.1) that 
cosh / = 1 + 2^(A^ - If. This finally leads to 

2 cosh ^ — tr^' — tr^. (3.4) 

Because of the one-to-one correspondence between conjugacy classes of hyperbolic elements 
7 e F and closed geodesies c on r\H we denote the length / of c by / = /(7), where 7 is 
some representative of the conjugacy class {7}r = {777~^; 7 £ T} related to c. 

The set of lengths of closed geodesies on F\7i is called the geodesic length spectrum 
£(F) of r\n, or briefly, of F, 

£(F) = {Z(7); 7 e F hyperbolic} . (3.5) 

Introducing the notation >C(F) = {Zi < ^2 < ^3 • • •} the counting function J\f{l) for the 
geodesic length spectrum is 

Af{l):^4{n; Z„ < Z} . (3.6) 

In general a length spectrum of a hyperbolic surface will be degenerate, i.e. there exist 
several closed geodesies on the surface of the same length. If e.g. F is strictly hyperbolic. 
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no 7 e r different from the identity is conjugate to its inverse, so that every length at least 
occurs twice. Interpreted in physics terms this means that the dynamical system possesses 
a time-reversal symmetry, since the geodesic corresponding to the inverse of a hyperbolic 
transformation is the original classical trajectory traversed backwards in time. In general 
the multiplicity of any / G C{T) will be denoted by g{l) G ]N. The counting function 
including multiplicities is the classical staircase function 

J\f(l) := # { {7}r; 7 e r hyperbohc and /(7) < /} . (3.7) 

Asymptotically, for I — > oo, the magnitude of Af{l) is universally determined by Huber's 
law [58, 51], also known as the prime geodesic theorem (PGT), 

J\f{l) - Ei{l) ~ - , l^oo . (3.8) 

t 

The remainder to this asymptotic relation will in more detail be dealt with in section 3.6. 
Ei{l) is related to the logarithmic integral li{x) — P Jq by Ei{\ogx) — li{x), see e.g. 
[39]. {P J dt . . . thereby denotes the principal value of the integral.) Comparing (3.8) with 
(1.5) yields that the topological entropy for geodesic flows on hyperbolic surfaces is r = 1. 

A closed geodesic is called primitive, if it is not a multiple traversal (> 2) of some 
other closed geodesic on the same surface. The corresponding conjugacy class {7}r is then 
also primitive, i.e. a 7' e {7}r is not a power (> 2) of some other (hyperbolic) element 
of r. The primitive closed geodesies give rise to the primitive geodesic length spectrum 
Cp{r) = {Ip i < lp^2 < lp,3 < ■ ■ ■}■ All quantities deflned for the full length spectrum 
can be introduced analogously for the primitive one as well. To distinguish the respective 
quantities one introduces an index p in all notations referring to the primitive length 
spectrum. The PGT (3.8) remains true also ior J\fp(l). Since 

[i/h] [i/h] 

Ml) = E 1 = E E 1 = E Mi^M , (3.9) 

ln<l r=l rlp^n<l ''=1 

and Ap(/) is positive and monotonically increasing, the counting functions for C{T) and 
£p(r), respectively, show the same asymptotical behaviour, 

Sf{l) ^ Afp{l) , 1^00. (3.10) 

In the following sections the multiplicities gp{l) of lengths of primitive closed geodesies 
will play a major role. As it seems to be impossible to determine the multiplicities explicitly, 
the only quantity one can lay hands upon appears to be an average < gp{l) >, and its 
asymptotic behaviour for / 00. In number theory [64] two functions f{n) and h{n), 
n G -2V, are said to be of the same average order, if 

N N 

E/H~EM^), A^^oo. (3.11) 

n=l n=l 
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Led by this definition a mean multiphcity < gp{l) > is introduced as a continuous function 
of I that is of the average order of gp{l), 

Kil) = E 9v{lv,n) - E < 9vilp,n) >= f^dMp{l') < 9p{l') > , (3.12) 

for Z — > oo. Substituting the multiphcities by their mean hence does not violate the PGT. 
Since N'pil) ~ Ei{l), I — > oo, one observes the asymptotic relation 

-1 

, Z — > oo , (3.13) 



< 9pil) >~ J 



~dr 



by differentiating (3.12). This result may be interpreted in more visual terms as follows. 
Np{l) is a step function with step-width A/„ := Ip^^ ~ lp,n-i and step-height Qpilp^n) at 
I — Ip^n- A mean of the staircase N'pil) should then have a slope of < gp{l) > / < 
Al >. As < Al is the mean density of primitive lengths it is asymptotically given 
by Thus ~< gp{l) > which is equivalent to (3.13). From this relation one 
concludes that one has to gain information on the asymptotics of the counting function 
N'p{l) of different lengths of primitive closed geodesies in order to be able to determine the 
asymptotic behaviour of the mean multiplicity < gp{l) >. Notice that the same procedure 
may also be carried through for the full length spectrum instead of the primitive one. Since 
both staircase functions N'il) and J\{1) share the same asymptotic behaviour as N'pil) and 
Mp{l), respectively, the mean multiplicities agree asymptotically. 

Next we are interested in relating the counting functions J\f^^\l) and for I — > oo. 

Thereby M^^\l) shall correspond to the length spectra £p(Fj), where Fi is a subgroup of 
index d > 1 in the Fuchsian group of the first kind F2. This means that F2 decomposes 
into d cosets of Fi according to 

F2 = FiUFi7iU...UFi7rf_i , (3.14) 

where 71, ... , ^d-i are elements of F2, but not of Fi. Let Cp{T2) be given as 

>Cp(F2) = {lp,i < lp,2 < lp,3 <■■■}; (3.15) 

Cp{Ti) shall now be determined in terms of the Zp,„'s. To this end one observes that if 
7 e F2, then there exists a A; e W, such that 7^^ e Fi. To see this take an arbitrary 7 G F2 
and form Ume^^il"^ C F2. The union on the l.h.s. cannot be disjoint, since according to 
(3.14) Fi is of finite index in F2. Therefore there exists a 70 G F17'' fl F17* for some pair 
r ^ s. Thus 7'^'"^ and 7*"^ lie in Fi. Choosing k = \r — s\ proves the assertion. 

Let now 7^ be a primitive hyperbolic element of F2, hence /(7p) e £p(F2). Then either 
i) 7p e Fi: 7p is also primitive hyperbolic in Fi and thus li'jp) G £p(Fi), or ii) 7^ ^ Fi. 
Then, by the above remark, there exists a /c G W with 7^ G Fi. The minimal such k 
takes care for /(7p) = kl{jp) to be an element of >Cp(Fi). Altogether, the primitive length 
spectrum of Fi can be characterized as 

>Cp(Fi) = {ki Ip^i, k-2 lp,2, ks ip,3, . . .} (3.16) 
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with positive integers kj. Since the /c^'s may take arbitrary values in W, this enumeration 
of elements of CpiVi) will in general not be an ordered one. 

The determination of the asymptotic behaviour of Mp^\l) in terms of that of N'^^\l) 
requires to know how often (asymptotically) a certain value of the fc^'s in (3.16) occurs. 
This question may be answered with the help of the decomposition (3.14) of r2 into cosets 
of Fi. On average, a fraction of ^ of the elements of r2 are also elements of Fi. Thus going 
through the conjugacy classes of hyperbolic elements of r2 and picking one representative 
from each class yields a fraction of ^ of the latter ones to lie in Fi. This statement is 
meant in the following sense: one chooses an ordering of the conjugacy classes according 
to the corresponding lengths / of closed geodesies (equivalently, in ascending order of their 
traces). Then, for / — > oo, the respective fraction of elements of Fi approaches Hence 
a fraction of ^ of the /c/s equals one. Proceeding further in the same manner yields a 
fraction of (1 ~ ^)^ to be two, and so on. 

By (3.15) and (3.16) obviously J\f^^\l) < N'^^\l). The above remark on the fraction of 
kj = 1 yields N'^^l) > iA/^^^)(/) in the limit / — > oo. This relation has to be an inequality 
rather than an equality because there might occur some kj > 2 with kjlpj < I. The 
number oi kj = k > 2 such that klp^j < I is, however, bounded by N'^^\^). But this is 
asymptotically dominated by the contribution coming from k — 1. One therefore observes 
the asymptotic relation 

•^i'nO-^A^-f (0, ^-oo. (3.17) 

Again, a corresponding relation also holds for the counting functions M^^\l) and 
of the full length spectra. 

Finally, the asymptotic relation of the counting functions for length spectra of two 
commensurable Fuchsian groups shall be dealt with, since this will be needed later when 
treating arithmetic Fuchsian groups. One recalls that two subgroups Hi and H2 of a group 
G are called commensurable, if the intersection Hi D H2 is a subgroup of finite index in 
both Hi and H2. Let therefore Fa and F;, be two commensurable Fuchsian groups of the 
first kind. Fq := F^ fl F^ shall be of index da < 00 in Fa and of index < 00 in F^. By 
(3.17) the counting functions are related through 

J^i'\l) - ^-^fW, (3.18) 
for I ^ 00. Prom this the asymptotic relation 

A^^HO-^^AAf (0 , 1^00, (3.19) 

easily follows. Recalling the asymptotic behaviour (3.13) of the mean multiphcities of 
primitive lengths yields 

<^nO>~^ <^i'nO>, ^-oo. (3.20) 
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Therefore, given two commensurable Fuchsian groups, the respective mean multiphcities 
of primitive lengths of closed geodesies are asymptotically, for / — > oo, proportional to one 
another. The factor of proportionality is given by the ratio of the indices with which the 
two groups contain their intersection as a subgroup. 

So far, the discussion of length spectra has been completely general in that it has 
been valid for all Fuchsian groups of the first kind. A major objective of the present 
investigation, however, is to stress the importance of distinguishing arithmetic Fuchsian 
groups from non- arithmetic ones. Concerning length spectra, the difference occurs with 
the behaviour of the mean multiplicities < gp{l) > as / — oo. For arithmetic groups 
r it shall be shown in section 3.4 that < gp{l) >~ cr / — > oo, where cr is some 
constant depending on the specific group F. This behaviour is exceptional and constitutes 
a major property by which the arithmetical systems excel among general chaotic dynamical 
systems. It is indeed known that for general hyperbolic surfaces gp{l) is always unbounded 
[57, 78, 36]. The proof of this proceeds algebraically; Buser [36], however, geometrically 
constructs explicit examples of degenerate closed geodesies, where the multiplicities of the 
respective lengths cannot be accounted for by symmetries of the surface. He can estimate 
the multiphcity of a length / to which the construction is apphed to be of the order for 
a = 1^ and / oo. Aurich [3] computed the lower parts of length spectra for several 
arbitrarily chosen compact surfaces of genus two and never observed a multiplicity that 
could not be explained by symmetries in the computed /-range. Thus it appears that 
high multiplicities (having values beyond the expectation based on symmetries) show up 
rather scarcely. All this information indicates that for non-arithmetic Fuchsian groups the 
multiplicities of lengths of closed geodesies do not grow exponentially. The latter seems to 
be characteristic of surfaces with arithmetic Fuchsian groups. 

Without giving the precise definition of arithmetic Fuchsian groups at this stage we 
would now like to try to illustrate the difference between arithmetic and non- arithmetic 
groups in order to give an intuitive understanding of the different properties of their re- 
spective length spectra. As described in section 3.1 a presentation of a Fuchsian group F 
may be yielded in terms of generator matrices (?!,... ,gn, where the gfj's define fractional 
linear transformations on Ti that identify pairs of edges of the fundamental domain of 
F. If ^ (and thus also F) is arbitrary, then gi, . . . ,gn will be matrices with "arbitrary" real 
entries. Consequently the set of matrix entries of all elements of F has no obvious alge- 
braic structure other than that the entries are composed of those of the generator matrices. 
Since concerning length spectra the traces of group elements are the objects of interest, a 
possible algebraic structure of the set of traces trT = {tr 7; 7 e F} of elements of F would 
influence jC{T). As an example of an arithmetic group the modular group F^od = SL{2, 
should serve. From the theorem proved in [66] one can conclude that tr T^nod — ^ (see also 
[81]). Therefore trT mod has a nice algebraic structure in that it is just the ring of rational 
integers. According to (3.4) one can hence enumerate the lengths of closed geodesies on 
^mod\H by n = 2 cosh(/„/2), n = 3, 4, 5, . . .. Since J\f{ln) = n — 2, one finds the asymptotic 
law N{ln) ~ e'"/^, n — > cxD, for the modular group. By (3.13) the multiphcities of lengths 
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therefore grow asymptotically as 



J/2 

<g{l)>r^2—, l^oo. (3.21) 

It was now observed in [32, 27] that such an exponential increase of the multiplicities is a 
common feature of arithmetic groups and may serve as a characterization of them. The 
example of the modular group shows most clearly that this law derives from the fact that 
the condition on the set of traces tr Tmod to be the "rigid" set Z forces the number Af{l) 
of different lengths to grow only like const, e'/^. One can transform this condition also into 
one which applies directly to trT. If Af{l) ~ ce'/^, / oo, then the number of lengths in 
an interval [/, / + 2e~'/^] is asymptotically given by c. If one returns to tr T, then the above 
interval corresponds asymptotically to the interval [t,t + 1], where t — 2cosh(//2) ~ e'/^. 
Thus the number of traces of elements of F in the interval [i, i + 1] is asymptotically given 

by 

#{trrn [t,t+ 1]} ~ c , t^oo. (3.22) 

This condition shows more explicitly that the condition on tr T enforces the traces to form 
a "rigid" discrete subset of the real line. The condition ^{trV fl [t, t + 1]} < c(r) for some 
constant c(r) has most recently been proven by Sarnak and Luo [80] to be fulfilled for all 
arithmetic Fuchsian groups F and has been called the bounded clustering property by them. 

Our plan now is to derive (3.21) for general arithmetic groups. Since this needs a 
lot of techniques from number theory, at first a brief review of relevant notions and facts 
concerning algebraic number fields and quaternion algebras shall be given. This then allows 
to state the precise definition of arithmetic Fuchsian groups. Only after this has been done 
it will be possible to study the traces occurring for arithmetic groups in more detail. 

3.3 Arithmetic Fuchsian Groups 

This section contains a collection of notions and facts from number theory that are needed 
to state the definition of arithmetic Fuchsian groups. More background on the algebraic 
foundations and on number theory may be found in [55, 98, 34]. As references to find 
additional information on quaternion algebras and arithmetic Fuchsian groups [86, 102, 71] 
may be consulted. 

An extension K of finite degree n of the field of rational numbers Q is a field that 
contains Q as a subfield and, viewed as a vector space over (Q (in the obvious manner), is 
of finite dimension n. Let (Q[a;] denote the ring of polynomials in a variable x with rational 
coefficients, a ^ K will be called algebraic, if it is a zero of some polynomial from (Q[a:]. 
The minimal polynomial of a is the (unique) element in (Q[a;] of lowest degree, and with 
leading coefficient one, that has a as a root. The field K is called an algebraic number 
field, if every a E K is algebraic. Every extension X of Q of finite degree is known to be 
algebraic. 

If M is some arbitrary subset of K, Q(M) is defined to be the smallest subfield of K that 
contains both M and Q. It is given by all values of all polynomials in the elements of M 
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with rational coefficients, and all possible quotients thereof. Q(M) is called the adjunction 
of M to Q. One can now show that every algebraic number field K of finite degree over 
(Q can be realized as an adjunction of a single algebraic number a & K to % therefore 
K=^{a). 

Since K isa vector space of dimension n over (Q, the n+1 algebraic numbers 1, a, . . . , a" 
have to be linearly dependent and thus have to obey a relation 

a„Q;" + • • • + Old! + oo = (3.23) 

with rational coefficients Oj and a„ ^ 0. Normalizing the leading coefficient to one leaves 
a as a root of an irreducible polynomial fa{'x) G Qf^] of degree n. fa{x) is the minimal 
polynomial of a. Since {1, a, ... , a"^^} may serve as a basis for K over Q, any x E K may 
be expanded as a linear combination of powers of a up to the order n — 1, 

X = 6„_ia"-^ + h 6ia + 6o , (3.24) 

with rational coefficients 6j. 

The polynomial fa{x) has n different complex roots ai, . . . , {ai = a). One can thus 
define n different homomorphisms ipi : K — > (D, i = l,...,n, that leave Q invariant, by 

iPi{x) := + • • • + + 6o , (3.25) 

(pi{x) = X. The ipiS are called the conjugations of K. If all images of K under these 
homomorphisms are contained in the real numbers, in other words if fa{x) has only real 
roots, then K is said to be totally real. 

On (Q the usual absolute value Uii^x) = a; G (Q, introduces a topology, which is, 
however, not complete. The n conjugations offer n distinct ways to embed K into M. 
Thus n different absolute values Ui arc given on K by i'i{x) := |9?i(a;)|, and these can be 
used to complete K to Ki,^ = M. The ViS are also called the (archimedean) infinite primes 
of K. 

All algebraic numbers in K whose minimal polynomials have coefficients in the rational 
integers Z form a ring TZk, which is called the ring of integers of K. An element x G TZk 
is also called an algebraic integer. In K at most n algebraic numbers can be linearly inde- 
pendent over (Q. (This is equivalent to saying that these numbers are linearly independent 
over Z.) Suppose now that Ai, . . . , A^, m < n, are linearly independent numbers from K. 
Then all linear combinations of the Aj's with integer coefficients form an additive abelian 
group. This is called a Z-module of rank m. A ^-module o C A' of the maximal possible 
rank n that at the same time is a subring of K is called an order oi K. Since we understand 
a ring to contain a unity, every order a <Z K contains the rational integers Further 
it is known that there exists a maximal order in K that contains all other orders, and 
that this maximal order is just the ring of algebraic integers TZk- An order o possesses a 
module-basis of n algebraic numbers uJi, . . . ,uJn that are linearly independent over Z and 
hence also, equivalently, over Q, 

= ^ui e • • • e . (3.26) 



27 



The discriminant of o is defined to be Dk/q{o) :— [det{(pj{uji))]'^ ^ 0. In complete analogy 
one can also define a discriminant for any ^-module of rank n in K. 

Another important notion, to be introduced now, is that of a quaternion algebra. In 
doing so, we will mainly follow [102, 86]. An algebra A over a field K is called central, if 
K is its centre; it is said to be simple, if it contains no two-sided ideals besides {0} and A 
itself. A quaternion algebra then is defined to be a central simple algebra A of dimension 
four over K. In more explicit terms A may be visualized as follows: the elements of a 
basis {1, a, P, 7} of A over K have to obey the relations 7 = afi — —(5a, = a, (3"^ = b; 
a,b & K\{0}. Any X & A may then be expanded as 

X — xo + xia + X2P + xzol(3 , (3.27) 

with Xq, . . . ,x^ G K. On A there exists an involutory anti-automorphism, called the 
conjugation of A, that maps X io X :— xq — x^a — X2P — x^a(3. Thus X — X and 
X - Y = Y ■ X. The conjugation enables one to define the reduced trace and the reduced 
norm of A, 

trA{X) := X + X = 2xo , 

nA{X) := X ■ X ^ xl-xla-xlh + xlah . (3.28) 

If A is a division algebra, i.e. if every X ^ Q m A possesses an inverse, nA{X) = implies 
X = Q. The inverse is then given by X~'^ = — jit^X. 

A .^-module O <Z A ol (the maximal possible) rank An that also is a subalgebra in A 
is called an order of A. The introduction of a module-basis {ti, . . . , T4„} turns the order 
into 

O ^ Zti®---® ^Tin ■ (3.29) 

We further introduce the group of units of norm one :— {e G O; G O, ua^s) = 1 }■ 
A well-known example of a (division) quaternion algebra is given by Hamilton's quater- 
nions 

^-■=[[-u, ^) ; ^'^^^} ■ (3.30) 

iff is a four dimensional iR-subalgebra of M(2,(D), the algebra of complex 2x2- matrices, 
characterized by the parameters a — h — —1. The subgroup of elements of reduced norm 
one is just S'f/(2,(D). An even simpler example of a (non-division) quaternion algebra over 
IR is M(2, Et). In fact, M and M(2, IR) are the only quaternion algebras over IR. 

A classification of quaternion algebras over K can now be achieved by looking at the 
corresponding algebras over IR with the help of the n completions Kj.^ = IR. Define Ai :— 
A ®Q Ki,^ = A ®Q IR, which is a quaternion algebra over M. Hence it is either isomorphic 
to HI (if it is a division algebra), or to M{2,IR) (if it is a non-division algebra). For the 
definition of arithmetic Fuchsian groups (see [102, 94]) we consider the case Ai = M{2, IR) 
and Ai IH ioT i — 2, . . . ,n. Therefore there exists an isomorphism 

p: A®qIR — > M{2,IR)®IH®---®IH , (3.31) 
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where there occur n — 1 summands of M. pj will denote the restriction of p to A followed 
by a projection onto the j-th summand in (3.31). The several reduced traces and norms 
for X e A in (3.31) are related by 

trpi(X) = trA{X) , 

detpi{X) = nA{X) , 

trjHPjiX) = <pj{trA{X)) ^ <pj{tr pi{X)) , (3.32) 

nmipjiX)) = ipj{nA{X)) = ipj{det pi{X)) , j = 2, . . . , n . 

The image of A under pi in M(2, M) may also be expressed in more explicit terms by using 
the basis {1, a, /3, a/?} for A, see (3.27): pi(l) is the 2x2 unit matrix; pi{a) and pi{/3) 
may be represented, by using the parameters a,b > 0, as 



- V _ r. . Pm - .ru V ■ (3-33) 





For X — xo + xia + X2P + x^a^ e A the matrix pi{X) in this representation takes the 
form 

(X) = ( ^° X3\/ab \ , . 

\ x^yfb — xj,\foh xq — x\^Ja J 

We are now seeking for a subset in A whose image under pi in M(2,1R) gives a Fuchsian 
group r. Therefore pi^{T) must be a discrete multiplicative subgroup of A. Furthermore, 
for pi{X) = 7 G r the condition det^ = 1 must be fulfilled. Thus by (3.32) n^(X) = 1 
has to be required. Hence we are led to look at groups of units of norm one of orders 
O C A. Regarding their images under pi one finds in [86, 94] the following 

Proposition: Let A be a quaternion algebra over the totally real algebraic number field 
K of degree n. LetOcAbean order and be its group of units of norm one. Then 
T{A, O) := pi(C^) is a Fuchsian group of the first kind. Moreover, T{A, 0)\H is compact if 
A is a division algebra. A change of the isomorphism p in (3.31) amounts to a conjugation 
of r(^,e)) in SL{2,1R). 

Remark: This proposition can be traced back to a more general theorem of A. Weil that 
deals with an adelic setting, see e.g. [103, 86]. For the proposition to be true it is essential 
that the quaternion algebra is such that on the r.h.s. of (3.31) there appears exactly one 
factor of M(2, iR) and n — 1 factors of iff. (3.31) therefore is an integral component of the 
definition of arithmetic Fuchsian groups. 

The proposition now tells us that we have found what we were looking for: a class of 
arithmetically defined Fuchsian groups. 

We are aiming at counting the numbers of distinct primitive lengths to gain information 
on the mean multiplicities in the length spectra derived from the Fuchsian groups under 
consideration. For this purpose (3.20) allows to enlarge the class of groups appearing in 
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the proposition a httle. 

Definition: A Fuchsian group T that is a subgroup of finite index in some T{A, O) will 
be called a Fuchsian group derived from the quaternion algebra A. (For T{A, O) the short- 
hand phrase quaternion group will also be sometimes used.) A Fuchsian group F that is 
commensurable with some V{A,0) will be called an arithmetic Fuchsian group. 

In the preceding sections the modular group Tmod always served as an example of 
an arithmetic Fuchsian group. It will now be seen that Vmod fits into the scheme just 
introduced. The number field to be considered is X = (Q, and hence n — 1. In Q the 
ring of integers is of course TZk — which is also the only order in K. The relevant 
quaternion algebra is the simplest one can think of, namely the one characterized by the 
two parameters a = 1 = h. Thus A is simply the matrix algebra M(2,(Q), which clearly 
is not a division algebra. This is in accordance with the non-compactness of the surface 
^mod\^- The order O <Z A (see (3.29)) that determines Tmod is characterized by the 
X-basis {ti, . . . , t^, 

^'^[\ o) ' ^^=(o I) ' ^'^[\ o) ' ^^ = (o ?) ■ ^^-^^^ 

Therefore, O = M(2, Z) and = SL{2, Z). The modular group is of course well-studied. 
A lot of information about this group and the spectral theory on the surface T^o^H can 
be found e.g. in [96]. 

As a second example let us introduce the regular octagon group F^eg. This is a strictly 
hyperbolic Fuchsian group that leads to the most symmetric compact surface Treg\H of 
genus two (see [17, 10, 4]). Freg is a subgroup of index two of the quaternion group T{A, O) 
that is defined over the number field K — (Q(-v^) of degree n — 2. The ring of integers of 
this field is TZk — ^[v^] = {m + n\f2\ m,n e Z}. The module-basis {101,102} = {1, \/2} 
for TZk may also serve as a basis for K over (Q. The only non-trivial conjugation of K is 
given by 

ip2{p + qV2) = p - qV2 , p, g G Q . (3.36) 

The quaternion algebra A necessary to define T{A, O) is determined by the two parameters 
a — 1 -\- \pl and h — \. The order O <Z A can be characterized by giving the .^-basis 
{ti . . . , Ts} (see (3.29)). In the present case this is {lo\ • 1, . . . , a;2 • Q;/3}, so that an element 
7 = px{X) for X e looks hke 

^^i a;o + Xi^l + V2 X2 + X3v^l + v^ 

\X2- X2,\Jl + Xq - Xi\jl + a/2 

with Xi — rrii + nj\/2, mi,ni e The quaternion group T{A,0) now consists of all 
matrices 7 of the form (3.37) with det^ — 1. The regular octagon group Freg is character- 
ized by the fact that mo has to be an odd integer. In [76] it is shown that by adjoining 
the additional (elliptic) matrix S — (i'lo) to Treg one obtains this quaternion group via 

r{A,o)^ 
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3.4 Multiplicities in Length Spectra for Arithmetic Fuchsian Groups 

The problem of this section is to characterize the set of traces trV of elements of an 
arithmetic Fuchsian group F, and then to determine the asymptotics of the number of 
distinct lengths corresponding to it. It, however, suffices to concentrate on quaternion 
groups, since every arithmetic group is by definition commensurable to a quaternion group. 
Then (3.20) can be used to obtain the desired asymptotics of the mean multiplicities for 
the given group from that for the related quaternion group. 

Let therefore F = F(A, O) be a quaternion group over the algebraic number field K 
of degree n as described in the preceding section. For X = xq + xia + X2P + x^aP G O 
denote Pi{X) = 7 G F. By (3.34) one sees that |tr7 = xq E 0\k- Since the coefficients xq 
of elements X G C will play a major role in the following, the set of these will be given a 
name. Let therefore be 

M := {xo; X = xo + xxa + X2(3 + xzal3 G O) . (3.38) 

Then 

trV^ ItaO^ C tTAO = 2M . (3.39) 

The inclusion ItaO^ C tvAO will in general be a proper one and we will return to this 
problem later. The aim now is to determine the number J\p{V) of distinct primitive lengths 
on T\H for / — 00. By (3.4) one hence has to count the number of distinct traces in F 
with 2 <tr-f <2R, R:= cosh(//2) 00. 

First we want to describe the set a little bit further. To this end we have to introduce 
some more notation. Let {cui, . . . ,a;„} be a basis for K (as a vector space) over Q. With 
the help of the basis {1, a, (3, a(3} of A over K then {xi, ■ ■ ■ , X^n} '■— {^1 ' ^, ■ ■ ■ i^n ' 
is a basis of A over Q. On the other hand the module-basis {ri, . . . , r4„} of O (see (3.29)) 
consists of 4n linearly independent (over Z as over (Q) elements of A and thus may also 
serve as a basis for A over Q. The two debases of A are therefore related by 

4n 

where (Mij) G GL(4n,(Q). The order O (Z A then takes the form 

4n 4n 

O^Zj2M,jXj®---®^J2 ^4n,j Xj (3.41) 
3=1 3=1 

after inserting (3.40) into (3.29). As the centre i^T of A is spanned by {xi,...,Xn} — 
{a;i, . . . ,u!n}, it turns out that 

n n 

M^^Yl ^ ^^Yl ojj . (3.42) 



(3.40) 
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Obviously, is a ^-module in K. Since (M^j) e GL(4n,(Q), out of the An algebraic 
numbers YJj=i Mi jUj, i — 1, . . . ,4n, n are linearly independent. One can therefore choose 
the module-basis {/ii, . . . , among them, 

M = ^//i ® • • • e . (3.43) 

In general, however, A4 is not a subring and hence no order in K, because the multiplication 
in it need not close. But in [94] one finds that trV = 2Ai is contained in the ring TZ^ of 
integers of K. Defining /ij := for i = 1, . . . , n then yields 

2M = Zfii ® • • • ® Zfin C TZk . (3.44) 

By (3.4) this means for the geodesic length spectrum C{T) that 2cosh(//2) e 2M. C TZr 
is an algebraic integer for every I e >C(r). 

So far, trAO = 2Ai has been described in an algebraic way. The object of interest, 
however, is trF = tr^O^. The problem thus is the following. Let an Xq € be given, 
which means that there is at least one X ^ O with X = XQ+Xia+X2P+X3a(3. There might 
be, however, several X^ e O all sharing the same first coefficient xq. In a more formal 
way this can be described using the bases (3.41) and (3.42). According to (3.42) xq can be 
represented as Xq — Yli=i Z)j=i i^iMijUj, where all rj e The choice of (ri, . . . , r4„) e Z^'^ 
is not imique, since Ai has only rank n. Thus one can vary (ri, . . . ,r4„) over a certain 
subset of Z'^"' without changing xq. X = J^ij^iTiMijcUj, however, does change under 
this variation. The so resulting (discrete) set oi X e O then is the set of X^s mentioned 
above. The question now is whether there exists some X e O amongst the ones all sharing 
the same first coefficient Xq that has a reduced norm n^(X) = 1. If the answer is in the 
affirmative, then also 2xq G tr^O^ = tr F. But in order to decide on this question one has 
to know whether there is a solution to 

—xfa — x\h + x\ah — 1 — x^ (3.45) 

for a given xq E M.ui the three variables {xi, X2, x^) so that X — xo + xia + X2P + x^ajS G 
O. The indefinite quadratic equation (3.45) defined over the given non-trivial domain of 
variables appears to be difficult to deal with in full generality. In the following this subtle 
problem shall be avoided, and the determinantal condition nA{X) = 1 for X E O shall 
be replaced by a weaker auxiliary requirement. The price to pay for this will be that a 
hypothesis will have to be introduced below without which no result could be obtained. 

Any X G (9^ is characterized within O by the condition n^(X) = 1. By (3.31) and 
(3.32) this implies that n]H{pj{X)) = (pj{nA{X)) = 1 for j = 2, . . . , n. Therefore Pj{X) G 
SU{2,€) and hence trmpj^X) = ^jitvy) G [-2, +2] for j = 2, . . . ,n. We will now call 

tr iT -.^ {2xo; xo e M, |(^,(xo)| < 1, j - 2, . . . , n } (3.46) 

the idealized set of traces of F. Instead of the inclusion trAO^ C tr^O we arc thus now 
considering trjT C tr^O. By the very construction of the idealized traces it is clear that 
trAO^ C triT. From this one obtains the chain of inclusions 

F = trAO^ C tnF C trAO . (3.47) 
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To avoid the problems of finding solutions to the quadratic equation (3.45), we are now 
going to count the number of idealized traces instead of the number of actual traces. Any 
2xq G tr/F that does not occur in tr T is referred to as a gap in the length spectrum 
corresponding to T. It is the number of these gaps we have to make a hypothesis on. 

In order to count the number of idealized traces up to a certain value, the counting 
function 

Afj{R):^^.#{xoeM; \xo\ < R, |(^,(xo)| < 1, j = 2, . . . ,n} (3.48) 

will be introduced with R — cosh(//2). The factor of | takes care of the overcounting by 
admitting both signs for xq. 

The determination of the asymptotics of A//(-R) for R ^ oo will now be achieved by 
investigating the number of certain lattice points in some parallelotope. The procedure 
we are going to follow uses some standard receipt from algebraic number theory, see [34] 
and [65]. At first K is being mapped to K,^^ x . . . x K,y^ = iR" by: x e K, x t-^ x — 
{xi, . . . , Xn) '■— {(pi{x), . . . , (pn{x)). In FT' we consider, given n linearly independent vectors 
ei, . . . , e„, a lattice 

L := Zei ® • • • ® Zen (3.49) 

with fundamental cell 

F := /ei ® • • • ® /e„ , (3.50) 

/ := [0,1). In M"^ we shall consider usual euclidean volumes. F then has a volume of 
vol{F) = det{eij), where {cij) denotes the n x n matrix formed by the n row vectors 
ej e M^. We further introduce the parallelotope 

— {x e iR"; \xi\ < R, \xj\ < 1, j = 2, (3.51) 

of volume voI{Pr) ~ 2" R. In a first obvious approximation, the number nL{R) of lattice 
points in Pr is given by voI{Pr)/voI{F). Corrections to this result are caused by contri- 
butions of the surface of Pr; this is of dimension n — 1 (in M"^), whereas Pr itself is of 
dimension n. One therefore expects the corrections to be of the order of voI{PrY'^~^^/^ . 
Indeed, in [65] it is shown that 

n4i?) = ^^^^ + c ■ voI{Pr)'-'/- , (3.52) 

with some constant c. The surface correction is therefore subdominant in the limit it! — > oo 
and one finds that 

on 

One can now construct an appropriate lattice L that allows to represent J\fi{R) as the 
corresponding ^niiR)- To this end one notices that the module-basis of 
A4 G K is being mapped to a set of n linearly independent vectors {ei, . . . ,e„}, with 
ej :— {(fi{iij), . . . , (fni/J'j))- The independence may be seen by [det{eij)]'^ — [det{(fi{iij))Y — 
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DK/qiM) 7^ 0. One can thus use these e^'s to define a lattice L as in (3.49). Its fundamental 

cell F has volume vol(F) — ^ Dk/q{M.). An element Xq — /ci/^iH Vkn^n & M., kj E 

is being mapped to xq = kiei + ■ ■ • + kn^n G L, and this relation is clearly bijective. 
One can thus embed as L in iR", and hence (see (3.48)) Afi{R) = jnL{R)- Using 
R = cosh(//2) ~ ie'/^ / ^ oo, and (3.53), one concludes that 

M(cosh(V2)) ~ 2"-2 [Dk/q{M)]-'/^ e'/2 , I ^ oo . (3.54) 

As already mentioned, we are not able to pin down the exact number of gaps that might 
occur for a general Fuchsian group of the type F = I'iA, O). We expect, however, that the 
following hypothesis holds true: 

Hypothesis: Asymptotically, for / ^ oo, 

^(0 -M(cosh(//2)) . (3.55) 

In other words, it is assumed that the number of gaps grows at most like 0{e^^~^^^), 5 > 0, 
I —>■ oo. 

We are now in a position to state our result as the 

Theorem: Let F be an arithmetic Fuchsian group, commensurable with the group r(A, O) 
derived from the quaternion algebra A over the totally real algebraic number field K of 
degree n. Denote by di the index of the subgroup Fq := F fl V{A^O) in F, and by ^2 
the respective index of Fq in F(A, O). Let Dk/q{M.) be the discriminant of the module 
M. <Z K that contains F(/l, O). Then, under the hypothesis (3.55), the number A/^(/) 
of distinct primitive lengths on V\H up to / grows asymptotically like 

Sf^{l) ^ 2"-2 ^ [Dk/q{M)\-''^ ■ e"^ , l^oo. (3.56) 

Proof: Assume the validity of the hypothesis (3.55) and recall the asymptotic relation 
Np{l) ~ M{1), I oo, from section 3.2. Therefore also A/'p(/) ~ A//(cosh(//2)). Using 
(3.19) and (3.54) then leads to the assertion. 

The main objective of this study was not the counting function N'pil) but rather the 
mean multiplicity < gp{l) >. As these two quantities are asymptotically related by (3.13) 
one can easily obtain from the theorem the following 

Corollary: The local average of the primitive multiphcities in the cases described in the 
theorem behaves asymptotically like 



< ^,(0 >~ 2^-" ^ ^DKiai{M) ■ — , / ^ oo . (3.57) 

To confirm the results of the theorem and the corollary the two examples of arithmetic 
Fuchsian groups introduced at the end of the preceding section shall be investigated now. 
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Again, at first the modular group will be treated. For this O — M(2, and — 
SL{2, Z). Expand X E O into the basis (3.35), X = km + ■ ■ ■ + hn, h G ^, from 
which one observes that ^tr^X = xq = \{ki + k^). This yields M. = \'Z and /ii = |, 
see (3.43). It is known [66, 81] that for the modular group trV = = 2A4 and therefore 
no gaps in the set of traces occur. The discriminant of M. now is trivially obtained, and 



Dq(M) — lii — \. As the modular group Tmod is the quaternion group T{A,0) itself, 
one concludes using di = ^2 = 1, 



1/2 



eV2 



< gp{l) > ~ 2 • — , Z ^ oo , (3.58) 

which agrees with the previously known result (3.21). 

The second example, the regular octagon group Fj-eg, can almost as easily been dealt 
with. For this one has to go back to the quaternion group V{A, O) described by (3.37). One 

observes that xq E M. = ^[\/2], for which one can use the basis {111,112} = {1, a/2}. With 
the help of the conjugation (see (3.36)) one obtains the basis {61,62} for the lattice 
L C as {(1,1), (v^, ~v^)}- This allows to determine the discriminant of Ad, leading 



to Dk/q{M.) = \det{eij)\ = 2^/2. Since Treg is a subgroup of index two in T(A, O), one 
has to choose di — 1 and (i2 = 2 in (3.56) and (3.57). Thus 



4^ 



eV2 



eV2 



< ^p(0 > ~ 8^ • — > 1^00 . (3.59) 

This is exactly the result obtained in [10, 4]. In [10] the fact that \tr^ = m + n^/2, m,. n E 
Z,foij E r,.ej, was found for the first time. The condition that \ip2{^t^l) \ = \fn—n\f2\ < 1 
was then observed empirically. In [4] the numbers Xq = m + n\/2 fulfilling \m — n V2\ < 1 
were called minimal numbers and the necessity of this condition for the regular octagon 
group was shown. Also, by a numerical computation of the primitive length spectrum up 
to Z = 18 it was demonstrated that gaps do exist for the regular octagon, but that their 
existence does not influence the numerically calculated mean multiplicity < (?p(/) >. 

A final remark on a more constructive approach to determine the asymptotics of 
< 9p{l) > for ^ given arithmetic group F will be added. First, one has to know the 
quaternion group r{A, O) which is commensurable with F, and the indices di and ^2 de- 
scribing F n T[A, O) as a subgroup in F and T{A, O). In [94] it is shown that one can get 
the relevant number field K by adjoining trV^A, O) toQ. Furthermore, the algebra A can 
be obtained as the linear span of T[A, O) over K. Analogously, the linear span of T{A, O) 
over TZk yields the order O <Z A. One then has to find the module-basis {ri, . . . , r4„} of 
O. This can be used to obtain the matrix (Mjj) appearing in (3.40). Given this one has 
to identify the module M. containing F(A, O) (see (3.42) and (3.43)) in order to deter- 
mine its discriminant Dk/q{M.)- One can now plug all this information into (3.57) to get 
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the answer to the problem. This procedure may, however, be quite formal and in special 

cases it may be more convenient to try a direct approach to determine < gp{l) >■ Never- 
theless, the above result is general and sometimes it will only be necessary to know that 

I /2 

< fl'p(0 const. without specifying the constant. The latter, however, is included 
in the expressions (4.36) and (4.50) for the model describing the statistical properties of 
the related quantum energy spectra. For a quantitative description a knowledge of the 
constants appears to be necessary. 



3.5 Pseudosymmetries 

Up to now those classical aspects of arithmetical chaos have been investigated that are 
connected with geodesic length spectra. There is a further peculiarity of the arithmetical 
systems, which is a property of the classical system also appearing in its quantum version. 
The phenomenon to be discussed in this section is the occurrence of infinitely many pseu- 
dosymmetries for hyperbolic surfaces r\H with arithmetic Fuchsian groups F. As these are 
closely related to the Hecke ring for F, which will be represented on the quantum mechan- 
ical wave functions by Hecke operators, the pseudosymmetries somehow mediate between 
the classical and quantum aspects of arithmetical chaos. Before going into the details, the 
description of symmetries of hyperbolic surfaces will be briefly reviewed in order to get an 
intuitive understanding of pseudosymmetries as some generalizations of symmetries. 

A symmetry g oi a hyperbolic surface r\H is an isometry of this surface, and therefore 
necessarily also an isometry of the hyperbolic plane Ti. Thus g is a fractional linear 
transformation on Ti. The corresponding matrix from SL{2, M) will also be denoted by 
g. The symmetry group S = {11, gi, . . . , g^^i} of the surface F\7i is hence a subgroup 
oi SL{2, IR). In order that & g e SL(2,1R) is a symmetry it has to commute with the 
Fuchsian group F. This may be seen as follows: z E H is being identified with 72; for 
all 7 G F, thus also gz with jgz. If 5? is a symmetry, then the identification of z with 
g'^'-fgz has to be the same as of z with 7' 2; for all 7,7' G F. This will be the case, if 
(and only if) g~^Tg = F. Defining the group F' := F U giT U . . . U gN-i^ this condition 
means that F is a normal subgroup of F', and E = F'/F. The fact that the surface r\H 
possesses symmetries can thus be formulated in an algebraic way. The Fuchsian group 
F is a normal subgroup of some other Fuchsian group F' and the symmetry group S is 
the factor group F'/F. This algebraic setting then allows to deal with symmetries in the 
context of the Selberg trace formula, see e.g. [100, 101] and appendix B. Since F' is also 
a Fuchsian group one can construct the surface F'\7^. This can be viewed as the result 
of a desymmetrization procedure. It is, loosely speaking, a fundamental domain for the 
operation of the symmetry group E on the surface F\7i. 

A certain generalization of this concept of symmetries is provided by the pseudosym- 
metries of arithmetic surfaces. The starting point of their discussion will be the algebraic 
properties satisfied by arithmetic Fuchsian groups. From this the geometric properties of 
the related surfaces will be studied. The algebraic side of the problem lies at the heart of 
the construction of arithmetic groups and may be found e.g. in [86, 100, 71]. It resulted 
from a generalization of Hecke's investigation [49] of automorphic forms for the modular 
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group. This is the reason why in number theory the notion of modular correspondences 
has been introduced for the algebraic setting. To the author's knowledge, it was Sarnak 
who recently introduced [79] the name pseudo symmetries for the geometric setting related 
to the modular correspondences. Since the geometric construction may be easier to un- 
derstand the problem intuitively, henceforth the notion of pseudosymmetries will be used 
throughout. It is the purpose of this section to explain the algebraic as well as the geo- 
metric structures accompanying the pseudosymmetries as explicitly as it seems possible. 
This, however, inevitably requires some algebraic notions that shall be introduced first. 

A major role will be played in the following by the commensurator F of a Fuchsian 
group F. It is a subgroup of G := GL+{2,R) = {g E GL{2,1R); detg > 0}. Then 
r :— {g E G; g^^Tg ~ F}, where "~" denotes the commensurability of two subgroups 
of a group. A g E r hence transforms the Fuchsian group F by conjugation into a group 
g~^Tg that is still commensurable with the original group F. Thus T'{g) := g~^Tg fl F 
is a subgroup of finite index in F and in g~^Tg. F clearly is a group, and it obviously 
contains F as a subgroup. If is a symmetry of the surface r\H, then g~^rg = F, leading 
to T'{g) — F, and thus g eT. Therefore, F contains F' (see above) as a subgroup. We will 
speak of a non-trivial commensurator, if F contains F' as a proper subgroup. The objects 
of F of interest are then the cosets of F'\F. We will henceforth tacitly assume to mean 
representatives of these cosets when speaking of the commensurator. These will give rise 
to pseudosymmetries. 

A commensurator F is defined for any Fuchsian group F and no use has so far been 
made of the arithmeticity of the groups of interest. The difference between commensura- 
tors of arithmetic and non- arithmetic groups is clarified by a theorem of Margulis [67, 93] : 

Theorem: If F is an arithmetic Fuchsian group, then its commensurator F is dense in G. 
If F is non-arithmetic, then F is commensurable with F. 

Remark: If F is commensurable with F, then FflF is of finite index in F. Now, FflF = F, 
since F C F. Thus F is a subgroup of finite index o? in F, F = F U F71 U . . . U T'jd-i- The set 
of non-trivial elements of the commensurator hence is {1, 71, ... , 7d-i}, which is finite, and 
therefore at most finitely many pseudosymmetries exist for non-arithmetic groups. Since 
in the arithmetic case F is dense in G infinitely many pseudosymmetries are then present. 
This criterion thus may serve as a characterization of arithmeticity of Fuchsian groups. 
Although in the non-arithmetic case finitely many pseudosymmetries might exist, to the 
author's knowledge no explicit example of a non-trivial pseudosymmetry is known for this 
case. 

By the definition of the commensurator, for every g E T the group r'{g) is a subgroup 
of finite index n in F, thus 

F = F'(5)7iU...UF'(5)7„ , 71 = 1. (3.60) 

The pseudosymmetry related to this g is then said to be of order n. A symmetry is in 
this notation a pseudosymmetry of order n — 1, since if g E 12, then F = F'(g'). For the 
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following it appears to be useful not to deal with g' G F, but rather with the double cosets 
TgT. By (3.60) these decompose as 

r^r = rg7i u . . . u rg7, = r^i u . . . u ra^ , (3.61) 

where the definition := g'-fi has been used. The double cosets TgT for G F are 
the quantities to construct the Hecke ring for F from. Its definition and discussion will 
be postponed to section 4.1 where the quantum aspects of arithmetical chaos will be 
investigated, since it directly leads to the definition of Hecke operators that are relevant 
for the quantum mechanical problem. 

The decomposition (3.60) of F into cosets of T'{g) now enables one to interpret the 
pseudosymmetry related to g geometrically We require g to he a pseudosymmctry of 
order n > 2 in order to deal with a proper generalization of a symmetry. Then T'{g) is 
a proper subgroup of F, and the surface r'{g)\H is an n-sheeted covering of the original 
surface r\H. To illustrate the situation one can draw the following commutative diagram 
(see e.g. [40, 86]): 

id 

n > n 

T\g)\n — ^ T\H 



if, (pg and / are the natural projections. For z ^ Ti denote by F^ := {7 G F; 72; = z} 
the subgroup of F that stabilizes z; the corresponding subgroup of T'{g) then is F'^ := 
TzHT'^g). Tz and F'^ consist of elliptic elements (or parabolic elements if one admits cusps 
as fixed points). For z E H denote its image under the projection on the surface F\7i 
by p = (p{z). The preimage of p under / on r'{g)\H then consists of the h < n points 
f~^{p) = {^15 • • • ) Qh}- If denotes the ramification number of / over qj, then Z]j=i = 
On H one then chooses points Wi,...,Wh such that qj = ipg{wj), i.e. p = f{ipg{wj)), 
j = 1, . . . ,h. Ramification numbers ej different from one can only occur at elliptic points of 
F\7i , since ej is the index of T'^. as a subgroup of Tyj. . These two groups are non-trivial only 
at elliptic points. The fractional linear transformations {cti, .... a^}, defined by ajz = Wj, 
mediate the mappings that interchange the sheets of the covering r'{g)\H F\7i, when 
projected down onto the surfaces by (fg and if. Since the above diagram is commutative 
one finds for all 7 G F that f~^{p) = f~^{ip{z)) = f~^{ip{'yz)) = ipgi^jz). Thus there exists 
a unique index j such that (pgijz) = qj] hence ipgl'jz) = qj = ipg{wj) = (pg{ajz). Since 72; 
and ajZ project to the same point on r'{g)\H, there is some S G T'{g) with 72; = SajZ. 
Therefore '~f~^5aj G F^, implying by inverting the l.h.s., that 7 G F'(gf)crjF;j. One hence 
obtains that 

F = T'{g)a^Tz U...U T'ig)a^,Tz . (3.62) 

It is not difficult to show that this decomposition of F is disjoint. Choosing z not to be 
an elliptic point, hence F;^ = {l} and h = n, and comparing (3.62) with (3.60) leads to 
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the conclusion that Gj — ^'fij-, for some 7^- e ^'{9)-, j — l,---,n. The transformations 
{71, . . . ,7n} appearing in (3.60) can therefore be given an interpretation as mediating the 
interchanging of the sheets of the covering outside branch points. If T is strictly hyperbolic, 
i.e. if T\H is a compact surface of genus g > 2 without elliptic points, then = {1} for 
ell z E H and hence all ramification numbers are one. In this case r'{g)\H — > T\H is an 
unramified n-sheeted covering. 

In conclusion one can give the following geometric picture of a non-trivial pseudosym- 
metry of order n > 2 (see also [79]): it leads to an n-sheeted covering of the surface r\H 
that is unramified, if F contains no elliptic elements. Otherwise it is ramified over the 
elliptic fixed points of F. The stabihzing group F^ of an elhptic fixed point z then deter- 
mines the ramification number at this point. In case the operation under consideration is a 
symmetry, it is a pseudosymmetry of order n = 1, and thus no non-trivial covering occurs. 
The generalization of symmetries to non-trivial pseudosymmetries in this picture consists 
in the spreading out of the coverings over the base surface. But still, by definition, these 
always have finite numbers of sheets. 

Another interpretation of a pseudosymmetry related to a non-trivial g & r can be given 
in terms of the effect of g on the closed geodesies on F\7i. The operation oi g on Ti is 
equivalent to its operation on F by conjugation. Thus its effect on a closed geodesic can 
be described by the mapping 7 1— > g~^jg for the hyperbolic 7 e F related to the geodesic. 
Since g~^Tg is commensurable with F, the fraction of 7 e F that are being mapped onto 
F is given by the index of F'(gr) = g~^^g fl F in F. By (3.60) this is just the order n of the 
pseudosymmetry. Thus a finite fraction of ^ of the closed geodesies on r\Tl are mapped 
by g again to closed geodesies on the same surface. Notice that this reasoning is similar to 
the one leading to (3.17). 

We are now going to consider the modular group Tmod = SL{2, Z) as an example to 
illustrate the above construction explicitly. As already mentioned, this was historically also 
the first case to be studied, and where modular correspondences have been introduced. A 
nice presentation of several facts about the modular group can be found in [96]; also [86, 71] 
are useful to be consulted. The first task in this context is to obtain the commensurator 
F of the modular group. This can be found to be F = GL+(2,Q), a fact that is relatively 
easy to prove. (See [86, 71] for details.) Notice that GL+(2,(Q) is dense in GL+(2,1?), as 
it is predicted by Margulis' theorem. Let now be Af„(^) := {g G M(2, ^); det g — n}, 
n G IV. This set is not a group, as the multiplication does not close in it. However, it may 
be decomposed disjointly as (see [96]) 

M„(^) = U F^„, f d ) • ^^-^^^ 

ad=n \ / 

0<b<d 

Now define the semi-group 

A := G M(2, Z); detg>0}= [j M,,{Z) , (3.64) 

new 

and take a G F. Then there exists a g G Q such that g' := qg G A. Since g~^Tg = g'~^Vg' , 
ones attention may be restricted from F to A. In [71] one can now find the following 
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result: let g & A, then there exist l,m ^ IN, /|m (this notation means that I is a divisor 
of m), such that rmodg^mod = ^mod[o^)^mod- It is then further shown in [71] that 
^mod (o V) ^rnod = U Fmod (q ^ , orf = Zm, < 6 < d, (a, b, d) = I. (The latter notation 
means that the largest common divisor of a, h and d is /.) Thus the {ai\ in (3.61) are 
given by {(q ^) ; ad = Im, < b < d, {a, b, d) = I}. Conventionally one looks at all e A 
with det g = n, n & ]N, simultaneously and finds that 



U ^modd^mod — Mn{Z) — U ^mod ( 

gGA ad=n \ 

det g=n 0<.b<d 




(3.65) 




The simplest non-trivial case is n = 2, where one has to choose g — (^^^^ with I — 1, 
m = 2. Then 

Tmod = Tmod ^ q 2 ^ ^mod ^ q 2 ^ ^ ^mod ^ q 1 ^ " (3.66) 

Since the r.h.s. of (3.66) consists of three cosets, the pseudosymmetry related to this g G 
A is of order three. The corresponding covering r'{g)\H — > T,nod\H will be found by 
studying the effect of g on an arbitrary 7 e Tj^od- Let therefore be G Tmod- Then 

s-'{:'d)9=iyi).Thus 

^'{9) = g'^^modg n Tmod = | ( c ) ^ ^' ^^^^} ^^"^^^ 

is the group that defines the desired covering over the modular surface Tmod\H. 

The modular group certainly provides the simplest example to consider. In [86, 71] also 
congruence subgroups of Tmod are treated explicitly. Miyake even goes one step further in 
that he deals with arithmetic groups that are unit groups of quaternion algebras over (Q. 
Already then he has to argue adelicly, which complicates the discussion considerably. To 
the author's knowledge there does not exist any explicit treatment of arithmetic groups 
derived from quaternion algebras over number fields K of degree n > 2. In such cases a 
first major obstacle is to identify the commensurator group F. An explicit knowledge of 
the relations (3.60) and (3.61) for a given arithmetic group is mandatory to construct the 
Hecke ring and the Hecke operators for that group explicitly. From the case of the modular 
group one learned that this knowledge helped a lot for the numerical determination of the 
energy eigenvalues and the eigenf unctions, see [91]. To apply the methods used for the 
modular group also to e.g. the regular octagon group, however, seems at the moment too 
hard a problem. 



3.6 Statistical Properties of Geodesic Length Spectra 

This final section of the present chapter contains a topic that might serve as a link to 
the following chapter dealing with the quantum mechanics of arithmetical chaos. The 
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reason for this is that the statistical properties of length spectra on hyperbolic surfaces are 
investigated using the Selberg trace formula and the Sclbcrg zeta function, which provide 
a duality relation between classical and quantum aspects. It is then possible to use results 
on quantum energy spectra to gain information on classical length spectra. Also, some 
analytic methods employed here using properties of the Selberg zeta function are quite 
similar to those that will be applied in chapter 4. 

In the same way as the spectral staircase N{E) plays a major role in the investigation of 
quantum energy spectra, the classical counting function A/^(/) is an important tool to study 
geodesic length spectra. Since the PGT (3.8) states that the respective counting functions 
J\fp{l) share the same asymptotic behaviour for all hyperbolic surfaces, it is the remainder 
to the asymptotic value that is the quantity of interest. It plays a similar role as Nfi{E) 
does for quantum energy spectra. The aim of this section now is to gain information on this 
remainder; and by what has been said so far about geodesic length spectra of hyperbolic 
surfaces it seems to be clear that arithmetic Fuchsian groups have to be distinguished from 
non- arithmetic ones. For the analytic part of the study an analogy to the Riemann zeta 
function appears to be constructive, since an important tool to obtain information 
on the remainder term to the PGT will be the Selberg zeta function Z{s) [84, 51, 100], 
which is in many respects similar to the Riemann zeta function. The theory of the latter is 
briefly reviewed in appendix A. A detailed analysis of C,{s) then leads to the prime number 
theorem (PNT) [97, 59], which states that the number of primes n{x) not exceeding the 
value X is asymptotically given by t:{x) ~ li{x) ~ x — > oo. Identifying the n-th prime 
Pn with e'" then clearly shows a close resemblance of the PNT with the PGT (3.8). It shall 
now be explained how far the analogy between primes and primitive closed geodesies on a 
hyperbolic surface can go. 

As already mentioned, the Selberg zeta function Z{s) is needed in order to proof the 
PGT and to estimate the remainder term. The notation already indicates that the Selberg 
zeta function is the dynamical zeta function (2.9) for geodesic flows on hyperbolic sur- 
faces. In the same way as the dynamical zeta function is derived from Gutzwiller's trace 
formula with a special test function (2.7), the Selberg zeta function is obtained from the 
Selberg trace formula [84, 51, 100]. The latter is an exact analogue of the smeared version 
(1.4) of Gutzwiller's trace formula, i.e. it is an exact identity and not only a semiclassical 
approximation. 

For ease of notation from now on only strictly hyperbolic Fuchsian groups will be con- 
sidered in this section. The general case of Fuchsian groups of the first kind yields nothing 
new regarding the present problem, since the contributions of the hyperbolic elements are 
the relevant ones concerning length spectra. But of course, elliptic and parabolic elements 
can also be treated, see [84, 51, 100] for details. 

The quantization of the geodesic flow on a hyperbolic surface V\H is determined by 
the stationary Schrodinger equation 

-AV'(^) = Ei^iz) . (3.68) 



41 



The hyperbolic Laplacian is given in terms of the coordinates of Ti. by 



A = y\dl + dl), (3.69) 

and the wave functions are required to be invariant under the operation of F on Ti, ipi'jz) = 
i/j{z) for all 7 G r, in order to yield functions on the orbit space r\H. —A can then be 
defined as a self-adjoint operator on L^(r\H) with a purely discrete spectrum = < 
El < E2 < . . ., En = p'^ + \. The scalar product on L'^{T\H) is derived from the hyperbolic 
metric as 

KA'py^ f ^W)^(^), (3.70) 

Jr\n 

for ip^ip E L?'(V\H). Quantum energies Ek that arc related to complex momenta p^, 
Q < Ek < J, are called small eigenvalues. Their existence depends on the geometry of 
r\7i, and it is known that only finitely many can exist on a single surface, see [36] for a 
review. On compact surfaces of genus 51 = 2 at most one might occur [83] . Small eigenvalues 
play a special role in connection with the PGT and have to be treated separately. 
The Selberg trace formula now reads [84, 51, 100] 

, area(jr) /•+00 ^ Kl) g{kl{-i)) 

E'^'^-^'^L ''''''''<^>*°°'''-''>%|:g„ 2sVnh?M(T)/2) - '^-"^ 

The outer sum on the r.h.s. of (3.71) runs over all F-conjugacy classes of primitive hy- 
perbolic 7 e r, thus equivalently, over all primitive closed geodesies on T\H. Comparing 
(3.71) with Gutzwiller's trace formula (1.4) shows that all Lyapunov exponents are identi- 
cal, = 1 for all 7, and hence are also identical to the metric entropy A = 1. An admissible 
test function h{p) has to be even, and to be holomorphic in the strip \Imp\ < | + e, e > 0. 
Also, h{p) = 0(|p|~^^^) for \p\ — > 00. g(x) = ^e*^^/i(p) denotes the Fourier-transform 
of the test function. The Selberg zeta function Z{s) arises from (3.71) as the dynamical 
zeta function has been obtained from (2.7), see (2.7)-(2.9). Its Euler product converges 
hr Res > T — 1, 

00 

^(«) = n n (1 - e-(^+")'(^)) , (3.72) 

where the variable s = | — ip is related to the energy variable through E — s{l — s). Choos- 
ing the test function h{p) — p2_^_(l_k^2 ~ p2_)_(j^_i)2 foi" Res, Re a > 1 yields a regularized 
trace of the resolvent operator for —A, which is the appropriate analogue of (2.7). From 
this relation one can obtain the analytic properties of Z{s). It is an entire holomorphic 
function with trivial zeros at s=0 of multiplicity ^^^^^^^ + 1, at s = 1 of multiplicity one, 
and at s = —k, k e W, of multiphcities ^^^^^^(/i; + 1). Its non-trivial zeros are related to 
the eigenvalues of —A through s„ = | it En = Sn(l — Sn); their multiplicities are given 
by the respective multiplicities of the eigenvalues. Small eigenvalues therefore correspond 
to zeros of Z{s) in the interval (0, 1). Leaving aside the latter ones Z(s) thus fulfills an 
analogue of the Riemann hypothesis (RH) for ({s) in that ReSn — 
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In the case of the Riemann zeta function the magnitude of the remainder term in the 
PNT is determined by the non-trivial zero with largest real part ctq, see appendix A. Since 
for the Selberg zeta function the RH is known to be true, and hence (Tq = |, once the 
contributions of small eigenvalues have been extracted explicitly, one would expect the 
remainder term Qr{1) in the PGT (see also (3.85)) to grow hke 

\QR{l)\ = e-^'-iuil) , (3.73) 

where uj{l) is a combination of powers and logarithms of /. Up to now it was, however, 
not possible to prove this. The analogy with ({s) will in the following be pushed as far as 
possible, adopting the strategy employed in appendix A for ({s) to construct a Dirichlet 
series whose abscissa of conditional convergence should act as a "detector" for o-q. It will 
then become clear where the obstacle comes from that prevents one from proving the 
assertion (3.73) about the magnitude of the remainder term. 

The Selberg zeta function has a simple zero at s = 1 due to the eigenvalue Eq = 0. 
Hence in the vicinity of s = 1 it behaves like Z{s) = Z'{1) (s — 1) + 0((s — 1)^). One 
could now think of using a Dirichlet series for Z{s) itself, since it has no pole at s = 1, 
as C(s) has, and thus there seems to be no obvious obstruction in pushing the domain 
of convergence of the Dirichlet series further to the left in the s-planc. However, Z{s) is 
entire holomorphic so that even on the critical line Re s = | there are no poles serving 
as such obstructions. Therefore the use of Z{s) would not yield the desired effect. To 
circumvent this problem one would like to use Z{s)~^ instead, which has the desired poles 
at the non-trivial zeros of Z{s), and then subtract the poles at s = 1 and at the e (|, 1), 
k = 1,...,M, corresponding to small eigenvalues. The problem one then immediately 
faces is that because of the product over n in the Euler product representation (3.72) of 
Z{s) there is no convenient Dirichlet series for the inverse of the Selberg zeta function. The 
way out of this problem will be to discard the n-product and to define a new Ruelle-type 
zeta function R{s) :— which has for Res > 1 the Euler product representation 

i?(s) = n(l-e-^'(^)) . (3.74) 

{7}p 

This is easily obtained by inserting (3.72) into the definition of R{s). (3.74) now is the 
exact analogue to the inverse of the Euler product for C{s). Using the analytic properties 
of Z{s) one can derive those of R{s). The latter is a meromorphic function of s G (D 
that is holomorphic for Res > 0. In this half-plane it has the same zeros as Z{s). Small 
eigenvalues lead to zeros at Si, . . . , sm in (|, 1) and at 1 — Si, . . . ,1 — sm in (0, |); the 
eigenvalue Eq = produces a zero at Sq = 1- If the small eigenvalues are not degenerate, 
R{s) behaves in the vicinity of Sk, k — 0, . . . , M, like R{s) — z{s^^+i) ~ Sk)+0{{s — s^)^). 
Its logarithmic derivative hence has a simple pole with residue one at the s^s. In order to 
obtain a meromorphic function with poles at the non-trivial zeros of Z[s) on the critical 
line that is holomorphic for Re s > | one can define 



43 



To simplify the notation it will henceforth be assumed that there do not occur small 
eigenvalues, hence M — 0. The only pole that has to be subtracted then is the one at 
So = 1. 

The next step now consists of finding a Dirichlet series representation for From 
the Euler product (3.74) for R{s) one finds that for Res > 1 

H^-E'Wy^^EE/We-*^'. (3.76) 

In view of Beurling's theory of generalized prime numbers (see e.g. [18]) one can identify 
primitive hyperbolic conjugacy classes {7}^ with primes p [15]. Comparing (3.76) with the 
Dirichlet series for — then leads to the definition of an analogue of the von Mangoldt 
function A(n), see appendix A. What is lacking so far is an analogue of the positive integers 
IV. Recently it has become customary to introduce such analogues as pseudo- orbits [69, 22]. 
These are the generalized integers in Beurling's theory and comprise of formal combinations 
of powers of primitive conjugacy classes, 

p:={7Mp®---®{7Mp- (3.77) 

On the surface r\H p corresponds to a formal combination of the primitive closed geodesies 
related to 71, . . . , 7^ that are traversed ki, . . . ,kn times, respectively. These are the objects 
that were named pseudo-orbits [22]. The corresponding "lengths" 

Lp = fci/(7i) + ■ ■ ■ + A;„/(7„) (3.78) 

are then called pseudo-lengths. Having defined analogues of integers it is now possible to 
introduce a von Mangoldt function for R{s), 

, otherwise " ^^'^^^ 

With the help of the above notions one can first introduce a Dirichlet series for i?(s)"^, 

-1 00 

(3.80) 

^y^j {7}j, fe=o p 

and then one for /r(s), 

fnis) = E e-^^ , Ap = AM - 1^ . (3.81) 

In appendix A the Chebyshev functions 9{x) and ip{x) were used to relate the analytic 
properties of ({s) to the PNT. Trying to carry this over to the case of the Ruelle-type zeta 
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function and the PGT one is led to define analogues of the Chebyshev functions as 



ML) := EA«(P) = E E Kl) (3.82) 

Lp<L k>l kl{j)<L 

k>\ 

(In the above notations the /(7)'s and the Lp's are counted with their respective multi- 
plicities.) As with the classical Chebyshev functions, see appendix A, an estimate of the 
remainder in iPr{L) = 9r{L) + TZr{L) gives 1Zr{L) = 0{e2^L'^), L ^ oo. Since the 
analogue of the prime counting function n{x) is the counting function J\fp{l) for primitive 
closed geodesies, exactly the same reasoning as in appendix A leads to 

m) -^ + 1% Ml') + 0{e¥l) , (3.83) 

where li denotes the length of the shortest closed geodesic on r\H. Using the integral 
(A. 11) then allows to express il^R^L) through by (6 > 1, L 7^ Lp) 

ML) = E Mp) ^ / - = 7^ - ■ 3.84 

p 2m Jb-ioo s 2m Jb-ioc s R{s) 

Since the analytic properties of the integrand on the very right of (3.84) are known one 
would expect to obtain an analogue of the explicit formula (A. 13) for the classical Cheby- 
shev function ip{x). This would involve a sum J2sn ^^^^ non-trivial zeros of Z{s) on 
the critical line. At this point, however, the present case differs from that of the Riemann 

zeta function in that the sum over the s„'s diverges and thus no explicit formula in the 
desired manner exists for ipR^L). The reason for this difference lies in the stronger growth 
of N{E) as compared to the counting function iV^(p) := {s„ = Pn+iln] < 7^ < p}- From 
(3.71) one can rederive Weyl's law (1.3) to yield N{E{p)) ~ "' ^^^'^'' p^, p — > 00, whereas it 
is known [97] that Ni;{p) ~ ^logp, p — > 00. There are hence "too many" terms in the 
sum over the s„'s that prevent it from converging. 

Hejhal [51] proceeds in defining an integrated version 'ijjR^i{L) := dl eJip r{1) of the 
Chebyshev function iIjr{L). Inspecting (3.84) one observes that one has to deal with the 
sum J2s„ s ts"+i) instead, which is conditionally convergent. A tedious analysis then leads 
to the PGT {Theorem 6.19m [51]) 

M 

ML) = + E e^'^^ + ^ii(^) , PR(L)=0(et^L^), 
fe=i 

M 

J\fp{l) = Ei{l) + E E^Skl) + Qr{1) , Qr{1) = 0(ei'/-5) , (3.85) 

k=l 
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where the contributions of small eigenvalues have been reintroduced explicitly. 

Another way to obtain the PGT is to restrict the contour of integration in (3.84) to 
a finite interval and then to estimate the resulting sum. Employing methods that can be 
found in [97], pp. 60, and in [60] one can show that (6 > 1, < T < oo, L 7^ Lp) 



ML) 



1 r^+iT ds 



R'i 



27ri Jb-iT s 



R{s) 



+ 0(L2T-ie^) + 0{T-\h - l)-^e^^) 



(3.86) 



The proof for the remainder terms on the r.h.s. is the same as for the lemma of section 
4.5, where Z{s) replaces We therefore postpone its explicit discussion to section 4.5. 

Deforming the contour of integration in (3.86) and using estimates for on the contour 
Iwaniec succeeded to estimate Qr{1) for the modular group as [60] 



QR,moci{l) = 0(e(i+^)') V£ > , 



(3.87) 



which is slightly better than the general result (3.85). This is also the best upper bound 
for the remainder term in a PGT available. It seems that it is very hard to "break the 
barrier" at e*' for the upper bound of Qr{1). 

This difficulty has a consequence for estimating the abscissa of conditional convergence 
for the Dirichlet series representation (3.81) of /r(s). Since /r(s) has been designed in 
complete analogy to the function /(s) introduced in appendix A to study the PNT, the 
analysis applied to the latter function will now be repeated for the former one. In appendix 
A also the convergence properties of general Dirichlet series have been reviewed, according 
to which the series (3.81) converges for Res > ac and converges absolutely for Res > aa, 
Ca > CTc- Arranging the pseudo-orbits in ascending order of their lengths, Li < L2 < L3 < 
. . ., one finds that 



1 



N 



aa = limsup -j—^ogY^ 



Afl(ri) - 



Z'{1) 



Z{2) 



(Tc — lim sup - — log 
— lim sup - — log 



^ ^ Z'{1)\ 



n=l 



(3.88) 



Z'il] 



Z(2) 



where J\f^^\L) := Lp < L} denotes the counting function for pseudo-orbits. In [15] 

it was shown that 



(3.89) 



with some constants Ci > and < o; < |. As for f{s) (see appendix A) one also 
concludes here, using (3.89), that E^=i \^R{n) - ~ const. ■ e^^ . The Dirichlet series 
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(3.81) therefore converges absolutely ior Re s > aa — 1. Inserting '^r{L) — + Pr{L) one 
observes that 



lim sup - — log 



(3.90) 



The only conclusion one can now draw from (3.90), and from the PGT (3.85), Pr{L) = 
0(e4^L2), is that Uc < 1. The analytic properties of /i?(s), however, suggest that ctc — ^ 
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since /r(s) is a holomorphic function for Res > and has poles on the critical line. The 
analogous function /(s) built from the Riemann zeta function indeed shows this behaviour: 
its Dirichlet series converges (conditionally) on the maximal half-plane Res > ctq where 
f{s) is still holomorphic. 

The weakness of the upper bound Uc < 1 hinges on the estimate Q^^\L) = 0{e^~'^^^"), 
which was obtained in [15] from the theory of generalized prime numbers, see e.g. [18]. 
However, this upper bound only requires the rather weak estimate Qr{1) = 0(6'"*^") for 
some constants c > and < f3 < 1. Since Qr{1) is known to be smaller than required, 
see (3.85), the true magnitude of Q^^\L) could be much smaller than just being below 
e^. To support this idea a function similar to /r(s) will be introduced. Starting with the 
Dirichlet series (3.81) for R{s)~^ one obtains for i?es > 1 through an integration by parts 



-sL 



J_ = dAA(^) (L) e-^^ = AA(^) (L) e-^^ ^ + s T dL AA^^) (L) e 
R{s) Jo Jo 



/•CO 

+ s dL Q^^\L) e-'^ , (3.92) 
Jo 



This yields an integral representation 

1 Z{2) 1 _ Z{2) 
R(s) ~ Z'{1) s - 1 ~ Z'{1) 

of a function that is holomorphic for Res > \ and has poles on the critical line. The 
estimate for Q^^\L), however, only permits to use the r.h.s. of (3.92) for it!es > 1. Again, 
the domain of holomorphy of the l.h.s. suggests that the integral might exist for i?es > |, 
yielding the estimate Q^^\L) = 0{e^^). Assuming now that indeed (Tc = \ for the Dirichlet 
series of fR{s), then by (3.90) one concludes that Pr{L) = 0(e*^5+'^)^) Ve > 0, leading to 
the estimate Qr{1) = 0(e*^2+'^)') > 0. The desired result (3.73) thus would follow if 
one assumed that the two representations (3.81) and (3.92) converged in the maximal 
half-planes where the functions they define are still holomorphic. 

Certainly, (3.92) gives the lower bound Q^^\L) — fi(e2^). This estimate already 
accounts for the lower bound (Tc > |, that arises from the fact that /r(<s) has poles 
on the critical line, by inserting Q^^\L) into (3.90). For Pr{L) Hcjhal could prove the 
lower bound Pr{L) = r2-|-(e2^(logL) 2 ) [51], implying for the remainder term in the PGT 
Qr{1) = r2±(e2'/^^(log/) 2). But unfortunately no further rigorous conclusion can be drawn 
from (3.90) and (3.92). The analytic theory seems to be stuck at this point. The question 
one might ask now is whether the problem is a technical one or whether there lurks some 
yet undiscovered phenomenon behind it. 
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One idea that could come to one's mind is that \Qr{1)\ — e^''u;{l) is indeed true for 
generic, i.e. non-arithmetic, Fuchsian groups. The arithmetic case should then be treated 
separately and might violate the expected behaviour of Qr{1). It seems to be quite natural 
to distinguish arithmetic from non-arithmetic groups, after having discussed the excep- 
tional structure of length spectra in the arithmetic case. Np{l) is a staircase function with 
steps of width AZ„ = lp,n+\ ~ lp,n a-^d of height gp{lp,n) I — lp,n- It is thus the inter- 
play of fluctuations of lengths and multiplicities that results in fluctuations of the staircase 
function Mp{l). This is in contrast to the fluctuation properties of the spectral staircase 
N{E), for which only fluctuations in the quantum energies En are responsible (since we 
require the systems under consideration to be completely desymmetrized and thus being 
void of degeneracies in their energy spectra). The interferences of the two contributions 
to Qr{1) describing the fluctuations of Np{l) are involved, but different multiplicities of 
lengths in the arithmetic and the non- arithmetic cases clearly lead to different kinds of 
fluctuations. For arithmetic groups the mean step heights are < gp{l) >~ c-t^j--, I oo. 
These alone give a contribution of the order /"^e'/^ to Qr{1), since the mean behaviour 
Ei{l) cannot follow the step structure of N'p{l). In addition to its mean behaviour fluc- 
tuations of gp{l) can give further contributions to Qr{1), let alone the fluctuations of the 
lengths themselves. For non-arithmetic groups the mean multiplicities, i.e. the mean step 
heights, do not give an exponential contribution to Qr(1). Since in any case the lower 
bound Qr{1) = f2±(e2'/~2 (log/) 2) requires exponentially large oscillations about Ei{l) (in 
the positive and the negative direction), in the non-arithmetic case these have to come 
from length fluctuations. 

In the following arithmetic and non-arithmetic Fuchsian groups will therefore be treated 
separately. At first for the generic, non-arithmetic case an approach of inverse quantum 
chaology will be employed. The latter notion stands for drawing conclusions on the classi- 
cal properties of a chaotic system from its quantum energy spectrum using the (Selberg or 
Gutzwiller) trace formula. As mentioned in chapter 2 Berry's theory of spectral rigidity pre- 
dicts for the energy dependence of the saturation value Aoo(-E') of the rigidity a logarithmic 
behaviour, Aoo{E) ~ ^ logE, E — ^ 00, if the classical system is chaotic and time-reversal 
invariant. It has also been mentioned in chapter 2 that this result cannot be applied to 
arithmetic systems because of their exponentially degenerate length spectra. Since Aoo{E) 
is related to Nfi{E) via (2.24), one can conclude that \Nfi{E)\ ~ -^^/logE, E ^ 00. 
Hejhal now proved a theorem that yields an estimate for the remainder term in the PGT 
depending on an upper bound for Nfi[E), therefore being truely a result in the spirit of in- 
verse quantum chaology. In the mathematical literature on the Selberg trace formula it has 
become customary to introduce the notation S{p) :— Nfi{E{p)) = ^ arg Z{^+ip). Hejhal's 
Theorem 14-18 [51] then states that for < 5 < 00 the estimate \S{p)\ — 0((logp)^) im- 
plies Pr{L) — 0{e2^L'^~^^). In our case, 5 = |, the upper bound Qr{1) — 0[e^H^) 
follows. Comparing this with the lower bound Qr{1) — 0±(e^'/~^(log/)^) one concludes 

that |Qii(OI — where the asymptotics of lie somewhere between and 

3 

/2. It should be stressed that this conclusion has not been drawn rigorously because it 
hinges on Berry's non-rigorous theory for A.^[L;E). It is not clear either what effect the 



48 



unboundedness of the multiplicities in the length spectra even for non-arithmetic groups 
will have on the applicability of Berry's theory. However, together with everything else 
discussed above, this gives a further hint supporting the expectation that the remainder 
term in the PGT grows like (3.73) for generic Fuchsian groups. 

Considering arithmetic groups the inverse quantum chaology reasoning does not yield 
the same result as in the generic case since the upper bound for S{p) assumed in Theorem 
14-18 of [51] is violated here. For a certain class of arithmetic groups (derived from quater- 
nion algebras defined over Q) Hejhal obtains in his Theorem 18.8 and Remark 18.14 of [51] 
the lower bound S{p) = logp). (Originally this was an unpublished result of Sel- 

berg.) Since the proof only requires the exponential increase of the multiplicities of lengths 
present for all arithmetic groups, the Q-result for S{p) applies to any arithmetic group. 
This lower bound is quite close to the general upper bound S{p) = 0{p/ logp), see [51]. In 
order to deal with the arithmetic case one then has to use Hejhal's Theorem 15. 13 [51] which 
states that if \S{p)\ = 0{p"), < o; < 1, then Pr{L) = 0{e^^^L^). The general upper 
bound for S{p) requires to take the limit a ^ 1, yielding Qij = 0(6*'/" 2), which is exactly 
like in the PGT (3.85). In section 4.4, however, it will be argued that |5'(p)| ~ 

p — > 00, so that « = I can be chosen. This yields the result Qr{1) = 0{e^H^^). Thus, 
using inverse quantum chaology in conjunction with reasonable assumptions on the fluc- 
tuations in the respective quantum energy spectra, the upper bounds on Qr{1) could be 
improved. But only in the non- arithmetic case the expectation (3.73) could be supported. 
In the arithmetic case it was merely possible to bring the bound down to e^H~3 . 

The methods of inverse quantum chaology did not work as effective as desired when 
applied to arithmetic Fuchsian groups. It will thus now be tried to support the assumption 
(3.73) on the magnitude of the remainder term in the PGT by numerical investigations 
of some arithmetic groups. Again, there are two possible approaches one could choose: a 
direct calculation of Qr{1) from length spectra, or an indirect one in the spirit of inverse 
quantum chaology. The latter way was employed by Aurich and Steiner in [14]. To explain 
this approach one has to go back to the starting point (3.86) of a modifled explicit formula 
for iI)r{L). Iwaniec succeeded in showing from (3.86) [60] that for 1 < T < e^^L""^ 

s L 

ML)^e'^+ E ^ + 0(r-iLV). (3.93) 

\Pn\<T 

The (finite) sum runs over pairs of non-trivial zeros s„ = ± ipn of Z{s). The optimal 
choice for T to resolve a given L then is T = Lei^ [60]. Inserting (3.93) into (3.83) thus 
yields 

Up{l) = Ei{l) + E E^Snl) + 0{T-He^) . (3.94) 

\Pn\<T 

The optimal choice for T results in a remainder term of 0{e^'') in (3.94). Aurich and 
Steiner, however, derived the formula 

^^p{l)^E^{l)-^E^{l/2)+ ^ Ei{sJ) + ... (3.95) 

\Pn\<T 
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by employing the Selberg trace formula with a special test function. They hence added an 
additional contribution of -~^Ei(l/2) to the r.h.s. of (3.94). In the modified explicit formula 
(3.93) this extra term docs not appear explicitly because it is hidden in the remainder. In 
the numerical evaluation of (3.95) this, however, appears to be needed. For the (arithmetic) 
regular octagon group Aurich and Steiner calculated the first 200 eigenvalues, leading to 
a cut-off at T = 14.2. Relating this to I through the optimal choice T — le*^ yields 
/ ~ 4.6. Indeed, the numerical calculations presented in [14] show a good approximation of 
A/^(/) in the interval 3 < / < 6. And furthermore, the r.h.s. of (3.95) reveals a reasonable 
approximation to the actual staircase even at larger values of /, which, however, cannot 
follow the step structure properly. But the approximation does not "leave" the staircase. 
This can be taken as an indication that the omitted remainder on the r.h.s. of (3.95) does 
not exceed the mean step height < gp(l) >~ 8-\/2^. Thus it seems that the complete 

remainder on the r.h.s. of (3.94) is of the order of magnitude of 

In order to investigate the fluctuation properties of arithmetic length spectra numer- 
ically in detail, first fiuctuations of the multiplicities will be studied. The first example 
chosen is given by Artin's billiard [68]. This is a system that is derived from a desym- 
metrization of the modular surface Tjnod\H. The latter possesses an orientation reversing 
symmetry commuting with the hyperbolic Laplacian. Dividing this symmetry out leads to 
a billiard problem in a triangle on the hyperbolic plane Ti first discussed by Artin [2]. For 
a detailed presentation of the desymmetrization procedure see [99]. The regular octagon 
group Vreg then serves as the second example. A third arithmetic group that is included in 
the numerical studies presented here is the Gutzwiller octagon group Fg- This is a group 
commensurable with Treg leading to a compact surface TgXH of genus two. The intersec- 
tion of Treg and Tg is a subgroup of index two in Treg and of index five in Tg- The mean 
multiplicity of lengths of closed geodesies on Tg\H is thus asymptotically smaller by a 
factor of I (see (3.20)) than that for the regular octagon group. For a detailed description 
of Tg see [73] . The numerical data of the lengths and multiphcities in the three examples 
have been kindly placed at disposal by the authors of [68, 4, 73]. 

As a numerical expression to calculate the mean multiplicities 

^ 9p{lp,n+k) ^ ^T^lKi^n+N) -Ki^n-N-l)] (3.96) 

^-"^ k=-N 

is taken. Asymptotically, for /„ oo, this expression approaches < gp{ln) >, compare 
(3.13). Analogously, < gpilnf > is approximated by i^'}2k=-N 9p{lp,n+kf ■ Employing 
Chebyshev's inequality, stating that (Z)^=i «n)^ < -^Z]n=i^n > 0, one merely 

ends up with the lower bound 

^s(0 :=< 9p{lf >-< 9p{l) >' > . (3.97) 

To find out about the asymptotic Z-dependence of the fluctuations c^(0 numerical cal- 
culations were performed in the three systems mentioned above. In fig.l a'^{l)Pe~^ is 
plotted for the three systems mentioned above. The respective primitive length spectra 
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are completely known up to Imax — 19.360 for Artin's billiard, Imax — 18.092 for the regular 
octagon, and Imax = 17.680 for the Gutzwiller octagon. The averaging has been performed 
over 2N = 100 lengths in all three examples. One notes that asymptotically (T|(/)/^e~' 
seems to fluctuate about a constant value so that the asymptotic behaviour of the fluctu- 

I /2 

ations appears to be |crg(OI ~ const. I — > oo. It thus turns out that the fluctuations 
of the multiplicities about their mean values do not give stronger contributions to Qr{1) 
than the mean multiplicities themselves. 

One might now be interested in seeing the behaviour of Qr{1) directly Since the 
contributions from the mean of the multiplicities and from their fluctuations are of the 
order of magnitude of flg.2 presents a numerical calculation of Q'{1) := \Qr{1) \ /e~'/^. 
As it is observed from fig. 2 that Q'{1) fiuctuates in some bounded strip one can conclude 

I /2 

that Qr{1) appears to be of the order of magnitude of ^j- in the computed range of l- 
values. It should, however, be noted that a slightly different power of / in the asymptotics 
of Qr{1) can hardly be excluded by the numerical results. Therefore, the expectation (3.73) 
is supported by numerical evidence to hold also for the arithmetic Fuchsian groups studied 
here. Numerical calculations for the modular group as well as for a non-arithmetic group 
also supporting (3.73) were performed in [29]. 

In summary, there are several reasons, both from an analytical as well as from a nu- 
merical point of view, that the remainder term Qr{1) in the PGT is asymptotically for 
/ — >• oo given by \Qr{1) \ = e2a;(/), where a;(/) is a combination of powers and logarithms of 
I, both for arithmetic and non-arithmetic Fuchsian groups. But since the mechanisms of 
fluctuations are different for these two classes of systems, the methods employed to support 
the expected behaviour of Qr{1) had to be different ones. It seems that the non-arithmetic 
case behaves like other "generic" chaotic systems (see also [87]). But for those the analytic 
tools used here are not applicable because the Gutzwiller trace formula is not an exact 
identity and hence no such detailed information about the analytic properties of the asso- 
ciated dynamical zeta functions are available like for the Selberg zeta function. Inspecting 
(1.4), however, there seems to be no fundamental difference between the general case and 
the case studied here. It is only that the powerful machinery of the Selberg trace formula 
is available for analytic investigations. The only peculiarity of the arithmetic systems that 
played a role was their exceptional structure of their respective length spectra including the 
exponential degeneracies, a property not shared by generic systems. One would therefore 
expect a behaviour like the one for non- arithmetic Fuchsian groups also for general, generic 
chaotic systems. 

4 Quantum Aspects of Arithmetical Chaos 

The item of this chapter is a discussion of the quantized versions of geodesic ffows on 
hyperbolic surfaces derived from arithmetic Fuchsian groups. The main objective thereby 
is to understand the statistical properties of arithmetic energy spectra in contrast to those 
of non-arithmetic systems. The Selberg trace formula will play a decisive role because 
it allows to trace back the exceptional spectral statistics of the arithmetical systems to 
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the peculiarities of their respective classical limits, namely to the exponentially growing 
degeneracies in their classical length spectra. 

The first section presents a discussion of pseudosymmetries in the quantum mechan- 
ical context and introduces Hecke operators. The discussion and interpretation of these 
leads to a discovery of constraints on the arithmetical quantum energy spectra that are 
taken as indications for exceptional statistical properties for the latter. Then the empirical 
observations about arithmetic quantum energy spectra arc reviewed and the role of the 
form factor for the spectral statistics is discussed. Sections 4.3 and 4.4 present a model to 
describe the level spacings distribution and the number variance, respectively. The final 
section of this chapter is devoted to an investigation of the convergence properties of the 
Selberg zeta function both for arithmetic and non-arithmetic Fuchsian groups. 

4.1 Hecke Operators 

In section 3.5 pseudosymmetries on arithmetic surfaces have been introduced as general- 
izations of symmetries. In quantum mechanics the latter ones manifest themselves as being 
represented unitarily on the wave functions. The Hecke operators are generalizations of 
the representation operators of symmetries to pseudosymmetries. For arithmetic groups 
they form an infinite algebra of self-adjoint operators commuting with the Hamiltonian 

As in section 3.5 we will now open the discussion of the realization of pseudosym- 
metries in quantum systems by briefly reviewing the case of symmetries. Let therefore 
be r a Fuchsian group of the first kind that is a normal subgroup of index in F'. 
F'/F = T, = {1, gi, . . . , gN-i} then is the symmetry group of the surface T\H. S shall be 
unitarily represented on the (finite dimensional) vector space by % e End(V^^). Since 
X may be decomposed into irreducible components, henceforth only unitary irreducible 
representation will be discussed; E* then denotes the set of these (the unitary dual of S). 
The relevant quantum mechanical Hilbert space is the space of K^^-valued square-integrable 
functions L^{T\H)®V^, which may be realized as follows. Let ip : 7i ^ transform under 
F' via the unitary representation x that is defined on F' by extending it trivially onto F, i.e. 
X is viewed as a representation of F' with F C ker(x) . Expressed in explicit terms: if 7' e F', 
then 7' e QiT for some gi e S, and, with some 7 G F, ^'{'y'z) = ip{gi'yz) — xi9i)'4'{^)- One 
is hence considering the spectral problem of the hyperbolic Laplacian on F'-automorphic 
functions transforming under unitary representations of F' that act trivially on F such that 
they yield irreducible representations of E = F'/F. The representation operators Tg. for 
g^i e E act as 

TgM^):^^{^'z)^x{9imz) . (4.1) 

On the level of the Selberg trace formula and the Selberg zeta function the desymmetriza- 
tion procedure can be carried through by decomposing the representation of F' that is in- 
duced from the trivial representation of F C F' into irreducible components, see [100, 101]. 
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The Selberg zeta function Zr{s) for the group F then factorizes as 

Zr{s)= n Zr>{s,xf^'^^ , (4.2) 

where for Re s > 1 

oo 

Zris, X) = n n det (l - xh') e-(^+")'(^')) (4.3) 

{7'}p "=0 

is the Selberg zeta function referring to T' and incorporating the representation %. (The de- 
terminant is defined on the representation space V^.) (4.2) clearly shows that an eigenvalue 
— {p^Y + I of the hyperbohc Laplacian related to the symmetry class x occurs with 
multiplicity dim V^, since this is the multiplicity of the corresponding zero s„ = ^ —iPn of 
Zr{s). The explicit treatment of an example for such a desymmetrization will be presented 
in appendix B. 

According to the algebraic setting of pseudosymmetries reviewed in section 3.5 the 
latter ones are related to (non-trivial) elements of the commensurator F of the arithmetic 
Fuchsian group F under consideration. For 5 G F the double cosets FgT are the basic 
objects the Hecke ring is constructed from, see [86, 71] as general references. Let now be 
A a semi-group, F C A C F, and form the free .^-module TZ{T, A) generated by all TgT 
for g E A, i.e. 

7e(F, A) := I ^ c{g) TgT; c{g) G X, c{g) + for finitely many (^j . (4.4) 

On 7?.(F, A) a multiplication will be introduced. To this end consider the decompositions 
(3.61) 

Tg\T = Tai U . . . U Tak , 

Tg^T = F/3iU...UFA, (4.5) 

corresponding to two elements gi,g2 G A. For every g & A. define the integer 

m{g):=i^{{t,j); Ta,(3,=Tg} . (4.6) 

The product of the cosets Fg^iF and Fg'2F is then defined as 

F^iF • F^2r := ^ m(^)F^F . (4.7) 
geA 

One can now show [86, 71] that this product does not depend on the choice of the repre- 
sentatives {cKj} and {Pj} in (4.5), and that it is associative. Furthermore, IL G A yields 
the identity element FIF for this multiplication. Extending the law of multiplication (4.7) 
linearly to all of TZ{T, A) turns this module into an associative ring with identity, which is 
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called the Hecke ring ofF with respect to A. Choosing A = F, TZ(r) :— TZ(r, F) is named 
the Hecke ring of F. In general, the Hecke ring need not be commutative, but there exists a 
sufficient criterion for its commutativity, see [86, 71]: if there exists a one-to-one mapping 
t : A A such that {gig2y = 5'2fl'i> TgiT = TgT for aU g e A, and F' = F, then the Hecke 
ring 71{T, A) is commutative. Several examples of arithmetic groups have been dealt with 
in the literature that actually do result in commutative Hecke rings, the most prominent 
ones, as always, being the modular group and its congruence subgroups, see [86, 96, 71]. 
Also the Hecke rings of arithmetic groups that are obtained from unit groups of maximal 
orders [85] or of orders of the Eichler type of level N [71] in indefinite quaternion algebras 
over (Q are known to be commutative. 

The Hecke ring 7?.(F, A) can be represented on various linear spaces; see [49, 86, 71] 
for examples not treated here. The one we are interested in is the Hilbert space L^(F\7i). 
Let therefore be -0 e L^(F\?i) and define for gf e A 

n 

T{g)i;{z):=Y.^{aiz) , (4.8) 

i=l 

where the fractional linear transformations ai, . . . , q;„ are yielded from the decomposition 
(3.61). The bounded linear operator T{g) on L'^(r\H) is called a Hecke operator, see 
[100] for proofs of its properties. Since T{g) acting on the function ip results in a linear 
combination of ip taken at points translated by operations of GL^{2, M), the Hecke operator 
commutes with the hyperbolic Laplacian. In [100] one also finds that T{g~^) = T{g)*, 
where "*" denotes the adjoint with respect to the scalar product of L^(F\7i). Thus, 
T{g) :— T{g) +T{g~^) is self-adjoint. Moreover, the T{g), g E A, form a representation of 
7e(F, A) on L'^{T\n), i.e. for ^1,^2 e A, TgiT ■ TgiT = EsgA^(^)r^r, one can show that 
T{gi)T{g2) ip{z)= EgeA m{g)T{g) -^{z). 

In case G F is a symmetry, i.e. a pseudosymmetry of order n = 1, it is an element 
of the group F' that contains F as a normal subgroup such that F'/F is isomorphic to the 
symmetry group. By (3.60) then only 71 = 1 occurs in the decomposition of F into cosets 
of r'(g) (— F). Thus the Hecke operator related to g only involves ai = g, see (4.8). One 
now observes that this Hecke operator is nothing else than the representation operator of 
the symmetry g, see (4.1), 

T{g)^{z)^^{gz)^x{9mz) ■ (4-9) 

Hecke operators therefore yield proper generalizations of symmetry operators, which is in 
accordance with the assertion that pseudosymmetries are generalizations of symmetries. 

At this point one might seek for a geometrical interpretation of Hecke operators. In 
section 3.5 pseudosymmetries of order n were demonstrated to result in n-sheeted coverings 
of the surface F\7Y. The 71, . . . ,7„ in the decomposition (3.60) of F into cosets of T'{g) 
were found to define the mappings that interchange the n sheets of the covering. Exactly 
these 7i's now occur in the definition (4.8) of the Hecke operator T{g) via ai — gji, which 
therefore in some sense averages the wave function ip{z) defined on F\7i over the n sheets 
of the covering r'{g)\H r\H. Looking at the diagram in section 3.5 one observes that 
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the (pgl'jiz) are the n points on r'{g)\H that he over (p{z) e r\H. Transforming these n 
points by g, T{g) averages over the resulting n images. 

As in section 3.5 the modular group Tmod = SL{2, ^) will serve as an example to 
illustrate the general procedure just discussed, see also [32]. An exhaustive treatment of 
this particular case may be found in e.g. [86, 96, 71] . The semi-group A needed to define 
the Hecke ring 7^(^, A) has been introduced in (3.64) as A = {g^ e M(2,X); det > 
0}. According to the convention (3.65) to look at all g e A with det g — n, n & W, 
simultaneously, the Hecke operators for T^od defined as 

Tn:=^ E T{g) , n G W , (4.10) 

dct g—n 

after normalizing them as it is commonly done in the literature, see e.g. [96]. The decom- 
position (3.65) then yields the explicit form 



for tp{z) fulfilling ipijz) = ip{z), 7 G Tmod- One can prove several properties for the Hecke 
operators (4.11), see e.g. [96]: they are self-adjoint on L'^{T\H) and form a commutative 
algebra. Their law of multiplication can be drawn from (4.6) and (4.7) as 

TnTm^ E ' n,meW. (4.12) 

d\{n,m) 

Since the T„'s, for n G IN, and the hyperbolic Laplacian all commute with one another, they 
can be simultaneously diagonalized. The square-integrable eigenfunctions of the Laplacian 
for Tmod, the so-called Maafi cusp forms, can be expanded on Ti as 

i^{z) ^NY, c{k)^K,,{2n\k\y) e^-^^ , (4.13) 

see [96]. K^{t) denotes the modified Bessel function, and p is the momentum related to 
the eigenvalue E, —AiIj{z) = EiIj{z), hj E = p"^ + j. The c(A:)'s arc the Fourier expansion 
coefficients, and is a normalization factor. The surface T,n,jdXH possesses one symmetry 
that can be realized in the coordinates employed to derive (4.13) a.s z ^ —z. Accordingly, 
the eigenfunctions of the type (4.13) can be distinguished as having positive or negative 
parity under this symmetry operation, ip±{z) = ±iIj±{—z). Thus the expansion of ip+{z) 
contains cos(27rA;a;) and the one of ip-{z) contains sin(27rA;a;) replacing e'^'^^^^ in (4.13). In 
addition, the sum then extends only over positive Ar's. The coefficients c{k) will now be 
chosen in order to yield simultaneous eigenfunctions of all the Hecke operators Tn, n E IN, 
Tn'^{z) — tnil^iz), tn G IR. Choosiug the normahzation N such that c(l) = 1, which is 
possible since it is known that c(l) 7^ 0, results in the identification c{k) = tk, see [96]. 
One thus observes the interesting property that the coefficients of the Fourier expansions 
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(4.13) of eigenfunctions of the Laplacian can be chosen to be the eigenvalues of the Hecke 
operators. This reahzation immediately yields from (4.12) the constraints 




(4.14) 



which are named Hecke relations, for the Fourier coefficients of the eigenfunctions. There 
are further properties of these coefficients and of their statistics either known, conjectured, 
or numerically investigated, see [51, 96, 53, 91] for further references. 

It is also possible to derive a trace formula for Hecke operators in close analogy to the 
Sclbcrg trace formula. The former already occurs in Selberg's original article [84] and has 
been further explored in [51]. The first numerical computations with this trace formula 
that yield the lower Fourier coefficients for some cusp forms (4.13) of definite parity have 
appeared in [28]. 

The Hcckc relations allow to express the Fourier coefficients c{k) as polynomials in the 
coefficients c{p) for primes p < k. Hence the coefficients are not independent, as it would 
be required by RMT for wave functions of classically chaotic systems [35]. There it is 
asserted that the expansion coefficients of wave functions in a generic basis are Gaussian 
random variables. The correlations induced by (4.14) are, however, known to be weak in 
the sense that [54] 

Cn(1) ^ E <k + I) c{k) = 0(7V-H-) , (4.15) 

-'^ k<N 

for all e > and every / G IV, leading to limTv^oo = 0. Although the expansion (4.13) 

uses a special basis, the Hecke relations manifest themselves in the wave functions irrespec- 
tive of any basis because the values of ip{z), being an eigenfunction of every T„, are related 
by Tn at different points (see (4.11)) for every n e ]N. Berry's random wave conjecture 
[19] now asserts that the values of wave functions for classically chaotic systems become, 
in the semiclassical limit, Gaussian random variables. The existence of Hecke operators 
inducing relations on the values of wave functions seems to indicate a violation of this 
randomness assumption. Hejhal and Rackner [54] , however, computed numerically several 
eigenfunctions of the hyperbolic Laplacian on r^ofiXH a-nd found good agreement with the 
random wave conjecture. It therefore seems that the correlations in the eigenfunctions of 
arithmetic systems induced by the presence of the infinitely many pseudosymmetries do 
not suffice to violate the conjectured random character for wave functions in generic clas- 
sically chaotic systems. This observation is in contrast to the case of symmetries, where it 
is known that in order to obtain a generic behaviour one has to desymmetrize the systems 
first. 

Although the presence of the Hecke operators seemingly does not cause an exceptional 
behaviour of the eigenfunctions, the pseudosymmetries they have been derived from could 
have an influence on the distribution of the eigenvalues. Therefore a flnal remark concerning 
the spectrum of the hyperbolic Laplacian on an arithmetic surface shall conclude this 
section. The infinitely many non-trivial pseudosymmetries of the surface r\H lead to 
infinitely many constraints for the eigenvalue spectrum ap{r) — {Eq < Ei < E2 < . . .} oi 
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—A on r\H. For every g the Fuchsian group r'{g) is a subgroup of finite index in F to 
which there is related the eigenvalue problem of —A on r'{g)\H. Denoting these spectra by 
(Tp{r'{g)), the inclusions T'{g) C F of Fuchsian groups result in inclusions o"p(F) C ap{r'{g)) 
of the respective spectra. In order to proof these inclusions one recalls that E G crp(F) is 
equivalent to the existence of a MaaB waveform ip : H fulfilling 

1. -Ai;{z) = E^{z), 

2. iIj{jz) = iIj{z) for all 7 G F, 

3. / ^ mz)\'<oo. 

The assertion that then also E e ap{y'{g)) is proved once one can show that '4>{z) obeys 
1. — 3. with T\g) replacing F for all G F: 1. is trivial, and 2. is also obviously true, since 
T'{g) C F. The requirement that T'{g) is of finite index n in F results in area(F'(5f)\7i) = 
n-area(F\7i) < 00, so that the integral in 3. gets multiplied by n when replacing F through 
T'[g) and thus remains finite. In conclusion, the discrete spectrum crp(F) for an arithmetic 
group F is a subspectrum of infinitely many spectra (Tp{T'{g)), 

<7,(F) c n ^p(r'(^)) . (4.16) 

The groups V {g) arc themselves arithmetic, since they are commensurable with the (arith- 
metic) group F. Arithmetic spectra are thus constrained by obeying infinitely many inclu- 
sions. The question now is whether these constraints are strong enough in order to yield ex- 
ceptional statistical properties of discrete energy spectra related to arithmetic groups. The 
following sections of the present chapter are devoted to trying to answering this question 
in the affirmative. The methods employed there, however, do not rely on these constraints 
but rather use the classical properties of arithmetical systems as starting points, since it 
seems to be difficult to formulate the constraints in a way that enables one to use them in 
a quantitative analysis. 

4.2 Spectral Statistics and the Form Factor 

Following the general belief on the existence of a universal classification for the statistical 
properties of discrete quantum energy spectra according to the characters of the corre- 
sponding classical systems it is expected that time-reversal invariant systems with chaotic 
classical limits possess spectra that can be described by the GOE random matrix ensem- 
ble up to a maximal scale L^ax- Berry's theory for the spectral rigidity supporting this 
assumption on scales 1 <^ L <^ Lmaxi however, presupposes that the multiplicities of 
lengths of periodic orbits asymptotically approach two for long orbits. The discussion of 
geodesic fiows on hyperbolic surfaces with arithmetic Fuchsian groups F in chapter 3, how- 
ever, revealed that for these systems the mean multiplicities of lengths grow exponentially, 
< gp{l) >~ cv^-j-i I 00. This property will certainly infiuence the spectral rigidity and 



57 



therefore the medium- and long-range correlations in the respective quantum energy spec- 
tra. Conjecturally, then also the short-range correlations, in particular the level spacings 
distributions, will be affected, since the universal behaviour of the spectral statistics for 
generic systems holds, according to empirical observations, on all scales up to Lmax- The 
question then remains what spectral statistics the arithmetical systems share? 

To get an idea what the answer might look like one can recall what is known about the 
function S{p) — - arg Z(| -|- ip) = Nfi{E{p)) for arithmetic as well as for non-arithmetic 
Fuchsian groups, and compare this with the corresponding results for classically integrable 
systems. For general cocompact Fuchsian groups the best known asymptotic upper bound 
is Nfi{E) = 0{y/E/logE), whereas the best lower bound is only Nfi{E) = 
see [51] and especially Theorem 7.10 m [50]. One hence observes a rather large gap between 
the upper and the lower bound. This can be understood once one consults the lower 
bound valid for arithmetic groups. In section 3.6 this lower bound was already employed 
and it was remarked there that although in [51] this was only proved for a certain class 
of arithmetic groups, the result extended to every arithmetic group, since it were the 
exponential degeneracies in the respective geodesic length spectra that were responsible for 
this lower bound of Nfi(E) — Q,±{Ei /logE). There hence remains for arithmetic groups 
only a much more modest gap to the upper bound. Berry's result on the spectral rigidity 
for generic systems, A^oiE) ~ ^logE, now suggests that in the non- arithmetic case 
\Nfi{E)\ ~ ■^^^ylogE, thus being close to the general lower bound. It therefore appears 
that the non- arithmetic groups obstruct the lower bound to be improved considerably, 
whereas the arithmetic groups are responsible for the upper bound. Certainly, the lower 
bound for arithmetic groups shows that Aao{E) > const. (_i,^y2 asymptotically for — >• oo, 
which clearly violates Berry's result, reflecting the fact that the presuppositions to apply 
the latter are not met. Actually, the discussion in section 4.4 yields that Aoo{E) ~ 17 1^^? 
which is only by a factor of log E larger than the lower bound. Therefore the arithmetic 
case reminds more of the saturation value of the rigidity for classically integrable systems, 
Aoo(^) ~ const. y/E. 

The integrable case also indicates that it is rather the lower bound for Nfi{E) that 
comes closer to its true magnitude than the upper one. Rigorous results for Nfi{E) are 
known for the quantization of the geodesic flow on a torus T = {]R/2nZy. The associated 
spectral problem is that of minus the euclidean Laplacian Ae — + dy acting on doubly 
periodic wave functions on Mi^, ijj{x,y) — il^{x + 27r, y) = i^{x,y -\- 27r). The spectrum is 
then given by Enm = + m^, n,m & Since N{E) is the number of points of 
inside a circle of radius oc ^/E, the estimation of N(E) is the classical circle problem, see 
[50] for a review. It is known that Nfi{E) = 0{Es-^) for some small values of 5. and 
Nfi{E) = Q±{Ei). The rigidity result is therefore in accordance with the lower bound for 
Nfi(E). 

Everything discussed so far has been concerned with A^lL; E) for L — > 00, that is with 
correlations in the spectra on very large scales. It seems that in this regime the arithmeti- 
cal systems behave more like classically integrable ones than like generic classically chaotic 
ones. It would now be interesting to learn whether or not this similarity to the integrable 
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case pertains also to smaller scales, especially for the level spacings. The first numerical cal- 
culations of quantum energies for arithmetical systems arc due to Schmit, who considered 
a special symmetry class for the spectral problem related to the regular octagon group. In 
[31, 17] he obtained a level spacings distribution that revealed a level attraction somewhat 
weaker than for a Poissonian but the -at that time expected- Wigner surmise was 

clearly ruled out. In addition it seemed that the computed P{s) would the more approach 
a Poissonian distribution the higher in energy one went. Aurich and Steiner [11] then 
calculated eigenvalues in all symmetry classes for the regular octagon group correspond- 
ing to one dimensional representations of the symmetry group for TregXH, and obtained 
the same findings as Schmit did. Hejhal was the first to compute a considerable number 
of eigenvalues for the modular group [52] and it was observed [92] that the corresponding 
level spacings behaved like the ones for the regular octagon group. Later, Schmit [82] could 
calculate more eigenvalues for the odd symmetry class on T„iod\H and further confirmed 
the results for the modular group. At that time, however, it remained unclear how the 
observed violation of the RMT hypothesis could come about and what class of systems 
would share ahke properties. Steil then computed [91, 32] 3167 eigenvalues for the even 
symmetry class and 3475 eigenvalues for the odd symmetry class comprising the complete 
spectrum up to p = 300, i.e. in energy up to = 90 000.25, revealing that P{s) can be 
rather well described for high energies, corresponding to 250 < p < 300, by a Poissonian 
distribution. In [32, 27] then the explanation was given that the arithmetic properties of 
the Fuchsian groups Freg and Pmod involved were responsible for the exceptional statisti- 
cal properties observed, and that consequently all arithmetical systems would share alike 
spectral statistics. In contrast, non-arithmetic systems were considered in [12, 14, 82, 27] 
and found to be in good agreement with a RMT behaviour of their energy spectra. Up to 
now, however, no rigorous argument or quantitative heuristics could be given that explains 
the observed phenomena. 

It is our aim to present in the following two sections a simple model that should account 
for the observed spectral properties of the arithmetical systems. The key quantity studied 
there in order to determine the level spacings distribution P{s) and the number variance 
E^(L;a;) is the spectral form factor K{t;x) [20] for the unfolded spectrum {xi}. It is 
defined as a Fourier transform of the pair correlation of the energy fiuctuations and will be 
introduced below. The pair correlation function g{t) is the two-point correlation function 
of the spectral density d{x), 



where < . . . > denotes a semiclassical averaging, as in section 2. g(t) then is the density 
function for the probability of finding an energy level in the interval (0,T), given one at 
t = 0. This can be used to construct the level spacings distribution P{s) as 



see e.g. Porter's contribution in [77]. Defining G{s) :— Jq dt g{t), one observes that P{s) — 



g{t) :=< d{x 



)d{x + -)> , 



(4.17) 




(4.18) 



59 



— ^ e and therefore 

r ds P{s) = e-«(°) - e-^(^) , T > . (4.19) 

Hence, as long as g{s) is intcgrable on any finite interval [0, T], but with dt g{t) = +oo, 
P{s) is a normalized probability density, dsP{s) = 1. 

The interesting contribution to d{x) = 1 + dfi{x) now comes from its fluctuating part 
dfi{x). Inserting this splitting into (4.17) for the pair correlation function, one is left among 
others with two terms of the form < dfi{x ± ^) >. It will now be argued that these vanish 
in the semiclassical limit. To support this idea one can go back to the regularization (2.20) 
for dfi{E). In section 2 it was demonstrated that averaging over an interval of length 
L, < dfi{E) > vanishes like L~^. Choosing now Ax — for the length of the interval 
involved in the semiclassical averaging, which meets the prerequisit 1 <C Ax <IC a: for a 
small enough power 1 > a > 0, < dfi{x ± |) > behaves like and thus vanishes in the 
semiclassical limit x oo. Hence, in this limit, 

git) ~ 1+ < dfiix - dfiix + l)> . (4.20) 

K{t; x) is now defined as the Fourier transform of the correlation function on the r.h.s. of 
(4.20), 

/+0O f f 

^ dt e^-*- < dfiix - -) dfiix + 2) > ■ (4-21) 

Since this definition involves only dfi{x), which is related via (2.20) to the periodic orbits 
of the classical system, the form factor is especially suited to be used in periodic-orbit 
theory. It only remains to reexpress the pair correlation function g{t) in terms of K{t\x) 
in order to be prepared for a periodic-orbit analysis of the level spacings distribution P{s). 
To this end one inserts 5{x) = dz e^'^^^^ into the r.h.s. of the identity 

t t \ xj y 

< dfiix - -) dfiix +^)>=Yt "^yy < ^/'(^ - 1) ^/'(^ +2^>[^iy-'t)- + 1)] , 

(4.22) 

yielding after some simple manipulations 

< dfiix - dfiix +^-) > 



Using (4.20) and (4.21) then leads to 

git) ~ 1 / dr sin(27riT) ^Kir; x) , (4.24) 

nt Jo or 

which is a fundamental expression to be used in the following section. 



X r+00 

— / dr sin(27rtr) 

Tit J -00 



2nt 

dy e'-^^y < dfiix - |) dfiix + |) > (4.23) 
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In order to derive a periodic-orbit expression for the form factor one has to rescale the 
unfolded spectrum by dfi{x) — ^-^Nfi{E) — =^^dfi{E). One can then use (2.20) and 

insert it into the definition of K{t; x) expressed in terms of E, 

K{r- E) = =^^ dX e^--'^(^)^ < dji{E - -) d,,{E + -) > . (4.25) 

Introducing the momentum variable p for convenience, (2.20) for the hyperbolic Laplacian 
on a surface T\H reads 

^ {In} k=l 

9p{ln) 

The outer sum runs over all distinct lengths 1^ of primitive closed geodesies, whereby 
their respective multiplicities have been incorporated in the amplitude factors >l„,fe. The 
possibility to include inverse hyperbolic orbits has been left open by allowing for an — 
— 1. As may be drawn from the Selberg trace formula (3.71) ordinary closed geodesies 
corresponding to hyperbolic 7 G F are hyperbolic orbits, i.e. possessing = 1. The 
consideration of inverse hyperbolic orbits is necessary in order to be able to treat Artin's 
billiard, which is obtained from TmodXH by dividing out the orientation reversing symmetry 
z 1-^ —z, thus yielding a billiard system on the hyperbolic plane Ti. Choosing the standard 
fundamental domain J-'mod = {z G H; \z\ > 1, — | < x < ^} for the modular group, 
Artin's billiard takes place on J^a = {z & ^mod', a; > 0}. Closed geodesies on T^odX^ 
that are invariant under the refiection z ^ —z then result in inverse hyperbolic orbits 
on the billiard domain, which are refiected an odd number of times on OTa, see [68] for 
further information. In addition, representations x '■ T ^ {±1} for the Fuchsian groups 
r are admitted in (4.26). For Artin's billiard it turns out that Xn = (^n- In (4.26) it 
has further been assumed that all gp{ln) geodesies 7 that share the same length /„ have 
identical X'y = Xn and = an, respectively. In Artin's billiard this requirement is indeed 
met [68]. 

Using the regularized fluctuating part of the spectral density (4.26) in (4.25), and 
performing the limit £ ^ at the end of the calculation, leads to 

see [20] and Berry's contribution in [41]. The inner sums extend over all non-zero integers, 
the prime indicating the omission of r, s = 0. 

The r.h.s. of (4.27) consists of two double sums, each running over the primitive closed 
geodesies and their repetitions. The diagonal contribution of these two double sums is 
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Naively one would expect the non-diagonal terms in (4.27) to be washed out by the semi- 
classical averaging. This to work would require large enough "random" phases e*^'^'''""*^'"^. 
However, by going to longer orbits different lengths /„ lie closer and closer due to their 
exponentially increasing density. This effect can be drawn from the exponential growth 
of Afp{l). Thus, for long orbits, the random phase argument fails, and the diagonal ap- 
proximation (4.28) is only reasonable up to some r*. For small values of r, r < r*, only 
short orbits contribute to (4.27), and these arc well separated in length. In [20, 41] Berry 
demonstrates that for t ^ oo the form factor saturates, K{t; E) 1. He further claims 
that even for r > 1 one obtains K{t;E) ^ 1. (Notice that by (4.28) already r = 1 
corresponds in the semiclassical limit to long orbits, I — Anpd{E), and thus r* <^ 1.) It 
should be remarked, however, that very recently Aurich and Sieber [8] found a violation 
of this saturation for the form factor in cases where there exists an eigenvalue Eq — for 
the hyperbolic Laplacian. They obtain an exponentially increasing K{t; E) for r ^ cxd 
instead. In such cases the contribution to (4.27) coming from the zero mode has to be 
subtracted using the trace formula, see [8]. 

For small values of r, r < r*, the diagonal approximation (4.28) indicates that 6- 
spikes determined by individual lengths In will characterize the form factor. In this regime 
K{t; E) therefore behaves non-universally. The intermediate range r* < r < 1 is then 
governed by the contribution of rather long orbits. According to Berry [20, 41] a sum rule 
of Hannay and Ozorio de Almeida [47], which exploits the uniform exploration of phase 
space by long orbits, yields K{t; E) ^ gr for t* < t < 1. Thereby it is assumed that the 
multiplicities of lengths of primitive orbits approach 'g = const, asymptotically for long 
orbits. For generic time-reversal invariant systems thus g = 2, whereas generic systems 
without time-reversal symmetry show g = 1. Obviously this sum rule cannot be applied to 
arithmetic systems because the exponentially increasing multiplicities of lengths prohibit 
< gp(l) > to approach a constant. One would, however, expect from the sum rule that 
K{r; E) grows much faster for the arithmetic systems than generically, sec also [27]. This 
observation lies at the heart of the model that will be presented in the next two sections. 

4.3 A Model for the Level Spacings Distribution 

It will now be attempted to set up a model describing the statistical properties of the 
eigenvalues of hyperbolic Laplacians on surfaces with arithmetic Fuchsian groups. The 
main tool to be employed will be the spectral form factor K{t] E). Then, as explained in 
the preceding section, the relations (4.18) and (4.24) allow to determine the level spacings 
distribution P{s) once K{t; E) is known in sufficient detail. Following Berry's reasoning as 
reviewed in section 4.2, the form factor can be substituted by its diagonal approximation 
Kd[t;E) (4.28) for r < r* ^ 1. From now on r* will be fixed at some value so that 
Kd{t;E) approximates the complete form factor sufficiently well. Then only Kd{t;E) 
will be used for r < r*. 

The semiclassical limit E ^ oo can for fixed r also be viewed as the hmit of long 
orbits, I — > oo, as can be drawn e.g. from (4.28). Thus the sum in (4.28) will be evaluated 
in the asymptotic regime Z — > oo. In particular the multiplicities gp{ln) appearing in the 
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in/2 

amplitude factors An k (4.26) will be replaced by their asymptotic values cr^-r— . Therefore 
one obtains 

An,k ^^Xn e''^'-'^'" [1 + 0{e-'^-)] , /„ ^ oo . (4.29) 

As always when dealing with periodic-orbit sums one observes also here that the {k — 1)- 
contribution to the sum over repetitions of primitive orbits is the leading one for Z„ — > oo. 
Thus, asymptotically in the semiclassical limit, one finds that 

Kj,(t; E) - i—— y 5iT 1— ) . (4.30) 

To get rid of the hardly tractable Dirac-5's one integrates (4.30), 

Introducing the asymptotic behaviour J\fp{l) ~ ^ e'/^, / — > oo, and differentiating the result 



with respect to r yields 

compare also [27]. Relation (4.32) shows that the diagonal approximation Kij{t; E) grows 
exponentially and already at a value of Tq := 2Trpl{E) ^°s(^pc?(-^)) it has reached the value 
one. To being a function of E decreases in the semiclassical limit and above some energy E — 
E it is smaller than the fixed value r*. Hence Kd{t; E) can be taken as an approximation 
for the complete form factor K{t;E) in the whole range [0, Tq] once the energy is chosen 
sufficiently large, E > E. In conclusion, the exponential behaviour (4.32) semiclassically 
describes the complete form factor in the mean up to tq. 

According to Berry's investigation of the form factor [20, 41] this approaches one for 
r oo. Since even for r = 1 mainly long orbits contribute, he argues that K{t] E) ^ 1 
for T > 1, and then takes K{t;E) = 1 to model the actual form factor in this domain, 
leading to his result for the spectral rigidity. In the sequel we will proceed analogously and 
define a model form factor 

1 , T > To 

In addition to Berry's reasoning Km{t; E) is defined to be one also on the interval [tq, 1]. 
This may be justified by the finding that Kd{t;E) already reaches the value one at 
To, whereas in the generic case Berry considers the form factor has to be interpolated 
on the interval [t*, 1] by the result obtained from the sum rule of [47], as described 
in section 4.2. Put in sloppy terms, the exponential increase (4.32) somehow exempts 
one from the need to discuss the interval [t*,1] separately. Inspecting (4.30) one notes 
that the definition (4.33) of Km{t;E) results in cutting-off the periodic-orbit sum at 
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Imax — 4:7rpd(E)ro — 2\og(^pd(E)). In the semiclassical limit the cut-off is therefore 



being removed automatically. Certainly, the actual form factor will oscillate about the 
mean value described by (4.33). This fine structure is left out in the model and thus one 
cannot really expect the resulting model distribution Pm{s) to describe the actual level 
spacings in full detail. However, the model should reproduce the facts at least qualitatively. 
Inserting (4.33) into (4.24) yields the model pair correlation function 



gj^(t) = i-^ r dT sin(27riT) e''^^"^^^)^ 

pd{E) 



1 - 



1 + 



\ogpd{E) . 



pd{E) 



Integrating this result, Gm{s) = dt gM{t), yields 



G„(.) ~ . - i {log(l^p3(E)) - 1 + j arctan(-=f^) . 

The model level spacings distribution is then obtained as Pm{s) — gM{s) e~'^'^^^\ 



(4.35) 



Pm{s) 



1 - 



pd{E) 



TT {pd{E)y + s2 
1 



exp < — s + 



TT 



arctan( 



pd{E) ■ 



(4.36) 



Consulting (4.19), and Gm{s) — >• oo for s — > cxd, one obtains that Pm{s) is a normalized 
distribution. Moreover, for £J ^ oo one finds that guis) — > 1 and Gm{s) — > s, leading to 
the observation that in the semiclassical limit Pm{s) approaches a Poissonian distribution. 
Since Pm{s) is designed to reproduce the actual level spacings distribution for E — > oo, 
also the latter is expected to converge to a Poissonian. 

The model shows a level attraction, which is for finite E weaker than that one of a pure 
Poissonian spectrum, since for s — > one finds 



Pm{s) 



gM{o) + [g'Mio)-9M{o) 
1 



7rpd{E) 



log(l^p3(B))-l + £L_J_^j 



(4.37) 



- 1- 



1 



7rpd{E) 



log{—pd{E)) 
cr 



cr 



^npd{E) 



s . 



Thus Pm(0) < 1 for finite E, and also > P'm{0) > —1. In the semiclassical limit the 
strength of the level attraction then approaches the Poissonian result P{s) ~ 1 — s. The 
model hence qualitatively reproduces the numerical findings in [31, 17, 11, 32, 27] correctly. 
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Using the modular group as an example it will now be studied how well the model 
describes the actual level spacings distribution of an arithmetical system quantitatively. 
The energy eigenvalues for the two symmetry classes occurring for the modular group have 
been kindly placed at disposal by Gunther Steil, see also [91]. The mean spectral density 
for the odd symmetry class of Artin's billiard reads [68] 

ME)-^-^^^-'-^^, (4.38) 
^ ^ 24 87r 87r ^ ' ^ ^ 

and for the even symmetry class 

1 3 1ogi? log2-41og7r 1 



I /2 

As a final input the constant cr = 1 of < gp{l) >~ is needed. Fig. 3 presents 

histograms of the level spacings in the odd symmetry class for the four momentum-intervals 
< p < 100, 100 < p < 200, 250 < p < 300, and 500 < p < 510, comprising of 341, 
1157, 1093, and 409 levels, respectively. The momenta p that have been used in the model 
(4.36), which is shown as the full curves, have been chosen in the middle of each interval, 
i.e. p = 50, p = 150, p = 275, and p = 505, respectively. One now observes from fig.3 
that the model Pm{s) approximates the actual distributions reasonably well, and that 
the quality of the approximation grows with energy. Especially, the strength P(0) of the 
level attraction is reproduced rather well, even in the low energy regime. Apparently the 
histograms as well as the model the more approach a Poissonian distribution the higher in 
energy one goes. Fig. 4 contains the same information as fig.3, but for the even symmetry 
class. Here the respective momentum-intervals contain 277, 1040, 1026, and 395 levels. 
Prom fig.4 one observes the same level of agreement of the model with the data as for the 
odd symmetry class. 



4.4 A Model for the Number Variance 

The discussion in the preceding section has revealed that the statistical properties of quan- 
tum energy spectra in arithmetical chaos can on small scales and in the semiclassical limit 
be described by those of a Poissonian spectrum. The present section now provides a con- 
tinuation of that study to medium- and long-range correlations in arithmetical quantum 
energy spectra. In chapter 2 the spectral rigidity A^i^L; x) has been introduced as a means 
to investigate spectral statistics on scales L > 1 that take several levels into account. For 
the following, however, it proves useful to study a different quantity that principally pro- 
vides the same information on spectral correlations; this is the number variance E^(L;a;), 
defined as the variance of the distribution of the numbers n{L; x) — N[x + L) — N[x) of 
levels in intervals [x, x + L], 

E2(L; x) :=< [n{L; x) - L]^ > , (4.40) 

where < ... > as usual denotes a semiclassical averaging over x. A truely Poissonian 
spectrum excels by a linear number variance, E^(L;a;) = L, whereas the GOE in RMT 



65 



possesses a E^(L;a;) that is asymptotically given by E^(L;x) ~ L for L — > 0, and by 
S^(L; x) ~ ;^[log(27rL) + 7 + 1 — ^] for L ^ 00; here 7 denotes Euler's constant. Further 
details can be found e.g. in Bohigas' contribution in [41]. 

In [21, 41] Berry presents a semi classical treatment of the number variance very much 
in the spirit of his considerations of the spectral rigidity [20] . In the semiclassical limit he 
expresses the number variance through the form factor by 

2 too rlj- 

T.\L- E)r.— — sin2(7rLT) K{t; E) . (4.41) 

TT^ JO T'^ 

Using the conclusions Berry draws for the functional form of K{t\ E) one obtains in the 
intermediate L-range 1 <^ L <^ L^ax that the number variance is given by the GOE result, 
whereas for L ^ L^ax it oscillates non- universally about a saturation value E^(£'). These 
oscillations as well as the value of T?^{E) are determined by the contributions of the short 
periodic orbits to the form factor. As can be drawn from the relation 

^3(§; E) = ^^ dr Y?{t- E) [L' - 2LV + r^j (4.42) 

of the spectral rigidity to the number variance, the saturation values of both quantities 
satisfy Ti^{E) = 2Aoo{E), see [12]. Hence the number variance can also be used to find 
out about the asymptotic energy dependence of < Nfi{Ey >= A^{E) = |E^(£^) for 
E^oo. 

The relation (4.41) now easily allows for an application of the model from the preceding 
section also to the number variance. It turns out, however, that it is possible to improve 
the form factor for the model a little in that the complete diagonal approximation is taken 
into account for r < Tq. Therefore 

A-„(r;B):={^^''([^^) ■ ^ ^ ;» , (4.43) 

with the diagonal term Kd{t; E) taken from (4.28), will be used in (4.41) to yield a model 

E|^(L;£') for the number variance in arithmetical chaos. Tj\^{L\E) then consists of two 
contributions Ti\j-^{L\E) and Tij^ 2{L', E) , the first one being derived from the integration 
in r along [0,ro], and the second one resulting from the respective integration along the 
remaining interval [tq, oo]. These integrals can be performed exactly, yielding 

Imax '■— 47rpd(£')ro, and Si (x) = J^dt^ denotes the sine integral, see e.g. [39]. The 
amplitude factor A^^k is given by (4.26). Because of the exponential vanishing of A^^k for 
k > 2 and Z„ — > oo (4.29) only the {k — 1) -contribution from the summation over k in 



1 cos(27rLro) 2L^.,^ ^ , ^ 

Si 27rLro + L . 4.44 
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the periodic orbit sums will be considered for the further analytic investigations; for the 
numerics, however, the complete double summation of (4.44) will be used. 

In the limit L ^ one can expand the sine in (4.44) and then use the asymptotic value 
(4.29) for An,k, yielding 

^mAL; E) ~ g^.^^l^^, N,{Anpd{E)T,) + 0{L') . (4.45) 

Employing the (semiclassical) asymptotics A^(/) ~ ^ e^/^, / oo, and expanding E|^ 2(-^j 
for L — > leaves one with 



npd{E) 



log{—pd{E)) - 1 
cr 



+ 0{L^) . (4.46) 



Thereby the definition tq = ^^^^^^^ log{^pd{E)) has been employed. Thus, for L — > 0, the 

model (L; E) follows in lowest order a Poissonian number variance. For finite values of 
L and large enough energies, however, it is smaller than the latter. 

Fig. 5 presents as the dotted curves a numerical evaluation of (4.44) for the odd sym- 
metry class of Artin's billiard, and compares these to the number variance obtained from 
the quantum energies computed by Steil [91]. The full curves in 5a) and c) refer to a 
sample of eigenvalues between the 1000*'* and the 2000*'* one, whereas 5b) and d) are ob- 
tained from the 2400*'* up to the 3400*'* eigenvalue. The respective momentum intervals 
are 164.92 <p< 229.70 and 250.86 <p< 296.83. The model (4.44) has been evaluated 
using momentum values from the middle of each interval, namely p = 200 and p = 270, 
respectively. The dashed and the dashcd-dotted curves provide a comparison with Pois- 
sonian and GOE number variances, respectively. One observes from fig. 5 that the model 
reproduces the small-L behaviour as well as the saturation of the actual number vari- 
ance reasonably well. It fails, however, to describe the oscillations properly; but due to 
the simplicity of the assumed form factor (4.43) no better agreement of the model with 
reality could actually be expected. A further conclusion that can be drawn from fig. 5 is 
that S^(i>; E) leaves the Poissonian form factor already for rather small values of L. The 
small-L asymptotics (4.46) of T,j^{L;E), however, reveals that the (negative) coefficient 
of vanishes for E oo. Thus the small-L behaviour of the model is Poissonian- like 
on the larger L-intervals the higher in energy one goes. This Poissonian behaviour of the 
number variance on small scales is in accordance with the findings about the level spacings 
distributions in the preceding section. Moreover, the rate of approaching a Poissonian 
distribution observed there is the same as the respective rate here, compare (4.46) with 
(4.36) and (4.37). 

A further question that e.g. arises in the context of inverse quantum chaology as em- 
ployed in sections 3.6 and 4.5 is the one for the energy dependence of the saturation value 
Aoo{E) = |S^(E). In the following therefore the limit L — > oo will be studied for the num- 
ber variance S|^(L; E) derived from the model introduced above. The second contribution 
^M,2(-^j £^) immediately yields 

El,^,iL;E) = ^ + 0{L-') , L^oo, (4.47) 
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when using the asymptotic behaviour Si (x) — ^ — + 0{-^) for x — > oo. The periodic- 
orbit term Yl\ji(L;E) contributes oscillations to the large-L asymptotics of the model 
number variance caused by the sink's. For fixed E the sum represents a superposition of 
finitely many oscillations of incommensurable wave lengths. In order to obtain the average 
value about which this superposition of oscillations fluctuates, one replaces each sin^ by 
its mean value h- Therefore 



2" 



<Ei„(L;i5)> ~ ^ E ^ 



4 f^^^^ du,{i) 



7r2 Ji, P ■ 

Since Imax — 'inpd{E)To — 2 \og{^pd{E)) ^ oo in the semiclassical hmit, one can introduce 
the asymptotics JVp{l) ~ ^ e'/^, Z — > oo, on the r.h.s. of (4.48), 

<Ei,,,(L;£)> . ^/^' f e-' 

C{li) :— ^[2^7^ — -E'i(^)] is an energy independent constant determined by the shortest 
primitive length li on r\H. The contribution (4.47) to the saturation value coming from 
E|f 2(-^; exactly cancels the second term on the r.h.s. of (4.49), 



^M,oo(^) - ^^^(log(%d(E)))+C(/l) 



2 pd{E) 
n \og{%pd{E)) 



E^oo. (4.50) 



The energy dependence of Aoo(-E') derived from the model thus is 

This result should be compared with the rigorous lower bound Nfi{E{p)) — S{p) — 
Q±(j^) [51] for arithmetic groups. Via (2.24) the latter yields 



This being a lower bound is well in accordance with the result (4.51) obtained from the 
model. Given the latter describes the actual saturation value of S^(L; E) correctly in the 
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semiclassical limit, this means that the lower bound (4.52) is off the true magnitude only 
by a factor of logE. The upper bound S{p) = O(j^), yielding A^{E) = 0{j^^^), is 
therefore much less sharp than (4.52). 

In order to test the model, the r.h.s. of the first line of (4.50) has been evaluated 
numerically and multiplied by ^ for the odd symmetry class of Artin's billiard. Fig. 6 
presents the result as the full curve. Steil has computed A^{E) from his numerically 
obtained quantum energies, see also [91]. These values arc given as dots in fig.6. One notes 
from fig.6 that the model appears to reproduce the functional form of Aoo{E) correctly, 
but that the constant C{li) = 0.075 seems to come out a little bit too large. 

Finally the above observation should be compared to very recent rigorous results of 
Sarnak and Luo [80] . They define the number variance as 

1 r"^^ o 
E|^(L; x):=- dx' [N{x' + L) - N{x') - Lf (4.53) 

X J X 

for a; — > oo, and obtain the following 

Theorem: In the range <^ L <^ Lmax-, Lmax ^/x, the estimates 

Sk(i;^) = n(j^) (4.54) 

hold for all arithmetic Fuchsian groups. 



The L-rangc these bounds refer to is the upper part of Berry's universal regime 1 -C 
L <S Lmax- The first line bounding the L-dependence is, if taken as sharp, in accordance 
with the numerical observations as well as with our model in that in any case S^(L; E) < L 
on the interval alluded to above. A direct comparison of the model (4.46) and the first line 
of (4.54) is not possible because (4.46) is only valid for L ^ 0. Since S|^(L; E) is certainly 
positive, an approximation using the r.h.s. of (4.46) up to 0{L'^) is restricted to 

TipdjE) 

^ ^ log(|^pd(^)) - 1 • ^'-^'^ 

Thus there is no overlap with the range of vahdity of the Theorem. Those two results should 
rather be viewed as being complementary. The second line of the Theorem estimating the 
-E-dependence at finite L is in accordance with the above findings about T^l^ ^{E) and, 
of course, with the rigorous estimate (4.52), although a direct comparison is also, strictly 
speaking, prohibited by the restriction L <^ L^ax iii the Theorem. 

Summarizing the results on the number variance of arithmetical systems obtained in 
this section one notices that E^(L; E) is for small values of L reasonably well approximated 
by a Poissonian behaviour. The L-range on which this agreement takes place grows in the 
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semiclassical limit. For larger values of L E^(L; E) deviates slowly from being linear as 
first described by the model (4.46) and then by the Theorem of Sarnak and Luo. At 
some Lmax oc \/E the number variance saturates and oscillates non-universally beyond 
that value. The oscillations and their mean value are determined by the short closed 
geodesies on the surface V\H. The number of distinct lengths contributing, however, is 
■N'pQ'max) ~ ^pd{E) and therefore tends to infinity for £■ — > oo. The saturation value 

< N^i{Ef >= Aoo(£;) = iS^(E) ~ f £; ^ oo, grows slightly less (by a factor 
of j^^) than the corresponding value for classically integrable systems. It is, however, 
certainly well beyond the one for generic classically chaotic systems, Aoo(£') ^ -^logE. 

Finally a remark on the class of systems that are being described by the findings 
of the present chapter will be added. The reason why the model works is provided by 
the exponentially increasing multiplicities of lengths of closed geodesies for arithmetic 
Fuchsian groups, since the exponential behaviour (4.32) of the form factor is caused by 
the compensation of the exponential damping present in the amplitude factor A„ ^ (4.26) 
through gp{ln)- In order this to work it has been assumed that all gp{ln) closed geodesies of 
the same length /„ shared alike factors x{l)- Otherwise the sum over closed geodesies could 
not have been rewritten as a sum over distinct primitive lengths, see (4.26). Assuming the 
simplest case of .^2-valued representations x : F {±1}, one can group the geodesies in 
classes of alike signs and define gp{ln) — gn'^9n-i where g^ denotes the number of geodesies 
of length In with x{l) — il; respectively. The diagonal term (4.28) then reads 

1 °° hi 

{^TipdiE))^ ^^^^ 47rpd{E) 

(4.56) 

The leading (k = l)-contribution thus contains the difference [gn~gn]^ of the multiplicities 
referring to ^(7) = +1 and xil) = Only if this difference grows like \g^ ^ 9n \ ^ 

const. ^tlH^ 00, the model of sections 4.3 and 4.4 is applicable. Once g^ and g^^ are 

of the same order of magnitude the leading contribution comes from k = 2. Since then 
[9n + {~^y9nV — 9p{^nf ~ Cp this term is of a similar form as the analogous one 
for A; = 1 in the non-arithmetic case. Thus the statistical properties are expected to be 
generic, i.e. the level spacings should be close to the GOE behaviour and the medium- and 
long-range correlations should be described by Berry's theory. 

As an example for a non-trivial representation % take an arithmetic Fuchsian group Fi 
leading to a symmetric surface Ti\H. Then Fi is a normal subgroup of index in another, 
also arithmetic, group F2. According to (3.20) thus < gf\l) >~ < g'^jPQ) >, I ^ 00. 
The symmetry group F2/F1 is represented via x '■ ^2 ^ End{V^) with kerx ^ Fi. Since 
therefore ^(7) = +1 for 7 G Fi, one concludes that g^ > g^\ln)- Again only .^2^valucd 
symmetry classes shall be considered for simplicity. Then x{l) = ~1 is only possible for 
7 e F2, hence < g- ><< gf\ln) >~ ^ < ^^H^n) >, yielding < ^+ > - < > > 
(1 — ;^) < g^^\ln) >~ const. for 00. Therefore the model is still apphcable to 
this case. 

The more general situation of arbitrary representations of a symmetry group can be 
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treated analogously. In view of the possible cancellations of multiplicities in the {k — 1)- 
term of (4.56) the case of a ^2^valucd representation, however, is the worst possible. Hence 
the Laplacian on a symmetric arithmetic surface is always expected to share statistical 
properties as discussed in the preceding sections. Artin's billiard may serve as an example 
that has already been studied above. Although the symmetry on the modular surface is 
orientation reversing and thus r2 is not a subgroup of SL{2,IR) but rather of the full 
group of isometrics of the above reasoning extends also to Artin's billiard, since inverse 
hyperbolic orbits occurring due to orientation reversing symmetries can as well be dealt 
with in the present framework, see [99, 100, 68]. 

There may of course exist more general representations of arithmetic Fuchsian groups 
r than the ones being derived from symmetries. The latter excel by their triviality on 
a subgroup of finite index, which leads to the observation just made that a large enough 
fraction of closed geodesies is equipped with positive x's. Examples for the former may be 
provided by the presence of Aharonov-Bohm type magnetic fluxes on arithmetic surfaces 
V\H. Depending on the strength of such a flux the spectral properties of the respective 
Laplacian are expected to deviate from the findings of the present chapter. Once the 
phases x{l) that arise when a wave function is carried along the geodesies related 
to the 7 G r and enclosing the Aharonov-Bohm fiux line "mix" sufficiently among those 
geodesies being degenerate in length, one could even retain generic spectral statistics like 
those for non- arithmetic groups. 

4.5 Convergence Properties of the Selberg Zeta Function 

The final topic of the present chapter now again deals with general Fuchsian groups of the 
first kind. However, it is observed in the course of the following discussion that arithmetic 
groups play a special role. As in section 3.6 methods of inverse quantum chaology are 
applied and once again it turns out that the strong fiuctuations present for arithmetic 
quantum energy spectra violate the prerequisits to apply the formalism developed below 
to the arithmetic case. Referring to the heuristic reasoning presented first in [6] one can, 
however, understand the reason for the obstruction occurring for arithmetic groups. 

The item of this section lies at the foundation of one of the major objectives of quantum 
chaology, namely the derivation of certain quantization rules that allow to determine the 
quantum energies of a classically chaotic system in a semiclassical approximation. Recently 
such quantization rules involving dynamical zeta functions have been introduced and suc- 
cessfully apphed to a variety of different chaotic systems, see e.g. [89, 68, 95, 26, 63, 23, 13]. 

All these methods make strong use of the fact that the semiclassical quantum energies 
are directly related to the zeros of the dynamical zeta function on the critical line. The 
problem one immediately faces when trying to compute the non-trivial zeros explicitly 
is that the Euler product (2.9) defining the dynamical zeta function in general does not 
converge on the critical fine. This phenomenon is also referred to as the existence of 
an entropy harrier, since it is the topological entropy r that determines the half-plane 
of convergence Res > r for the Euler product. Concerning geodesic flows on hyperbolic 
surfaces and the Selberg zeta function it was mentioned earlier that the topological entropy 
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universally is r = 1, and that the critical line is located at Res — ^. The entropy barrier 
to be overcome hence has a width of |. 

In order to find a quantization rule one therefore has to develop a consistent procedure 
to calculate the non-trivial zeros of the zeta function other than searching for the zeros 
of the Euler product. It was McKean [69] who apparently first noticed the possibihty to 
rewrite the Euler product (3.72) of the Selberg zeta function as a Dirichlet series, but made 
no use of this. Berry and Keating [22] then were the first to introduce this Dirichlet series 
in order to obtain a quantization rule from it, leading to their Riemann-Siegel look alike 
formula. Sieber and Steiner [89], Matthies and Steiner [68] and Aurich and Steiner [13] 
then investigated the convergence properties of the Dirichlet series in several examples and 
used it to calculate non-trivial zeros. After that in [6] a statistical model was developed 
that predicts the domain of conditional convergence for the Dirichlet series of the Ruelle- 
type zeta function (3.74), which as well can be used to set up a quantization rule. This 
model can also be extended to the Selberg zeta function, see e.g. [5]. 

The above remarks stress the importance of investigating the convergence properties of 
the Dirichlet series representing the Selberg zeta function. Very much alike the discussion 
in section 3.6 methods from analytic mimber theory will be used below in conjunction with 
inverse quantum chaology to derive statements about the abscissa of convergence for the 
Dirichlet series. Again, as in section 3.6, for simplicity only cocompact Fuchsian groups F 
will be considered. But as also has been stated earher, the results do not depend on this 
restriction, since it is only the contributions of the hyperbolic 7 e F that are relevant. 

Before entering the detailed studies just announced the model introduced in [6] will be 
reviewed. In order to convert the Euler product (3.72) into a Dirichlet series one has to 
transform the product over n e IVq with the help of Euler's identity 
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n (1 - yx") = 1 + 5: ^ ' 

n=0 m=l llfc=l(.-'- ^ ) 

first; |a;| < 1 and y G (D suffice for the product and the sum to converge. One then obtains 
that for Res > 1 

Z(s) = 5:^,6-^" , (4.58) 

p 

where the sum extends over all pseudo-orbits p, see (3.77), and the Lp's denote the pseudo- 
lengths (3.78). The coefficients Ap of the generalized Dirichlet series (4.58) are given by 

Ap = n«» ' «i ■= nrA(i-e-''^o ' ^^'^^^ 

for p — {'fT^jp © ... © {'yr'"'}p- Recall that a pseudo-orbit p is a formal combination of 
finitely many (not necessarily primitive) closed geodesies, the integers mi , . . . , denoting 

the numbers of respective repetitions of the primitive closed geodesies corresponding to 
{71 }p, . . . , {7r}p. The asymptotic relation A/'(/) ~ A/^(Oi ^ ~^ C)0, which is an analogue of 
(3.9), expresses the fact that the number of all closed geodesies of lengths not exceeding 
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I is in the asymptotic regime I — > oo already given by the respective number of primitive 
geodesies. In the hmit L — oo therefore the number jV^^^ (L) of pseudo-orbits with Lp < 
L is dominated by those p that are completely composed of primitive closed geodesies 
traversed only once, i.e. mi — ... — rrir — 1. Their respective coefficients Ap look like 

^P = flr^~(-1)'^ (4-60) 

where \p\ := r denotes the number of primitive geodesies the pseudo-orbit p consists of. 

Since already the lengths of primitive closed geodesies are at least twice degenerate 
because of the time-reversal symmetry, the pseudo-lengths L„ will in general also possess 
multiplicities exceeding one. Assuming that these multiplicities g^^\Ln) are exclusively 
caused by degenerate primitive lengths and not by different combinations of primitive 
lengths yielding the same pseudo-length, degenerate p's have identical coefficients Ap. Thus 

Z{s)= Y.9^''\Ln)A^e-^^- (4.61) 

for Re s > 1. Numerical calculations now hint at an exponential increase of the multiplic- 
ities according to [6, 5] 

< ^(^)(L) >~ rfe"^ , L^oo; (4.62) 
d and a are constants to be determined numerically. Thus 

^(^)=E^^n^^e-(-)-", (4.63) 

{Ln} 

which can be considered as a generalized Dirichlet series in the variable s — a. The theory 
of Dirichlet series [48] as e.g. also briefly reviewed in appendix A now allows to determine 
the convergence properties of the Dirichlet series for Z{s) from (4.63) once the pseudo- 
orbits have been arranged in ascending order of their respective pseudo-lengths, = Lq < 
Li < L2 < . . .. According to this theory (4.58) converges for Res > and diverges for 
Res < Gc\ furthermore it converges absolutely for Res > aa, (Ja > Oc- The abscissa of 
convergence Oc and of absolute convergence a a are determined by the formulae (A. 16). The 
asymptotic growth (3.89) of the number J\f^^\L) of pseudo-orbits with pseudo-lengths not 
exceeding L then fixes the abscissa of absohite convergence to be iTa = 1. In [6] a statistical 
model was established that yielded (Jc = ^-j^. The model assumes that after applying the 
approximation (4.60) the coefficients An in (4.63) represent random signs. This conjecture 
is based on the observation that by (4.60) Ap ~ ±1, the sign depending on whether 
the number of primitive closed geodesies comprising the pseudo-orbit p is even or odd, 
respectively. Arranging the pseudo-lengths in ascending order and taking into account a 
supposed irregularity in the distribution of primitive lengths should make the numbers 
and |pn+i| modulo two independent of one another. 
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Taking this randomness hypothesis for granted one can obtain the result for ac also 
slightly differently, although in the same spirit as in [6]. For it is known [61] that a series 

CXD 

Y^i-ir^c, (4.64) 

k=l 

of positive coefficients and random signs (—1)'^'= either converges or diverges, depending 
on whether the series J2k'=i c\ converges or diverges, respectively. Recalling (3.89) and 
(4.62), one observes that the number of distinct pseudo- lengths up to a value of L grows 
asymptotically proportional to e^^~"^^^. Thus the criterion for convergence stated after 
(4.64) requires Res > for the Dirichlet series (4.63) to converge, reproducing hence 
the outcome of the model in [6] . In conclusion, one learns from this model that the Dirichlet 
series representing the Selberg zeta function is not expected to converge on the critical line 
Res — \. The distance of Uc to the critical fine is determined by the growth of the 
multiplicities of pseudo-lengths. 

For arithmetic systems oc is expected to be large since already the multiplicities of 
primitive lengths grow exponentially. Indeed, by numerical calculations of pseudo-length 
spectra up to some cut-off value Lmax it was observed in [6] that a = 0.4658 for the regular 
octagon group, and a — 0.279 for Artin's billiard. However, non- arithmetic and completely 
desymmetrized systems should possess a considerably smaller value for a. In [5] an example 
for such a system was studied numerically and a = 0.0572 was found. It even may be that 
in these cases the multiplicities will not really show an exponential behaviour, but will 
rather follow a power law, leading to an effective vanishing of a. Then the Dirichlet series 
would converge for Res > ^ and diverge for i?e s < |, but it would not be known whether 
it converged on the critical line. In any case it would be possible to evaluate Z{s) close to 
the critical line and to obtain the non-trivial zeros as minima. One could even try to go 
onto the critical line and hope that a divergence would not show up when using the finitely 
many available pseudo-lengths in (4.58). An example of such a procedure is presented in 
[5]. 

A remark on the (arithmetic) case of Artin's billiard seems to be in place now. It is 
expected that will in this case be well above |, keeping in mind the rather large value 
of a = 0.279. However, a numerical evaluation of the formula (A. 16) for cTc yields a value 
below I [68, 6], at least in the finite range of available pseudo-lengths. An explanation for 
this observation is that the randomness hypothesis the statistical model is built upon is 
apparently violated [6] in the computed range of the pseudo-length spectrum. Thus Artin's 
billiard can apparently not be understood by the above considerations. The question 
whether or not this phenomenon pertains to higher values of L remains open. 

In the following the question for the location of the abscissa of convergence ac of the 
Dirichlet series (4.58) will be approached from a different side, employing methods from 
analytic number theory and inverse quantum chaology. As guiding references [97, 51] may 
be consulted. The idea to be pursued below is similar to the one how to obtain the PNT 
from the analytic properties of (see (A. 12)), or how to obtain the PGT from (see 
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(3.84)). To this end define the function 



p,Lp<L 

where the above notation should indicate that the pseudo-lengths have to be counted with 
their respective multiplicities. The abscissa of convergence is then according to (A. 16) 
given by 

CTc = limsup — log \ipz{Lp)\ . (4.66) 

Lp— »oo J^p 

If it were possible to derive an 0-estimate for ipziL), this would yield an upper bound for 
(Tc, i.e. ipz{L) = 0(L"e^^), a & M, b > 0, results in the bound < b. Accordingly, an 
f2-result for ipz{L) would give a lower bound on ac- Thus it will be attempted to estimate 
ipz{L) for L — > oo. The principle tools to be employed have already been used in section 
3.6 and in appendix A. 

Using the Dirichlet series (4.58) and the integral (A. 11) one easily obtains 

1 rc+ioo fj Q 

^P,{L)^— -e^^Z(.), ol, (4.67) 

Am Jc-ioo s 

in analogy to (3.84) and (A. 12). However, it does not proof particularly useful to extend 
the contour of integration in (4.67) from c — ioo to c+ioo. In the following the integral will 
therefore be restricted to the finite interval [c — iT, c + iT], T > 0, and the remainder that 
has been left out in comparison to (4.67) will be estimated. This is achieved by the following 

Lemma: With the notations introduced above, and c> 1, L > 0, L ^ Lp, 



Proof: The proof is a standard calculation in analytic number theory, see e.g. [97], pp.60. 
Due to the importance of the result for the further considerations the main ideas shall, 
however, be reproduced here. 

Inserting the Dirichlet series (4.58) into the integral on the r.h.s. of (4.68) one is left 
with an integrand of ^e^^^^^''^ Depending on the sign oi L — Lp one has to choose different 
contours to render the following integrals finite: 

1. L>Lp: 

-i. / +/ +/ pe^(^-^^) = l, (4.69) 

ZTll yJ-oo-iT Jc-iT Jc+iT J S 

because the pole of the integrand at s = is enclosed by the contour. 
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2. L<Lp: 



2ni [Joo-iT Jc-iT Jc+iT J S 



(4.70) 



A typical integral to be estimated now can be treated by an integration by parts (L > Lp) , 



''-iT ds 



/c—ii 
-oo—i 



oc-iT S 



siL - Lr, 



+ 



1 



-oo-iT Lp 



^-^T ds 



da^ ^ . (4.71) 



{c-iT){L- Lp) L-Lp J-oo {(J-iT) 

To obtain an upper bound for the absolute value of the above expression one extracts the 
maximal value of the exponential under the integral on the r.h.s. and obtains 



'^-i'T ds 



AL-Lp) 



oo-iT S 



AL-Lp) 



< 



= o 



iT\(L-Lp) L-L 

^c{L-Lp) \ 



^c{L-Lp) 

+ — / d(T 



T{L-Lp)^ 



(4.72) 



The remaining three integrals from (4.69) and (4.70) are of the same type and obey the 
same bounds. Thus 



1. L>Lp: 



2. L<Lp: 



1 [<^+iT ds L) . ^ / e 
/ — e^'^^'^"' = 1 + 

27rZ Jc-iT s 



1 



ds 



AL-Lp) _ Q 



c(L-Lp) - 

T{L-Lp)^ 

^c{L-Lp) \ 



2m Jc-iT s V-^(-^ ~ Lp) J 

Using these bounds one finds with the help of (4.65) that 

1 pc+iT d<i I P^^ p~<^^P 



(4.73) 



(4.74) 



(4.75) 



One is therefore left with the task of bounding the sum on the r.h.s. of (4.75). Somewhat 
tedious but straightforward calculations that may e.g. be found in [97], pp.60, and will not 
be reproduced here yield 



\L-Lr. 



o 



c-1 



+ 0(Le(^-^)^) , 



(4.76) 



finally proving the lemma. 
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The next task in order to derive a bound for ipz{L) is to estimate the integral that is 
left on the r.h.s. of (4.68). Because of the factor of e*^ under the integral it proofs useful to 
move the contour from Re s = c > 1 to the left in the complex s-plane as far as possible. It 
turns out that this can be achieved up to directly before the critical line. What is missing 
yet is an estimate of Z{s) on the contour. It is at this point where an inverse quantum 
chaology argument enters the game. Namely, Hejhal can prove an estimate for Z{s) in the 
half-plane Res > ^ depending on an upper bound for S{p) — Nfi{E(p)) — ;^argZ(| + ip). 
Define 

A(p) := p^'i^logpYiloglogp)^ , (4.77) 

where the exponents jj,, i^, A are chosen such that A(p) tends to infinity for p ^ oo. Using 
the notation s = a + ip, a e IR, p > 0, Theorem 10.10 in [51] then states that |<S'(p)| = 
0(A(p)) imphes 

log Z{s) = O (a(p)2'"^[°'^-'^1 log A(p)) , (4.78) 

ior (7 — Res >\ + i^gAip) ^^'^ P large enough, p > po{A). It will henceforth be assumed 
that po(A) is chosen that large that A(p) is monotonically increasing for p > po(A). Prom 
now on a shall be restricted to the domain | + i^gA^p) < i + logiogA(p) < cr < 1, on 
which 2max[0,l - a] = 1 - i^^i^fAp- Thus logZ(s) = C>(A(p) e"i^^^M+'°siogA(p)^ ^ 

_ logA(p) , 1 

0(A(p) e i°gi°gA(p) ^_ Since therefore | logZ(s)| < i] A{p) logiogA(p) fQj. gome constant 77 > 
and p > po(A), which again must be chosen large enough, one obtains that | logZ(s)| < 
eA{p) for all £ > and p large enough. Thus in the domain alluded to above 

|Z(s)| < e"^^^') , V£>0. (4.79) 

The integration contour in (4.68) can by Cauchy's theorem now be moved to Re s = ^ + S, 
5 :— log log A ; without loosing control on the magnitude of Z{s), 



e^^Z(s) = — { + + ^ — e^^Z(s). (4.80) 

2m Jc-iT S 2m [Jc-iT Jl+S-iT J^+S+iT) s 

The first and the third integral on the r.h.s. behave alike, 

■h+S±iT ds 



— e^^Z(s) < / da -——\Z{a±iT) 



o(^Le-HT)ecL'^ . (4.81) 



The second integral, however, can be bounded according to 



J^+S-iT S J-T 



,1 + 5 + 

O (^e(5+^)^ e^^(^) £ dt ^ j (4.82) 
C>(e(5+5)^e^'^(^)logr) . 
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Combining the estimates (4.81) and (4.82) with the lemma one obtains the following 

Proposition: The function ipz{L) can be estimated for L ^ oo using the notations 
introduced above as 

ML) = O (^e^^(^)e^^) + O (e(H^)^e^^(^) log t) + O (^e^^) + O (^e^) . (4.83) 



Recall that one seeks for an upper bound of the type '4'z{L) = 0(L"e''^), 1 > 6 > 0, in 
order to bound the abscissa of convergence by a"c < 6 < 1. The third and the fourth term 
on the r.h.s. of (4.83) therefore require to take T = e'^^ for some appropriate o? > c — 1 > 0. 
Then, however, A(T) = e^^'^^{dLy (log dV)^, and the first two terms prohibit to obtain the 
desired form of the estimate unless // = 0. Once, however, A(p) = (log p)'' (log log p)^ for 
i/ < 1 or 1/ = 1, A < 0, and thus \S{p)\ — O(logp), one observes with the choice d — c>l 

iJz{L) = 0{e''^) + 0{L e(i+'5+^^)^) + 0{L e-^^-^''^) (4.84) 

for L — oo and for all e > 0. Since 6 = — , '^^,rr.^ vanishes for L — cxd and e can be made 

log log A(T) 

as small as required, one draws from (4.66) determining the abscissa of convergence the 
bound 

ac<^ + e' for all e' > . (4.85) 

Therefore, the Dirichlet series (4.58) for the Selberg zeta function converges (conditionally) 
for all s with Res > ^, since one can then always choose e' as small as desired. 

This being a conditional result, the question for the vahdity of the input \S{p)\ — 
\Nfi{E{p))\ — O(logp) immediately arises. For arithmetic groups the lower bound already 
employed in sections 3.6 and 4.2, S{p) — ^±(i^), forces to chose > | in (4.77), therefore 
ruling out an application of (4.83) to obtain an upper bound for ac- This negative observa- 
tion comes in accordance with the result obtained from the statistical model, cTc = > |, 
where a, describing the growth of the multiplicities of pseudo-lengths, is rather large for 
arithmetic groups, and thus ac violates the lower bound (4.85). The lower bound for gen- 
eral (cocompact) Fuchsian groups, S{p) — ^±{\/^^^)-, still allows for expecting /j, — 
and u <1. Now suppose that S{p) = 0((logp)^), i.e. Nfi{E) = C>((log£;)^). In order to 
obtain the saturation value Aoo(-E') of the spectral rigidity in the semiclassical limit (see 
(2.24)) one has to evaluate 

— j^J dE' [Nf,{E')f = OiilogEr) + O i^{logEr-'—j (4.86) 

in the limit L — > cxd and E ^ oo. The interval of length 2= to be integrated over has to 
be kept small compared to E. Choosing L — E"-, < a < 1, then yields the semiclassical 
asymptotics 

^^{E) 0{{\ogEf^) . (4.87) 
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Berry's semiclassical theory for the spectral rigidity yielding for generic classically chaotic 
systems A^{E) ~ ^ \ogE now implies z/ = | < 1. Once one believes in the applicability 
of this heuristic theory to non-arithmetic Fuchsian groups one has to draw the conclusion 
that the Dirichlet series (4.58) for the respective Selberg zeta functions converge condition- 
ally for i?es > |. The above reasoning can be supported by the numerical evaluation of 
the formula (A. 16) in an example of a non- arithmetic group in [5], which yielded a result in 
accordance with the bound (4.85). From the part of the pseudo-length spectrum calculated 
in [5] one can, however, not draw a clear-cut conclusion on the precise value of Uc because 
the plotted curve still oscillates rather strongly. 

A further confirmation of the above results may be provided by considering the Riemann 
zeta function C(*)- Assuming the Riemann hypothesis its non-trivial zeros are given by 
Sn — \-^i^n-i In > 0. Supposing that either 7„ or 7^ correspond to quantum energies of a yet 
unknown physical system, the spectral statistics show a behaviour as if the classical limit of 
this system were chaotic without time-reversal invariance, i.e. the level spacings can be well 
described by the GUE random matrix ensemble, and the spectral rigidity and the number- 
variance saturate for L — > 00, see e.g. [74, 21]. The function S{p) — ^ argC(|-|-ip) is known 

to obey the same lower bound S(p) — ^±{\f^^^) [72] as the analogous quantity for a 
generic Selberg zeta function. Its upper bound, however, is given by S(p) — 0(j^|^), see 
e.g. [50]. Berry's theory of the rigidity stating that A^{E) ~ log E would predict that 
l^d?)! is asymptotically given by y/\ogp, possibly times some power of loglogp, therefore 
clearly being within the rigorous bounds. This is also in accordance with the belief that 
the f2-estimate is "likely to be best possible" [72] and possibly sharp. The upper bound 
for S{p) leads to an analogue of (4.78), namely logC(s) = 0((logp)^°^^[°'^~°'l loglogp) for 
a — Res > ^ + loglogp - The estimate corresponding to (4.79), ({s) — 0{p^), \/e > 0, is 
then equivalent to the Lindelof hypothesis, see [97] for details. 

In conclusion it appears that as long as the geodesic flow on a hyperbolic surface T\H 
is generic in the sense that the spectral statistics can be described by RMT on small scales 
and follow Berry's prediction for the rigidity, especially if the prescribed saturation occurs 
with the predicted energy dependence, the Dirichlet series for the Selberg zeta function 
converges for Re s > |. The arithmetic systems once again drop out of this general scheme 
by reasons that seem to be understood: sticking to the statistical model they are provided 
by the exponential degeneracies of the (pseudo-) length spectra, whereas consulting inverse 
quantum chaology it is the exceptional spectral statistics, showing much less correlations, 
that cause the trouble. Hence the arithmetic case once again exemphfies the duality of 
classical and quantum properties and demonstrates the mechanisms of their interplay. 

It is tempting now to express the expectation that generic systems with a chaotic 
classical limit might show an analogous behaviour regarding the convergence properties 
of their dynamical zeta functions. As long as their phase spaces are compact and Pesin's 
theorem [75] on the equality of metric and topological entropy holds, the Dirichlet series 
for the dynamical zeta functions should converge conditionally for Res > f = f , i-e. up to 
immediately before the critical line. The expectation is based on the fact that although the 
above inverse quantum chaology argument required detailed knowledge about the analytic 
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properties of the Selberg zeta functions, the mechanism seems to be universal as it dwells 
on the magnitude of Nfi{E). One only needs a "rigid" spectrum with a saturating spectral 
rigidity that can be described by Berry's general scheme. For generic chaotic systems all 
evidence is for this to be satisfied, and thus, although the knowledge about the associated 
dynamical zeta functions is much poorer than for the Selberg zeta function, at least there 
seems to be no obvious obstacle around against the expectation to hold. One is, however, 
far from proving this since the technical problems are enormous in the general case. 

5 Summary 

This investigation contained a discussion of the semiclassical quantization for a class of 
strongly chaotic systems. The relevant aspects of classical and quantum mechanics for the 
unconstrained motions of single particles on hyperbolic surfaces with arithmetic fundamen- 
tal groups were studied. The main body of this text consisted of two major parts: chapter 
3 discussed classical mechanics, i.e. the geometry of the surfaces the systems are defined 
on; chapter 4 then was devoted to an investigation of the quantum mechanical energy 
spectra of the arithmetical systems. It was worked out that the arithmetic nature of the 
fundamental groups involved had consequences for the geometry of the respective surfaces. 
In particular the length spectra of closed geodesies reveal high degrees of degeneracies. 

In the context of the semiclassical quantization of classically chaotic systems the philoso- 
phy of studying "generic" systems includes the requirement of only dealing with completely 
desymmetrized systems. Symmetries may lead to unwanted effects that superimpose the 
structures one tries to explore, and in many cases they can rather easily be removed. Once 
the systems have been desymmetrized one can compare them irrespective of their differ- 
ences in detail. Those quantities that appear to share common properties can then be used 
to characterize the class of chaotic dynamical systems. One manifestation of a symmetry 
in a classical system is the presence of degeneracies in the spectrum of periodic orbits with 
respect to their actions. After removing symmetries two classes of generic systems remain. 
The first one comprises of time-reversal invariant systems, whereas the second one consists 
of systems without time-reversal invariance. The philosophy referred to above continues in 
assuming that then periodic orbits generically can at most be twofold degenerate in action 
due to a time-reversal invariance. Further multiplicities would be considered as accidental. 
The discussion of the arithmetical systems, however, revealed that there exist perfectly 
chaotic Hamiltonian dynamical systems with multiplicities of lengths of periodic orbits 
that even grow exponentially with length. These multiplicities are not really accidental 
since they can be traced back to the structure of the set out of which the geodesic lengths 
are allowed to be taken. The arithmetic structure inherent in this set forces the lengths of 
closed geodesies not to cluster too strongly for Z — > oo. Since the total number of closed 
geodesies with lengths up to I has to grow according to the universal prime geodesic theo- 
rem, N{1) ~ y , / — > oo, the low number of distinct lengths up to I, Jv{l) ~ ^ e'/^, / — > oo, 

I /2 

has to be compensated by exponentially increasing multiplicities, < g{l) >~ cr^, I — > oo. 
In chapter 3 details on the mechanism resulting in the exponential law for the mean mul- 
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tiplicities occurring for arithmetic fundamental groups have been worked out. Since there 
exist infinitely many arithmetic Fuchsian groups these exceptional systems form a whole 
class of strongly chaotic dynamical systems that cannot be neglected. In addition, there is 
a further (discrete) dynamical system known with similar properties. This is the so-called 
cat map, whose classical and quantum properties have been discussed in detail by Keating 
[62]. 

It has long been known that the arithmetical systems excel by a further property, 
namely by the existence of infinitely many pseudosymmetries. Although their defini- 
tion includes geometric symmetries, non-trivial pseudosymmetries cannot be removed in 
a kind of desymmetrization procedure. It might appear that the impossibility to "de- 
pseudosymmetrize" the arithmetical systems suffices to consider them as generic, but the 
discussion of multiplicities in their length spectra showed that the algebraic and geometric 
structures induced by pseudosymmetries are important enough to result in considerable 
effects. On the quantum mechanical side of the problem, which was discussed in chapter 4, 
these structures affect quantities that are commonly be considered as characteristic for a 
distinction of classically chaotic and integrable systems. In particular statistical properties 
of quantum mechanical energy spectra were discussed. Regarding the latter, symmetries 
manifest themselves as independent superpositions of spectra referring to individual sym- 
metry classes. Thus the total quantum energy spectrum of a system possessing a discrete 
and finite symmetry group contains finitely many subspectra, which can each be viewed as 
generic, since they are spectra of desymmetrized systems. It was observed that regarding 
non-trivial pseudosymmetries the situation is somehow reversed. The eigenvalue spec- 
trum of an arithmetical system is a subspcctrum of infinitely many other, also arithmetic, 
spectra. If it were an independent superposition of infinitely many spectra, one could im- 
mediately identify the result as showing Poissonian fluctuations, see e.g. [70]. However, the 
fact that it is only a subspectrum in such an inflnite superposition complicates the use of 
this point of view to draw conclusions on the spectral statistics for arithmetical systems. 

This is one reason for chapter 4 to proceed differently in its investigation of arithmetic 
energy spectra. It appeared to be more convenient to employ the exponential growth of 
multiphcities in the geodesic length spectra. The spectral form factor turned out to be 
a useful means for a periodic-orbit investigation of spectral statistics. The two quantities 
that were picked out to be studied were the level spacings distribution and the number 
variance. The former yields information on short-range correlations, whereas the latter 
takes medium- and long-range correlations into account. The exponential increase of mul- 
tiplicities of lengths allowed for the development of a simplified model for the form factor. 
This model was essentially only based on the obtained exponential increase for small r and 
the saturation for r — > cxo of the form factor. Applied to the level spacings distribution 
and to the number variance the model was found to describe the numerically observed 
phenomena qualitatively correctly. Quantum energy spectra of arithmetical systems arc 
reminiscent of those for classically integrable systems. Their fiuctuations are much stronger 
than those for generic classically chaotic systems. They show a level attraction that grows 
with increasing energy and the level spacings approach a Poissonian behaviour for E ^ oo. 
This finding is in contrast to the integrable case that seems to yield stationary distributions 
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already at finite energies. On larger scales the correlations in arithmetical spectra appear 
to be slightly stronger than those observed for integrable systems. This is reflected in the 
energy dependence of the saturation value of the spectral rigidity. The latter was found 
by Berry [20] to be Aoo{E) ~ const. \fE, E oo, whereas the model for the number 
variance of arithmetical systems yielded A^oiE) ~ ^ i^^; E ^ oo. However, the spec- 
tral statistics in arithmetical quantum chaos are much more similar to those of classically 
integrable systems than to the ones of generic classically chaotic systems. 

Sections 3.6 and 4.5 on fluctuations in geodesic length spectra and on convergence prop- 
erties of the Selberg zeta function, respectively, had to take the different spectral statistics 
for the hyperbolic Laplacian on arithmetic and non-arithmetic surfaces into account. It 
turned out that in both sections an application of inverse quantum chaology proved useful. 
The desired results, however, could only be obtained in the non-arithmetical case. There 
Berry's observation on the saturation value of the spectral rigidity, Aoo(-E') ~ ^logE, 
E oo, sufficed as an input to apply Hejhal's theorems [51] of inverse quantum chaology. 
In section 3.6 the remainder term to the leading asymptotics in the PGT thus followed 
to be of the form Qr{1) — e'^^ uj{l), with uj{l) denoting some unknown function containing 
powers and logarithms of I. The result of section 4.4 on A.^{E) for arithmetical systems 
could only bring down the upper bound of e'i^ for Qr{1) to e^^uj{l). Numerical evi- 
dence obtained from three arithmetic groups, however, suggested that the exponent for the 
remainder term Qr{1) should also be | for arithmetical systems. Thus the inapplicability of 
the inverse quantum chaology reasoning for arithmetic groups rather seems to be of a tech- 
nical nature than of a fundamental one. Regarding convergence properties of the Dirichlet 
series for the Selberg zeta function, however, the difference between the arithmetic and the 
non-arithmetic case seems to be not void of consequences. Using Hejhal's Theorem 10.10 
[51] of inverse quantum chaology and Berry's result on Aoo(-E) for generic systems revealed 
a convergence of the Dirichlet series at least until directly before the critical fine, which 
is the physically interesting domain. Again, the strong spectral fluctuations present for 
arithmetical systems prevented an application of this method to the latter. The statistical 
model for the convergence properties that was introduced in [6] now hints at the reason for 
this obstruction. The exponentially increasing multiplicities of geodesic lengths yield an 
exponential growth of the multiplicities of pseudo-lengths. The exponent a describing the 
latter increase is a measure for the distance of the domain of convergence to the critical 
line. 

In conclusion, the arithmetical systems that were studied in the present text excel by 
properties of important classical and quantum mechanical quantities that distinguish them 
from those of strongly chaotic systems that are commonly considered as generic. Rather 
convincing heuristic reasons for the exceptional spectral statistics could be derived from 
the classical properties of arithmetical chaos. It was mainly a combination of heuristic 
reasoning with the intuition gained from numerical observations that could be used in con- 
junction with rigorous results. This amalgam of different methods proved to be particularly 
fruitful. The arithmetical systems thus turned out to provide a convenient test-ground for 
the ideas and methods developed in the framework of periodic-orbit theory. 
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A The Riemann Zeta Function 

In the main body of the present work the Selberg trace formula and the Selberg zeta 
function have been extensively used. Historically, Selberg introduced his formalism in 
close analogy to the theory of the Riemann zeta function and the distribution of prime 
numbers. Since many techniques appearing in the latter theory can be carried over to the 
case of the Selberg zeta function, some important tools that were developed to study the 
Riemann zeta function will be briefly introduced in this appendix. Those details that will 
be omitted can be found in e.g. [97, 59]. 

The Riemann zeta function ((s) is a meromorphic function for all s G (D that has a 
simple pole at s = 1 with residue one. Por i?e s > 1 it is defined by 

OO -1 

c(«) = Erj=n(i-p-r^ (A.1) 

n=l "' P 

The Dirichlet series for ^(s) extends over all positive integers n, whereas its Euler product 
runs over all primes p. The analogy to Selberg's zeta function Z{s) leads to an identification 
of the primes p with e^^^\ where /(7) denotes the length of the primitive closed geodesic 
related to 7 e F, and an identification of the integers n with the pseudo-orbits p. These 
identifications describe what for the Selberg zeta function is the "classical" side of its 
theory. The "quantum" side is missing in the theory of C,{s) in that there is no self-adjoint 
operator known, whose eigenvalues are related to the non-trivial zeros of the Riemann 
zeta function. The trivial zeros of ai'c, however, explicitly known to be located at 
Sk = —2k, k e IN. If a self-adjoint operator related to the non-trivial zeros were known, 
the Riemann hypothesis (RH) would be true, since then (depending on how the operator 
is defined) either 7^ or 7^ is a real eigenvalue of it. This means that the non-trivial zeros 
s„ = I ± ijn lie on the critical line Res = |. But, the RH still being unproven, the 7„ 
can be complex. It is only known that < i?es„ < 1. There is, however, tremendous 
evidence in favour of the RH from extensive numerical computations of non-trivial zeros 
[74]. Therefore it seems to be justified to assume the vahdity of the RH throughout, if not 
stated otherwise. 

The importance of the Riemann zeta function derives not only from its connection to 
the RH, being one of the most famous unsolved problems in mathematics, but also from its 
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decisive role played in the proof of the prime number theorem (PNT) and for estimating the 
remainder term appearing in the PNT. Thus the analytic properties of C(s) are essential for 
describing the distribution of prime numbers, which is a central issue of number theory. It 
will now be explained how the validity of the RH influences the magnitude of the remainder 
in the PNT. Let therefore 

7r(x) := p<x} (A.2) 
be the counting function of prime numbers. The PNT now states that 

7r(x) = li{x) + Q{x) , Q{x) — o{x/ logx) , x — > oo . (A. 3) 

Tr{x) may also be expressed by two other functions that have been introduced by Chebyshev. 
Using the von Mangoldt function 

{log p , Ti ^ p^ 
, otherwise i-^-^) 

these are defined as 

n<x k>lp''<x 

= E^(^^)- 
fc>i 

When < 2 is reached the last series brcakes off, i.e. when k > Writing ^(x) = 

9{x) + 7l{x) one can easily estimate that TZ{x) = 0(y^(logx)^), x ^ oo. From the 
definition of 9{x) one obtains 

de{x) = ^ logp 5(x -p)dx , (A.6) 
p 

and thus 

r^-i: ^^sp r '^(^ - = E 1 = ^(^) ■ (A-7) 

J2 logt ^ J2 logt ^ 

An integration by parts yields 

T^{x)= T^ = T^+ dt n \o +0(1) . A.8 
^ ^ J2 logt logx J2 t{logt)^ ^ ^ ^ ^ 

Using 9{x) — ijj{x) + 0{^/x{logx)^) leads to 

nix) = p^+ Tdt + OiV^logx) . (A.9) 

logx J2 t{logt)^ 
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Proving the PNT is thus equivalent to determining the leading asymptotic behaviour of 

the Chcbyshcv function ip{x) for x oo, and estimating the remainder term Q{x) can be 
achieved by knowing the remainder to the asymptotics of ip{x). 

At this stage now the Riemann zeta function enters the game. Using its Euler product 
for Res > 1 one observes that 

Employing the Cauchy integral theorem one can derive that (for 6>l,a>0, Ot^I) 

J_ r^+''^ ^ f 1 , a > 1 

2m Jb-ioo *s\0,a<l' K ■ ) 

This result may be used to show that {x ^ IN) 

,/ X w X 1 f'>+"^dsfxY 1 fb+ioo x'C'(s) , 

= E A(") ^ - (-) = 7^ . (A.12) 

a relation that together with (A. 9) clearly shows how the PNT is related to the analytic 
properties of One can now use the WeierstraB representation of C,{s) as a product 

over its zeros to obtain from (A.12) the explicit formula of Riemann and von Mangoldt, 

i^ix) = ^ - E ^ - ^ log(l - X-') - ^ , (A.13) 

where the (conditionally convergent) sum runs over all non-trivial zeros Sn — + ijn 
of C(s)- Denoting (Tq := sup s„ = Pn + iln}, the RH is equivalent to ctq = |. An 

estimate for the sum over the non-trivial zeros in (A.13) can be found in [59], J2s„ Y~ ~ 
0(x'^"(loga;)^), leading to 



V'(x) = x + P{x) , P{x) = 0{x^°{logxf) . (A.14) 

X 



Since tt-^ = li{x) — + 0(1), one observes from (A. 9), using (A.14), that 



7r{x) = li{x) + Oix"" logx) , (A.15) 

which is the PNT with an estimate for the remainder term Q{x). This relation shows the 
influence of the validity of the RH on the PNT. Notice that in order that Q{x) — o{x/ log x) 
one has to show that o"o < 1, i-e. C(s) must not have a zero on Res — 1, nor a subsequence 
of zeros with real parts accumulating at 1. To show this was the main achievement of 
Hadamard and de la Vallee Poussin in proving the PNT in 1896. 

The point of view that will be taken now introduces a connection of the RH (equivalently 
the PNT) to the convergence properties of a certain Dirichlet series. The general theory of 
Dirichlet series is presented in [48] , where one can find the following results. A generalized 
Dirichlet series is a series of the form F{s) — J2'^=i ctne"^"*, a„ e (T, Ai < A2 < A3 < 
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If A„ = logn, it is called an ordinary Dirichlet series. In any case there exists a number 
(Tc G iRU{±oo} such that the Dirichlet series converges for Res > Uc and diverges for 
Re s < Oc- Since the series in addition converges uniformly on compact sets, Fi^s) is 
a holomorphic function in the domain of convergence. There also exists a number (Ja, 
Ca > c"c, such that the Dirichlet series converges absolutely for Res > Oa- The abscissae 
of convergence Uf. and a a are determined by 



(7c — lim sup - — log 



N 
n=l 



lim sup - — log ^ |a„ 



(A.16) 



Applying (A.16) to the (ordinary) Dirichlet series for (^(s) one obtains (a„ = 1, n G IV) 
that (Jc = 1 = (Jq. The location of the abscissae of convergence follows from the existence 
of the pole of C(s) at s = 1, which prevents the series to converge for i?e s < 1. 

Our goal now is to define a Dirichlet series that yields a meromorphic function for s G (D 
and whose abscissa of convergence is given by o^. i.e. by the non-trivial zero of (^(s) 
with largest real part. In the vicinity of s = 1 the Riemann zeta function behaves like 
C(s) = +7 + *^(("5~ 1))' 7 being the Euler constant. Thus = — ;^ + regular terms, 
s — > 1, is a meromorphic function with poles at s = 1 and at the zeros of C(s)- Then 



/(«) C(«) + 



as) 



(A.17) 



defines a meromorphic function that is holomorphic for Res > a^. In the critical strip 
< i?e s < 1 its poles are located at the non-trivial zeros s„ of C{s). Inserting for Res > 1 
the Dirichlet series' for ({s) and yields 



On 



= 1 - A(n) . 



By (A.16) the abscissae of convergence for this Dirichlet series are 

1 



N 



Since A(n) ^ only for n 
one concludes that Y,n<N 1 1 



lim sup 

N-*oo 

lim sup 
lim sup 
lim sup 

N-^oo 



logN 
1 

logTV 
1 



logE|l-AC 



n] 



n=l 

TV 



log 



E(i-AN) 



n=l 



logiV 
log|P(iV)| 



log I AT -^(A^) I 



logA^ 

p'^, and n{x) = ij^{p^ 
A(n)| ~ 2N for N 



(A.18) 



(A.19) 



oo. Therefore the Dirichlet series for 
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f{s) converges absolutely for Re s > aa — 1. The interesting observation one makes with 
(A. 19) is that the abscissa of conditional convergence is determined by the first pole that 
is encountered when moving with the axis Res = const, to the left. Namely, since cTc has 
obviously to fulfill ctc > ctq, but the asymptotics P{x) = 0(a;'^"(loga;)^) gives cTc < ctq, one 
concludes ac — o-Q. The Dirichlet series for f{s) therefore converges in the maximal possible 
domain. In sloppy terms one could call this Dirichlet series a "detector" for the non-trivial 
zero of with largest real part, indicating through its convergence properties. 
The condition Uc > ctq now also gives a lower bound for P{x), namely 

P(x) = n{x''°-') V£ > . (A.20) 

Thus the true magnitude of the remainder term Q{x) in the PNT is asymptotically bounded 
from below by x'^°~^ for all £ > and from above by x'^^' log a;. Therefore Q{x) = x'^^' -u^lx), 
where u!{x) is some combination of logarithmic functions. It is thus the "leading" term x'^° 
that determines the fine structure in the PNT. If now the RH were true, the remainder 
term in the PNT would have the "minimal" asymptotic behaviour oc -^/x • u}{x). In 

this case the best lower bound available is a;(x) = [59]. 

B Desymmetrizing the Hyperelliptic Involution 

The general procedure of desymmetrizing the quantum problem of a particle on a hyper- 
bolic surface possessing symmetries has been reviewed in section 4.1. This appendix now 
contains an explicit application of the general formulation derived in [100, 101] to a rather 
simple case, namely the so-called hyperelliptic involution, emphasizing the point of view 
employing the Selberg trace formula. A discussion of the example used here can be found 
in [5]. The following presentation, however, differs a little from the one given in [5] by 
being closer to [100, 101] in order to serve more explicitly as an example for the general 
situation. 

The symmetry under consideration is the hyperelliptic involution (see e.g. [40]) present 
on all hyperelliptic (compact) Riemann surfaces. The latter ones can be realized as two- 
sheeted coverings of the sphere. If the surface r\7^ is hyperelliptic and of genus g, this 
covering is branched at the 2g + 2 WeierstraB points. The operation that interchanges the 
two sheets of the covering is an involution (i.e. the symmetry group is isomorphic to ^2 if 
no other symmetries are present, as will be assumed henceforth), called the hyperelliptic 
involution. The hyperelliptic surfaces form a (2gf—l) -dimensional subvariety of the {3g—3)- 
dimensional moduli space of compact Riemann surfaces of genus g; for g = 2 all compact 
surfaces are hyperelliptic. 

For the following it proves useful to change the model of hyperbolic geometry and 
to pass to the Poincare unit-disc V = {w E <C; \w\ < 1} by mapping the upper half- 
plane H via z I— > w = Cz, C = ^ (Ii)' z E H. SL{2,1R) is then mapped to 
SU{1, 1) = {(if) ; - = 1} by conjugation with C. kg e SU{1, 1) operates on 
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T> via fractional linear transformations of the form 



aw + (3 



w eV . 



(B.l) 



gw = 



l3w + a ' 



The conjugation of a Fuchsian group V C SL{2, M) by C yields a discrete subgroup of 
SU{1, 1) that will also be denoted as F by abuse of notation. The geodesies of the appro- 
priately transformed hyperbolic metric are the half-circles and straight lines perpendicular 
to dV = {w E<C; \w\ = 1}. [17] may serve as a reference to find further details concerning 
the geometry of this model. 

Any Fuchsian group F leading to a hyperelliptic surface r\D may be obtained in the 
following manner (see also [12]). One can choose a fundamental domain !F d T> whose 
boundary dT consists of A:g geodesic segments. The corner points wi, . . . , w^g are enumer- 
ated in ascending order when going counterclockwise along dJ^ and starting with wi on 
the positive real axis. • • • ^W2g are placed in the upper half of T> [Im Wi > 0). The 
remaining corner points are obtained as W2g+i = —Wi, . . . ,W4g = —W2g- T is therefore 
symmetric under the operation w i— > —w. Having fixed W2, • • • ,tt'2g one has to vary w\ 
(hence also ^2^+1 = on the real axis until the constraint area(J^) = ^T^{g — 1) is 

fulfilled. The fractional linear transformations of the form (B.l) that identify opposite 
geodesic segments comprising dT then serve as generators for a strictly hyperbolic Fuch- 
sian group F possessing as a fundamental domain. This identification of pairs of edges 
of T is compatible with the symmetry w ^ —w, so that F\P is a compact hyperelliptic 
surface of genus g. The above construction yields every such surface by choosing the corner 
points W2, ■ ■ ■ , W2g appropriately in the upper half of T>, therefore clearly showing that the 
subvariety of hyperelliptic surfaces in the moduli space of compact surfaces of genus g is 
of (complex) dimension 2g — 1. 

The generators 61 , . . . , will be enumerated such that bi identifies the geodesic segment 
connecting Wi and Wi+i with the one connecting —Wi and — Wj+i, i — 1, . . . , 2g. By their 
very construction then 6i+2p = hj^ for i = 1, . . . , 2g. The generators obey the constraint 



The hyperelliptic involution is realized as the mapping w —w, which can be represented 
by the matrix S = i) ^ SU{1, 1), S ^ T. Since the 6j's identify opposite edges of JF, 

one notices that SbiS = 6~\ where 5"^ = 1 as an identity in PSU{1, 1) = SU{1, 1)/{±1} 
has been used. The 2g + 2 fixed points of S on T\D, i.e. the Weierstrafi points, can now be 
identified. Obviously, i^i = is fixed by S. The second WeierstraB point U2 is represented 
by the corner points of J- that are all being identified to one another by F. The remaining 
fixed points M3, . . . , ^23+2 are given by the mid-points of the 2g pairs of opposite edges of 
as determined by the hyperbolic metric. 

A new Fuchsian group F' is now introduced by adjoining S to F, i.e. F' consists of 
all words in the generators 61, . . . , 643, 5", subject to the constraint (B.2) and fulfilling 
b^^ = &i+2g- = Ij SbiS = bj^. It is thus possible to rewrite any word in those generators 
as being of the form S^bi^ ■ ■ - h^, e e {0, 1}. One therefore obtains that F is a normal 



6162 63^4 



,-1 



• • • b2g-lb2g &l • • • &2g-3&2s-2&25-1^29 = 1 • 



(B.2) 
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subgroup of index two in F'. The latter decomposes disjointly according to 

r' = r u 5r , (b.3) 

thus reproducing for the hyperelhptic involution S the algebraic setting of symmetries 
reviewed in section 3.5. The symmetry group S is yielded as S = {1,5'}, S*^ = 1. Thus 
S = and also for the unitary dual E* = Explicitly, S* is given by the two 

representations x+ X- of T'; x+ denotes the trivial representation x+il') — 1 foi" 
7' e r', whereas x_ is defined as 

. (B.4) 

Since S e SU{1,1) is elliptic dtrS"! < 2), T' is not strictly hyperbolic. Fortunately, 
it is possible to identify all elliptic conjugacy classes of V explicitly by their fixed points 
on the surface r'\D. An elliptic R E V has one fixed point zr in the interior of T>, and 
its conjugacy class {-R}r' fixes the set of points T'zr that are identified under V. To each 
elliptic conjugacy class there hence corresponds the point T'zpt on the surface T'\D; but 
this must be one of the WeierstraB points ui, . . . ,U2g+2- The Fuchsian group F' therefore 
contains 2g + 2 eUiptic conjugacy classes {Rjv, all of them of order m{R) = 2 (meaning 
the minimal positive integer with = ±1). One can easily determine representatives 

for the elliptic classes, 

• Ui is fixed by S, 

• U2 is fixed by 5*6162 ^&3&4^ • • • ^29-162^, 

• Ui is fixed by Sb^_}2, ior i — 3, . . . ,2g + 2. 

All other conjugacy classes in F' are hyperbolic ones. 

The general receipt of [101] now proceeds in constructing the unitary representation of 
F' that is received as being induced from the trivial representation of its subgroup F. To 
this end define the one dimensional representation (on (D) 

«^') - { J ; y ^ sr • 

The induced representation p : F' — > End((D ®(D) is then obtained as 

p{^'){v(Bw) := [xb')v + x{l'S)w] © [x{S^')v + x{Sj'S)w] , (B.6) 

for 7' G F' and v,w E^. On F this representation operates trivially, p{'y){v (Bw) = v (Bw, 
7 e F; on ^F, however, it acts according to p{S'y){v © w) = w © v, 7 e F. As a matrix 
representation on (D^ the induced representation is hence given by 
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Later, in the Selberg trace formula, trp{'^')^, k e W, is needed. For 7 e F, clearly 
tr pi^Y — 2) but for 7' G ST one obtains that tr pi^'Y = 1 + ("l)*^- Altogether, one can 
reformulate this as tr pipi'Y = X+{.l'Y + X~{.l'Y f*^^ l' ^ ^' ■ 

Venkov and Zograf now demonstrate [101] that the hyperbolic and the elliptic contri- 
butions to the Selberg trace formulae for F endowed with the trivial representation and for 
F' endowed with the induced representation p coincide, respectively. The hyperbolic terms 
thus yield 

- /(7) g{kim ^ trp{ifl{i)g{kl{i)) 
4fe2sinh(|/(7)) ^^^,ti 2sinh(|/(7')) ^'^ 

^ l{i)g{kl{i)) ^ X-{if Ki) 9{kl{i)) 

m„ h 2 smh(|/(7')) {,Vr, h 2 sinh(|/(y )) ' 

where the outer sums extend over the respective primitive hyperbolic conjugacy classes. 
This result yields a factorization formula for the Selberg zeta function, 

Zr{s) ^ Z+{s) ■ Z:^,{s) , (B.9) 

where for Re s > 1 

zMs) := n n (1 - ^-^^^"^'^^'0 , 

{7'}j„ n=0 

00 

ZAs) := n n(l-X-(7)e^^^^"^'^^'0 , (B-IO) 

{7'}r' "=0 

which should be compared with (4.2) and (4.3). The factorization (B.IO) comprises the 
desymmetrization with respect to the hyperelliptic involution on the level of the Selberg 
zeta function because the non-trivial zeros of Zyi{s) correspond to the eigenvalues of —A 
referring to the two symmetry classes of S. The reason for this is that the respective 
wavefunctions transform under x± according to '4^±{'^'z) = X±(7')V^±(-2), 7' £ F'. Hence 
the ip±{z) are invariant under F (x±(7) = 1 for 7 e F), and transform under S as ip±{Sz) — 

Finally it should be remarked that the exphcit knowledge of the eUiptic conjugacy 
classes of F' allows to determine the elliptic contributions to the Selberg trace formulae for 
F' endowed with the representations x± ^ 

_^g + l h{p) 



7-00 cosh 7rp 



cosh(7rp) 

see [5]. There one can also find a numerical evaluation of the Selberg zeta functions Zp,{s) 
on the critical line for a specific surface of genus g — 2. 
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C Estimates of Remainder Terms 



In analytic number theory several estimates are used to describe remainder terms to the 
leading asymptotic behaviour of functions of interest, the most prominent example being 
the remainder Q{x) in the PNT, see appendix A. Since in other fields, like in physics, 
some of these estimates are not so commonly used, their definitions will be supplied in this 
appendix. 

Let therefore f{x) be a function of the real variable x > that shall be estimated for 
X — > oo and compared to the positive and monotonic function g{x). 
The first estimate is 

I f(x)\ 

f{x) = 0{g{x)) :^ limsup ^^f^ < oo . (C.l) 

x—^oo g\x) 

Asymptotically |/(a;)| is thus bounded by g{x), which might therefore also be referred to 
as an upper bound. 

Another, stronger, upper bound is 

f{x)^o{g{x)) :^ limsup^^^^O. (C.2) 

A sort of asymptotic lower bounds is provided by Jl-estimates, which are given as 

f{x) = n{g{x)) :^ limsup^^^>0, 

x-^oo g[x) 

f{x) = n+{g{x)) :^ limsup44>0' (^-3) 

»oo g\x) 

fix) = ^-igix)) :^ liminf 44 < ■ 

The O-estimate can also be formulated slightly differently: f{x) — 0{g{x)), if there 
exist some Xq > Q and M > such that < M g{x) for x > Xq. Analogously, 

fix) = n{g{x)), if \f{x)\ > Mg{x) for x > Xq, etc. 
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Figure 2: The quantity Q'{1) = \Qr{1) \ / e 2 is shown for a) Artin's biUiard, b) the regular 
octagon and c) the Gutzwiller octagon. 
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Figure 3: The model Pm(s) applied to the odd symmetry class of Artin's billiard is shown 
as the full curves compared to the actual level spacings distributions in the intervals a) 
< p < 100, b) 100 < p < 200, c) 250 < p < 300, and d) 500 < p < 510. The dotted 
curves show Poissonian distributions. 



94 




Figure 4: The same as in fig. 3, but for the even symmetry class. 
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Figure 5: The number variances obtained from the eigenvalues are given as the full curves in 
comparison with the model, which is shown as the dotted curves. The Poissonian (dashed 
curves) and GOE (dashed-dotted curves) results are added, a) and c) refer to the interval 
164.92 <p< 229.70, b) and d), however, to 250.86 <p< 296.83. 
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Figure 6: The saturation value ^^M^ooip) the model number variance for the odd sym- 
metry class of Artin's billiard is shown as the full curve. The dots mark iS^(p) obtained 
from the eigenvalues, see [91]. 
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